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A MATHEMATICAL DERIVATION OF THE EXCHANGE OF 
A LABELED SUBSTANCE BETWEEN A LIQUID FLOWING 
IN A VESSEL AND AN EXTERNAL COMPARTMENT* 


W. C. SANGREN AND C. W. SHEPPARD 


MATHEMATICS PANEL AND BroLocy Division, OAK RmMGE NATIONAL LABORATORY 
Oak RipGE, TENNESSEE 


A situation is considered in which a fluid containing a substance flows through a vessel at 
a constant rate, the substance being permeable to the vessel wall. In the region outside the 
vessel there is supposed to be rapid mixing in the direction perpendicular to the axis of the 
vessel but no mixing longitudinally. The solution for the spatial distribution at any time is 
given for the case of an arbitrary initial distribution along the vessel length in the absence of 
an input. The solution is also given for the case of a single impulsive input, the concentration 
being initially zero everywhere. 


In isotope studies, substances are often found to disappear rapidly from 
the circulation (Sheppard ef al., 1951) or rapidly move back and forth 
between vascular and extravascular compartments (Walker and Wilde, 
1952). In such cases the approximations involving labeled substances 
leaving uniformly mixed pools cannot be considered adequate (Sheppard 
and Householder, 1951). 

The general problem of the nonuniformly mixed circulatory pool is a 
complex one since it requires a detailed consideration of the dispersing 
action of the variable circulatory paths, the process by which the label 
moves within the blood vessels, the mechanism by which it is removed 
from the circulation, and other matters whose details are considered in 
part elsewhere (Sheppard et al., 1953). One problem remains, however, 
which will be of increasing interest to isotope workers. This is the descrip- 
tion within a vessel of the transient behavior of an initially intravascular 
isotope which exchanges with an initially non-labeled external compart- 
ment. This problem is of possible broader physiological interest since its 
application to exchange across the wall of the renal tubule or of the in- 
testine awaits future consideration. In general, then, given an initial 
distribution of activity along the vessel, how will this distribution alter 
with time, what will be the temporal and longitudinal spatial relations of 


- * Work performed under Contract W-7405-eng-26 for the Atomic Energy Commission. 
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its appearance in the external compartment, and how will these depend 
on the system parameters such as flow velocity and exchange rate across 
the wall? We shall not consider here how this initial distribution is to be 
measured. However, if a special initial condition is given by a sharp pulse 
of given area, the pulse would represent the sudden injection of a small 
slug of labeled material, a situation frequently encountered. 

Because of the analogy between the interfusion of an isotope and the 
transport of heat, related problems are found in the theory of heat ex- 
changers and regenerators (Carslaw and Jaeger, 1947; Goldstein, 1938). 
However in the present problem the boundary conditions and resulting 
solutions are different. It should be noted, however, that as in the previous 
cases certain restrictive assumptions must be made through mathematical 
necessity. These will be discussed at the concluding point of the paper. 


p 
A G(x,O) =O | 
p F (x, O)= a (x) —>v 


ceonreranmmenens aa spy aie eo x 


FIGURE 1 


The problem under consideration is illustrated in Figure 1. 
The following notation is used: 


x = the distance along the vessel from some fixed reference point. 
V = the constant velocity of flow inside the vessel. 
= time. 
G(x, t) = specific activity in the external compartment. 
F(x, t) = specific activity in the internal compartment or vessel. 
p’ = fraction of contents of external compartment exchanging per unit 
time. 
p = fraction of contents of internal compartment exchanging per unit 
time. 


In mathematical language the problem is given by 


G, = p' (F—G), (1) 
G(x, 0) =0, (2) 
F,= p(G—F) — VF., (3) 


F(x, 0) =a(x), (4) 
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where the subscripts denote the partial derivative with respect to that 
subscript and where one assumes no source or the equivalent condition of 
a bounded solution for all x. Further let us assume that the given initial 
distribution a(~) be continuous and vanish when x < 0. The quantities 
p, p', and V are assumed to be constants. 

The problem may be conveniently solved by using the Laplace trans- 
form (Churchill, 1944). Applying the I aplace transformation formally to 
equations (1) and (3) and using equations (2) and (4) one obtains 


Of 


sg=p'(f—g) and as = Ole p (Ge 


where 
g(x, s) =L{G(x, 1} = [eG (x, 1) di and -f(#)ss) =I ROMO. 
0 


Eliminating g between these two equations gives 


OF a. ae Ss _ pp 
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x 


pl 
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The solution of this first order linear equation is 


f(x, s) =e-# [fe erdst+ A], 


where A is an arbitrary constant. Since F(x, ¢) and therefore also f(x, s) 
are to be bounded when x tends toeither + ~ or — ~, and since we may 
assume that Re(c) > 0; i.e., the real part of the complex parameter s is 
a large positive number, it follows that 


—o a(z) oe 
A = — if is é az , 
Now a(z) = 0 for z < 0 and therefore 
= a(z) 
— pcr ced gz, (5) 
j(%, s) =e I 7 eds 


This transform must now be inverted in order to obtain F(z, t). In general 
the inversion process may be difficult. However the following rather gen- 
eral case indicates the character of F(s, ?). 
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Assume that a(x) = 0 when x > L > 0. By invoking the second law of 
the mean for the integral calculus (Townsend, 1928) equation (5) becomes 


¥ (45.8) =90 when nce ee 


. pg (1s fie See ies 
f(x) =e vf eda = : 


whereO<#<2<L; 


ee A) ee ee RS Pe et ae 
HED) = G V fe 43 = V : . 


whereO <4 <L <x. 

It is clear that depends upon x while x does not. 

The inversion can now be carried out with the aid of standard Laplace 
transform inversion formulas (Churchill, 1944). If L-1{k(s)} = H(2), i.e., 
L{H(t)} = h(s), then 


oes aes + P e— @tr)t 


VC + panache 
since 
POS s+ p’ coke P: + 1 p 
VGO* sts-Fipt se) spp Sp eee 
Further 
1 x x 
1) pcr) = p—(2/V) 25 ES = a 
Ls ee e rH (x,t >) when 0<7<!, 
= 0 when 0<i<Z, 
where 


, 
H(x, t) = Aes er m/v yti/iet ei b 
pt+p s 
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IL {h(s+a)} = e-*L-{h(s) } = e-*H (2), 


ifn f= f'n0—nar, 
and 
Lf err = [)(2 Vht) , 
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a modified Bessel function, it follows that 


[ e- PZ (av zee: ') 
a fe ener) Ty (avEe@ es) —+) at] 
+ a 7 [ ener ChE pf en wterre-wtmn 
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The specific activity in the internal compartment, therefore, is given by 


Boge}. A) CREA G® 
7 p [Ese Ea 
F(x, 1) =0(2)[P +P wrent], (0<:<4), i 
(Ole <a <a) 
F (x,t) = a(a)[—, fotos eter | in 


x x 2 
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Fiz, i) =0, (0<1< ree Oe (10) 


x—L 
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(11) 
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x—L 


F(x, t= aa) ee 2/6 a —L,t- 


nT: (x) e-6/"WH( x, 


et), (12) 


(QO 2 = bere 


If desired, the specific activity in the external compartment can be ob- 
tained in a straightforward manner by using the above results and equa- 
tion (3). It should be kept in mind that a() is a function of x. 

The study of the above solution is quite informative. The six ranges for 
the solution are, upon reflection, quite logical and, incidentally, indicate 
that the Laplace transform method is very appropriate to this type of a 
problem. The appearance of the variable combination ¢ — «/V indicates 
wave propagation in the positive direction. This is not surprising consider- 
ing that the elimination of G between equations (1) and (3) leads to an 
equation of hyperbolic type for F. Also from the physical point of view the 
fluid should carry a propagated disturbance with it. Other similar re- 
marks can be obtained by the reader. 

A special case of some interest occurs when a(x) is a square wave, i.e., 
a(x) = K for 0 < «x < Land a(x) = 0 otherwise. In this case a(#) = 
and a(x) = K. Although it is somewhat tedious, it is not difficult to show 
in this case that the solution obtained above formally satisfies the bound- 
ary value problem given by equations (1) through (4). In this situation 
the rectangular distribution initially has a vertical leading and a vertical 
trailing edge. The various portions of the solution may be recognized as 
follows: 


(a) Solution for negative values of « (zero). 


(b) Solution at a point initially inside pulse before the trailing edge 
arrives. 


(c) Solution at a point initially inside pulse after the trailing edge 
arrives. 


(d) Solution at a point initially beyond the pulse and before its arrival 
(zero). 


(e) Solution following the arrival of the leading edge. 
(f) Solution after the passage of the trailing edge. 


If one takes the limiting case of an infinitely narrow and high pulse of 
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finite area M the important contributions are those of (11) and (12). 
These become 


F (x, t) =Me-(/V) 


F (x, t) = Me e-lo#/¥)—ett-@/V 1] (2 ( ae -+)) 
p ; \ Pea ap 


px aif aangisie —1/2 
«(fe y"C-F) 


These expressions are simplified by adopting the natural coordinates 
p = px/V and gq = p’(¢ — x/V). Equation (13) describes the propagation 
of the pulse downstream with uniform velocity and progressive attenua- 
tion as the label moves out of the internal compartment. Equation (14) 
describes a progressively growing tail which follows the pulse and declines 
monotonically when observed as a function of time at a fixed point. This 
part indicates the return of activity from the external compartment. 
The natural coordinates # and g might have been used from the start. 
With these coordinates equations (1) and (3) assume the simple form 


F,=G-F and G,=F-G. 


Fae 
—_ 
ies) 

a 


(14) 


However, the solution of the problem in this form appears less direct since 
the Laplace transform techniques cannot be used as conveniently. 

We recognize that the problem as stated here cannot be more than an 
imperfect first approximation to the description of such physiologically 
important phenomena as the exchange of isotopes across the walls of renal 
tubules, capillaries, or intestinal mucosa. The assumption of uniform 
velocity cannot be reconciled with the variations in vessel caliber and 
ignores the known existence of laminar flow in the circulation and pos- 
sibly elsewhere. The assumption of uniform mixing and absence of longi- 
tudinal movement in the external compartment also is unrealistic. A 
realistic treatment of longitudinal movement in the external compart- 
ment appears to be quite difficult. Nevertheless the problem which we 
have discussed here would appear to be the best point of departure in 
considering these additional complexities. 

If p, p’, and V are considered to be functions of x, equation (5) becomes 

—(*%dz rz — (*edy 
fH) = eo! (“ay gS dz 
In particular cases it may be possible to carry out the integrations and 
invert the transform. 
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ON THE PROPAGATION OF A DISTURBANCE THROUGH A 
VISCOUS LIQUID FLOWING IN A DISTENSIBLE 
TUBE OF APPRECIABLE MASS 


ROBERT B. Jacons* 


NATIONAL BUREAU OF STANDARDS 
BOULDER, COLORADO 


The velocity of propagation of a disturbance wave in a liquid flowing in a distensible tube 
is computed. The mathematical model is more general than those used in previous analyses: 
the tube wall properties are realistic; the convective part of the axial inertia forces is taken into 
account; radial inertia forces of both the fluid and tube wall are present; viscous stresses are 
present. 

Four parameters influencing the velocity of propagation are obtained and discussed. 
Curves are plotted illustrating the effects of the parameters. Contrary to the results of previous 
analyses, viscous effects are shown to be appreciable in blood flow. It is also shown that radial 
inertia effects can be important in laboratory set-ups. 


Introduction. The velocity of propagation of a disturbance wave through 
blood vessels has been of interest to biological scientists for at least 150 
years. Thomas Young (1808) did some work on the problem while he was 
preparing some lectures on blood circulation. He derived an expression for 
the velocity of propagation of a very small disturbance through an in- 
viscid liquid, which was at rest, in an elastic tube; he neglected the radial 
inertia forces of both the tube wall and liquid, and considered a tube wall 
which is unreal. 

A. I. Moens (1878) rederived the same formula, not knowing of Young’s 
work. Since then, the formula has been identified with Moens. Although 
verification through empirical data has not been realized, Moens’s formula 
is still used. The velocities of propagation predicted by Moens’s formula 
are always too large. In some experiments a correction factor of 0.6 has 
been required while in others, the deviation between theory and experi- 
ment is small. 

Since Moens, several investigators have been interested in the problem 
and have attempted to account for the effects of viscosity on the velocity 


* The material presented in this paper was adapted from the Ph.D. thesis written by the 
author at Harvard University. 
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of propagation. The shortcomings of these investigations are due to 
three factors: 

1) In linearizing the equations of motion, it was assumed that the 
inertia forces were small. This assumption is unfortunate because it 
eliminates the effects of the mean fluid velocity from the model. It will 
be seen that these effects, which are two in number, are very important. 

2). The effects of the radial inertia forces due to the masses of both the 
liquid and tube wall were neglected. It shall be seen that these effects can 
be important also. 

3) All previous work has been done on models which are unreal. One 
model [as an example, refer to Karreman (1952)] has tube walls which 
exhibit a constant modulus of elasticity and a constant wall thickness. For 
an elastomeric substance, these two specifications are mutually incom- 
patible. The second model [as an example, refer to Frank (1926)] has a 
tube of constant specific dilatability. The difficulty with both models 
is that they do not represent any actual system, much less a blood vessel. 


In the present work a generalized relation (elastic equation) is used 
which can be made to represent real systems accurately. The results of 
the present work will therefore be applicable to blood vessels, rubber 
tubes, and probably most actual systems. 

The physical model. The problem considered herein consists of an in- 
compressible Newtonian fluid of constant viscosity flowing through a 
distensible tube. All motions are such that symmetry always exists with 
respect to the axis of the tube, the tube axis is a straight line which is fixed 
with respect to the earth, and there are no tangential velocities in the 
fluid. 

The fluid velocities consist of a steady mean flow with an unsteady 
perturbation superposed upon it. It should be noted this model contains 
the convective part of the axial inertia forces which previous models have 
neglected. The tube wall itself is assumed to move in the radial direction 
only, and in such a manner that the slope of the tube wall is small com- 
pared with unity. This means that the radial velocity of the tube wall 
must be small relative to the velocity of propagation. 

The mass of the tube wall may be appreciable. The problem is formu- 
lated so that the mass of all tissues (not only the mass of the blood vessel 
wall) which affect the velocities of propagation is taken into account. 


THE MATHEMATICAL MODEL 
Nomenclature: 


a = radius ratio (=R/R>). 
a) = the mean radius ratio, a constant. 


a 
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“~r 


a quantity defined by equation (24). 

a quantity defined by equation (25). 

the elastic constant of the tube. It measures the “stiffness” of the tube. 
Bessel’s function of zero order. 

the unstressed length of the tube. 

the elastic exponent. It measures the variation of the elastic modulus 
of the tube. 

a “Mach” number (=qo/v»). 

fluid pressure. 

the mean fluid pressure, a constant. 

the mean fluid velocity, a constant. 

the axial component of velocity in the 7, ¢ plane. 

the radial component of velocity in the 7, ¢ plane. 

the radial coordinate in physical space. 

inside radius of the tube wall. 

the unstressed inside radius of the tube wall, a constant. 

time. 

a quantity defined by equation (23). 

axial component of velocity in the physical plane. 

radial component of velocity in the physical plane. 

the velocity of propagation. 

a generalized Moens’s velocity of propagation, called Moens Velocity: 


the viscosity factor, defined by equation (26). 

the axial coordinate in physical space. 

the frequency ratio which accounts for the effects of the frequency and 
the Moens velocity on the radial inertia forces (=w/Vo). 

a quantity which accounts for the mass which is effective in the radial 


inertia forces 
2 
ap) Dodie Ore | 
l ge. 0 ay 2 a é 


a time coefficient, defined in equation (20). 


B = a space coefticient, defined in equation (20). 
6 = the radius-length ratio of the unstressed tube (=Ro/L), a constant. 


oF 
v 


the eigenvalues which determine the modes of oscillation of a disturb- 
ance in the system. 

a dimensionless axial spatial coordinate (=x/Ro). 

the kinematic viscosity of the fluid, a constant. 


p = the fluid density, a constant. 
- p, = the “effective” mass of the tube wall material per unit length of tube, 


divided by 27. 
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¢ = a dimensionless radial spatial coordinate (=r/R). 
y = a stream function, defined by equations (17) and (18). 
w = the frequency in time. 

Equations. Three sets of equations are required to describe the mathe- 
matical model: the elastic equation, the equations of motion, and the 
continuity equation. 

The elastic equation expresses the relation between the radius of the 
tube, the pressure exerted on the tube by the fluid, and the properties of 
the tube wall. Under conditions for which the mass of the tube wall is 
negligible, the elastic equation is assumed to be 


ei ppt (as) 
a 


The values of the constants for four common tubes are given in Table 1. 


TABLE 1 
Tube C.(IN.H20) m 
Constant modulus of elasticity (E) and constant wall thickness (é)..| Et/Ro —2 
Constant speciiie dilatability (PD) asa = see eee 2/D —1 
Artery, data from Best,& Paylor'(1945) 232 eee ee 60 —1.632 
2” Penrose drainage tubing (data taken by the author)............ 645 —7.46 


The data for the artery and the Penrose Drainage Tubing are satisfied by 
(1) to closer than 2%. When the inertia effects of the tube wall (and sup- 
porting structure) become important the elastic equation is 


1 da Oa 
& chet aang ae meee (2) 


(It will be shown later that the wall is at ¢ = 1.) 

As the determination of the elastic constant (C.) and the elastic expo- 
nent (m) is usually performed statically, the effects of tube wall slope, on 
the elastic equation, have been neglected. Note, also, that the effect of 
the viscous drag of the liquid on the tube wall, with respect to the elastic 
properties of the tube, has been assumed negligible. 

The equations of motion are simply the Navier-Stokes equations sub- 
jected to the restrictions of constant density, constant viscosity, and axial 
symmetry 


ou 1 0p ou , 1 du , d% 
wt" 3 ate Pr (He 4t sat sa) (3) 
ov Ov dv 1 0p 1dv, d%v 9 af 
rats eB atte, sats at aa) (4) 
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Equations (3) and (4) are expressed in terms of the velocities and co- 
ordinates in physical space. It is now opportune to discuss the reasons for 
introducing the n, ¢ space. 

In order to specify the mathematical model uniquely, radial boundary 
conditions must be specified. Assuming that there is no slip between the 
fluid and the wall, the conditions in the physical plane are 


dR 


“(x,R, f) =0, v(x,R, MU iimteet 


kee 0;é)7=0, 0(%,:0, 4) = 0. 
As R(q, #) is unknown, the problem is unspecified in two ways: The loca- 
tion of the boundary (wall) is unknown, and a boundary condition 
(dR/dt) is unknown. 
The problem becomes unique if two transformations are used. Trans- 
posing from the physical plane to the n, ¢ plane by means of the relations 
x =Ron (5) 


and 

r =Ro¢ =RKoagd (6) 
fixes the tube wall at ¢ = 1. Note that when a partial derivative with 
respect to a new variable is indicated, the other new variables are the 
ones which are to be held constant. When the velocities are transformed 
through the relations 


UuU =Roqn (7) 


and 
v= R46 (5 +4 Lorem so) + a4], (8) 


the boundary conditions on the velocities become unique. If the differ- 
entiations are confined to particle paths, the velocities in the 7, ¢ plane 
may be written 


_dy 
Way 
and 
_d¢ 
eerie 


Note that the above derivatives (of one independent variable with respect 
to another independent variable) have a very special and limited meaning, 


and must be used cautiously. 
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The boundary conditions in the n, ¢ plane are 


Qn(n, 1,4) =0, qo(n, 1, ¢) =0, 


94 ss Oe pr=0n 
ad (n, 0, t) OF ge (n, ) ) 


There is no need to write the transformed equations of motion at this 
time..They are complicated; however, under the assumptions used in 
setting up the present model, they will be greatly simplified. 

Subjected to the restrictions which have been imposed upon the 
model, the continuity equation in the 7, ¢ plane is 


O9n da ,aa 
an | "3! 6 a6 

The major simplifications of the equations of motion are the results of 
four assumptions: 


1) The tube is “long”’; that is, if 


($4) +22 =0. (10) 


=— then Ole 


2) The radial velocities of the tube walls are small compared with the 
propagation velocities. This leads to 


52) = 068) and 9242 0(8).. 
én on 
3) The velocities are “nicely behaved” functions, as in laminar flow. 
4) The Reynold’s number, based upon the axial component of velocity, 
is restricted such that 
2Roa > 
ener ee 


1% 


If the above assumptions are employed in an order of magnitude 
analysis similar to that used by Prandtl (1934) in his boundary layer work, 


and if all terms of the order 5? and smaller are neglected, the desired simpli- 
fications are attained. 


The solutions to the problem may then be written 


m= IMs", = p> qa", 


n=0 


p= > pa, a= > aln)sn, 


n=0 n=0 


(11) 
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The coefficients gi, g, p™, and a™ are to be determined. They are 
not to be of larger order than the terms with the zero superscripts. The 
solutions are expressed in the form of (11) primarily because that form is 
compatible with, and follows from, the primary assumption that & < i. 
Notice that the series converge rapidly if the explicit assumptions are 
fulfilled; the errors may be made arbitrarily small by using a sufficient 
number of terms. Also, if 5 = 0, the original mathematical model is un- 
changed. In general, the calculations can be made by the method of small 
perturbations. 

In accordance with previous assumptions, the first terms in the series 
(11) express the steady mean flow; hence it will be assumed that 


g©) = q, = constant , gq =0, p©=p,=constant , a) = a) = constant . 


The assumption that g‘ is a constant is the only radical statement; be- 
cause of it, the mean flow does not satisfy the “‘no-slip” boundary condi- 
tion at the wall. As it is the average effect of the mean flow on the non- 
steady components which is important, the above assumptions are war- 
ranted. All the effects of the radial distribution of the mean velocity have 
thus been neglected. The convective effects, caused by the mean flow, 
remain in the model. 

The non-steady parts of the solutions are expressed by ¢, g/)), p"?, and 
a”), The complete solutions are therefore to be 


y= yt 690) 9, = 89), P= P,+ 62%, and a=a,+ 5a). (12) 


Notice that flow which is nearly ‘parallel’ in the 7, ¢ plane can be far 
from parallel in the physical plane. Refer to (8). This means that, although 
gs may be small, the radial velocity component in the physical plane (2) 
may be large. Without restricting the present analysis severely, it will 
therefore be assumed that the velocity component g is not larger than 
the ratio 5. Because the largest terms which were neglected in the series 
(11) are of the order 6’, the model is basically accurate to within 100 8%. 
Substituting (12) into the elastic equation, the equations of motion, and 
the continuity equation, and utilizing the aforementioned assumptions 
leads to the equations which define the mathematical model. 
(cee 7p icae ec aly fe es En ae) 
at an | pR? dn | Ril dg? '¢ 06 
do , 1 3p _4 (14) 


$a, ae pR? d¢ t 
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dp 1 da Pan (ae 
eee "9, So ee ee ; (15) 
Belt cow Las aPan t ° an 
a? (1) 
gts do OQn 
= (Y). 16 
reed re 5 ae (16) 


Solution. The solution to the system [(13), (14), (15), and (16)] is 
straightforward. The solution and all discussion will be confined to the 
n, ¢ space, which is dimensionless. The deviation of the dimension of any 
quantity from the corresponding dimension in physical space is thereby 
explained. In order to transform results to physical space (5), (6), (7), 
and (8) must be used. The pressure (p™) is eliminated from (13) by solv- 
ing for (p™) through (14) and (15). The term which represents the radial 
inertia force of the fluid is a function of radial position (¢). As the prob- 
lem is one-dimensional as far as the velocities of propagation are con- 
cerned, a radial inertia force averaged across the tube is used. The 
continuity equation (16) is satisfied by the stream function (yw) defined by 


oy 


77803 yes SE 
on 


qa —— [9 La ws]. (18) 


(17) 


and 


Substituting the expression for p™, (17), and (18) into the axial equation 
of motion (13) and simplifying leads to the single equation which describes 
the system 


oty Ory o*y _ 
2 Faq of 7 Sigg t Vo(l—M) os 
—v ( 0? ree Sah 
Sapte alt ae 9 3m . 


The mass which is effective in determining the radial inertia forces is 
represented by sz. It is defined by 


Pt 
ay 2 au) 
A generalization of the velocity of propagation obtained by Moensis Vo. 
It will be called “Moens velocity,” and is defined by 


Phe 


dp 
ar ACTS me Sat 
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Note that the Vo, calculated from the above expression, is in the n, ¢ 
plane and has the units: (unstressed radii/time). 

We call Mo a “Mach” number because its definition is similar to that 

of the Mach number defined in compressible flow theory. It is defined by 

M,=—.., 

0 Vs 

The solution to (19), for finite velocities, is made up of components of 

the form 


Y =constant e('+6n) J, (Ag), (20) 
where 
2 42 9 
> 2a B —a”— 2qaB+ Vo(1 — M2) B’ 
At = Ser OE (21) 
Reai * FMP) 


and Jo(A¢) is the Bessel function of zero order of the argument (A¢). 
The eigenvalue for the problem is A and its values are determined by 
the radial boundary conditions. If it is assumed that there is no slip at the 


Ficure 1 


tube wall, the values of A are merely the roots of the equation J, (A) = 0, 
which are shown in Figure 1. Note that the boundary conditions at the 
tube axis are also satisfied. The various eigenvalues determine the modes 
of motion of the system. These modes are shown in Figure 1. For the first 
mode (A = 2.40), the radial velocity distribution is essentially parabolic; 
the second mode (A = 5.52) has one reversal of flow, etc. 

As our present interest is in propagation velocities, the other details 
of the solution need not be considered here. It is sufficient to note that 
in a system with a positive displacement pump like the heart, probably 
the first mode is excited, and the eigenvalue (A) will be taken to be 2.40. 
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For each mode, (21) gives the functional relationships which must exist 
between the time factor (a) and the space factor (8). It is assumed that 
the disturbances are periodic in time; therefore a = iw, wherei +~W/—1 
and w is a real number called the frequency. 

The velocity of propagation (V) is obtained directly from (21); it is the 


negative reciprocal of the real part of (8/a). The result is 


V _ IM— (1—y"2)1 (22) 
Va M+T 
where 
T= [4(Vb?+ c?+56)]'4, (23) 
- 2 

b=| (1—9%s) MiG], (24) 

=|-"(1—5*2) |, (25) 
@ 

Ae 
— w - 
=T.: (27) 


The preceding six equations express the solution to the problem. 

Discussion. The first point which must be made is that, even when the 
system is completely specified with regard to the fluid, tube, and mean 
flow conditions, there is a double infinity of possible values for the veloci- 
ties of propagation of a disturbance. This is because the velocities of 
propagation depend upon the character of the disturbance itself. Each 
mode (indicated by values of the eigenvalue A) present in the disturbance 
has its own set of propagation velocities, each element of the set having 
two propagation velocities which are specified by a component frequency 
of the mode. For the remainder of this paper, it will be assumed that only 
the first mode is excited. 

Equation (22) gives the factor by which the Moens velocity (Vo) must 
be multiplied in order to produce the actual velocity of propagation (V). 
There are four parameters, in addition to the eigenvalue (A), which de- 
termine the velocity of propagation; the Moens velocity (Vo) accounts 
for the influence of the elastic properties of the tube wall, the “Mach” 
number (Mo) accounts for the mean fluid velocity, the quantity (y*s), 
which is called the “force ratio,’ brings the effects of the radial inertia 
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forces into the model. The viscosity factor (w/w) brings the effects of 
viscosity into the picture. 

The effects of the various parameters will be illustrated graphically, 
but first a few general remarks can be based on a perusal of (22), (23), 
(24), (25), (26), and (27). 

The Moens velocity does not affect the velocity of propagation in a 
simple manner; besides occurring in the left side of (22), it enters through 
the force ratio (yz). When the effects of the force ratio are discussed in 
detail, it is well to bear in mind that a change in force ratio can be wrought 
by an alteration in the Moens velocity as well as by a change in the radial 
inertia forces. 

The occurrence of the ““Mach” number (M) indicates that it is not the 
mean fluid velocity which is important, but the ratio of the fluid velocity 
to the Moens velocity. Note that when the effects of the ‘““Mach”’ number 
are discussed, the effects can be due to changes in the Moens velocity 
(elastic properties of the tube) as well in the mean fluid velocity. The in- 
fluence of the “Mach” number on the propagation velocity is quite com- 
plicated, but one of the more obvious effects is seen in a model with no 
radial inertia forces and no viscosity. For this simple model, the velocities 
of propagation are merely the superposition of the mean flow velocity and 
the Moens velocity. Equation (24) shows that M, enters into a viscosity 
term also. 

The force ratio (yz) may be considered as a ratio between the radial 
inertia forces and the elastic forces. The assumptions which were made 
in setting up the mathematical model necessitate that the force ratio be 
less than unity. Another way of approaching this limitation is to note that 
y’z can also be considered as a ratio of the radial velocities to the Moens 
velocity; it was stated previously that this ratio must be less than unity. 

The influence of viscosity on the velocity of propagation is introduced 
through the quantity T. The appearance of the ratio of viscosity to fre- 
quency in the viscosity factor is physically sound. It is obvious that the 
viscosity effects should increase as the viscosity increases. Concerning the 
frequency, with higher frequencies the boundary layer does not have as 
much time to build up as with lower frequencies, and so the viscous effects 
are less. Consideration of Figure 1 leads to the reason for the appearance 
of the eigenvalue (A) in the viscosity ratio: for higher modes, the radial 
velocity gradients, and hence the viscous stresses, are greater. The greater 
influence of viscosity in the smaller tubes, which is seen from (26), is 
predictable from knowledge of flow in rigid tubes. 

A final piece of information can be obtained from (22). Once all of the 
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parameters of the system are fixed, two velocities of propagation exist. 
From (23), (24), and (25), it can be derived that if 


M2> (1—y?z), 


then My > T, and if 
Me< tea yay; 


then M, < T. In the former case, a disturbance propagates as two waves 
traveling in the direction of the mean flow. In the latter, one wave travels 
in the positive 7 direction and one in the negative 7 direction. The positive 
sign in (22) is associated with the wave moving at the slower speed. 

In order to perceive more clearly the effects of the various parameters 
on the velocity of propagation, some computations should be made. 
Figure 2 is a plot of the velocity of propagation described by (22). In it 
the effects of all of the parameters can be seen. If actual systems can take 
on the values of the parameters shown in the figure, then the theory pre- 
dicts large viscosity and radial inertia effects. 

It is now necessary to ascertain whether or not actual systems take on 
values of the parameters which yield pronounced effects. Blood flow will 
be considered as an example. 

An average kinematic viscosity of blood may be taken to be 10 
ft.2/sec. Figure 3 gives the viscosity ratio (w/w) as a function of vessel 
diameter. The curves are plotted for two values of the mean radius ratio 
(ao)[a9 = 1 and a) = 1.414], and for two values of pulse frequency 
(w)[w = 10 rad./sec., and w = 20 rad./sec.]. For vessels under 2 m.m. 
diameter, viscosity effects are certainly important. Data given in Best and 
Taylor (1945) indicate that the fundamental frequency of the heart beat 
is less than 10 radians/second. For example, Figure 3 gives a viscosity 
ratio of 1.3 for a blood vessel of 1.25 m.m. diameter. Figure 2 shows that 
the velocity of propagation ratio is therefore altered by 12%, in this case, 
due to viscous effects. 

From the values of mean lineal velocity and propagation velocities 
given in Best and Taylor (Joc. cit.), it seems likely that the “Mach” num- 
ber (Mo) in humans is less than 0.1. Inspection of Figure 2 shows the 
effect of My on pulse propagation in this case. 

The final item to be considered is the effect of the disturbance recitehny 
on the propagation velocity. Calculation shows that the effect of radial 
inertia forces on propagation velocities in blood flow is negligible (for the 
fundamental frequency at least). It is conceivable that, if a pulse in a 
blood vessel must move a considerable mass of tissue, the effective mass 
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(p:) could be large enough to produce important radial inertia effects. 
Also, in experimental set-ups where large diameter Gooch tubing and high 
frequency pulses may be used, the effects are appreciable. A simple com- 
putation using Figure 2 and the definition of the force ratio (yz) shows 
that a tube 1 inch in diameter with a Moens velocity of 5 meters/second 
a radius ratio (a) of 1/2, a specific gravity of unity, and a thickness at 


VISCOSITY RATIO: 


=V/V, 


12 4 
FORCE RATIO=y 


VELOCITY. OF PROPAGATION RATIO 


FiGurRE 2 


1 inch hasits velocity of propagation affected by 5% if the disturbance has 


a frequency of 30 cycles/second. 

Notice, in Figure 2, that the curve for a “Mach” number of 0.2 and a 
viscosity ratio of 3 has a “dip.” This peculiar behavior is due to instabil- 
ity in the mathematical model caused by violations of some assumptions. 
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Something must be said about the difference between the results of 
this analysis and those of previous analyses. The latest work is by 
G. Karreman (1952); perhaps the errors made in setting up his model 
account for the aforementioned differences. 

The differences between the models have already been mentioned. The 
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difference in the results is the effect of viscosity. The results of previous 
analyses are doubtful because they are physically unsound: They state 
that the velocity of propagation goes to infinity as the viscosity goes to 
zero. It is obvious that as viscosity goes to zero, the velocity of propaga- 
tion should approach the results of an inviscid analysis. 
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A general theory of the drying of frozen tissue is developed and applied to the measurement 
of the drying rate of frozen guinea pig liver. It is shown that for a given temperature of the 
subliming ice crystals the minimum drying time of a piece of guinea pig liver is greater than 
the minimum sublimation time of a piece of ice of the same size and shape by a factor of the 
order of one thousand. This fact has many implications in the design of freeze-dry apparatus 
which will be developed in a following paper. 


Introduction. Rapid freezing of various biological materials and their 
subsequent dehydration at a low temperature are widely used techniques. 
For example, dried plasma is prepared in this way, as are penicillin and 
other biologicals. There is a possibility that foods might be so processed. 
The great advantage of the method is that once the material is dry it will 
keep at room temperature. The method is also widely used in cytology 
and histochemistry. If the freezing is rapid enough, chemical action is 
stopped almost immediately. After a piece of tissue is dry it can be di- 
rectly infiltrated with paraffin, sectioned with a microtome, and stained in 
various ways to localize the cellular site of different compounds (Gersh, 
1932); or the sections can be used for spectrometric studies. The method 
also shows some promise as a method of preparing cells for electron 
microscopy (Sjéstrand, 1951). 

Because of the importance of the method, there has been a great deal of 
interest in it and a good deal of literature has appeared. That work which 
has been directed toward the drying of plasma and materials with a very 
small solid component has been based on well established physical prin- 
ciples (Flosdorf, 1949; Greaves and Adair, 1939; Greaves, 1946; Kramers 
and Stemerding, 1951; Carman, 1948; Zamzov and Marshall, 1952; Ede, 
1949). However, there has been an attempt to carry this work over into 
the drying of tissue without necessary modifications. The principal omis- 
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sion in previous work on the drying of tissue (Packer and Scott, 1942; 
Scott and Hoerr, 1950; Malmstrom, 1951) has been that no adequate 
account has been given of the effect of the dry shell, which develops as the 
tissue dries. In this paper, a theoretical and experimental analysis of its 
effect is given. In a following paper, the application of this analysis to the 
design of drying apparatus will be developed. 

Theory. When tissue freezes, ice crystals form which separate from the 
organic matrix of the protoplasm. This organic matrix is compressed by 
the developing ice crystals and so forms a “‘wall’’ around each ice crystal. 
After the tissue is dried, embedded and stained, the honeycomb formed by 
these walls can be observed and is called the “‘ice-crystal artifact.’’ The 
ice crystals vary in diameter from the limits of optical resolution (0.2) 
to several micra. 

As the tissue dries it does so from the exterior inward. Each ice crystal 
sublimes from its particular space and the water molecules arising from 
its sublimation must traverse the spaces formerly occupied by exterior ice 
crystals in order to reach the exterior and be removed from the system. 
Together, these spaces and their walls form a dry shell which increases in 
thickness as drying proceeds (Fig. 1). Sublimation occurs from ice crystals 
lying at the interface between the dry shell and the frozen exterior. 

After a water molecule sublimes, it either (1) passes through the dried 
shell into the drying chamber (Path 1 in Fig. 2) and eventually reaches 
the vapor trap or vacuum pump where it is permanently removed from the 
system, or (2) it returns to the interface to be recaptured by an ice crystal 
(Path 2 in Fig. 2). 

Both paths may have variations. Initially a molecule which escapes 
from an ice crystal enters the tissue space partially occuped by that ice 
crystal (space A in Fig. 3). After making an indeterminate number of 
collisions with the walls of this space, it is recaptured by the ice crystal or 
passes through one of the holes in the walls to a neighboring space. This 
space may be another at the interface, such as D, or it may be one closer 
to the exterior surface of the tissue, such as B. Only in the latter case has 
the molecule made any progress toward escape from the system. A mole- 
cule in space B will eventually enter a neighboring space which may be 
nearer to the exterior C, nearer to the interface A, or similarly located, E. 
A molecule which returns to A may be recaptured, may return to B, or 
may enter still another space. In general, a molecule which reaches the 
surface does so only after making a very large number of collisions within 
the dry shell. However, once a molecule reaches the exterior surface, al- 
though it may re-enter the dry shell, its chances of eventually being 
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permanently removed from the system are greatly improved. As we shall 
show below, the dry shell is usually the principal obstruction. 

That the dry shell is a formidable obstruction can be demonstrated 
without resorting to experiment. We will neglect flow parallel to the inter- 
face and consider flow through a tube perpendicular to the interface, com- 
prised of ice-crystal spaces which will be assumed cubical and identical 
with end faces of area S, and holes of area s in these end faces. The space 
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Ficure 1. Diagram showing progressive increase in the thickness of the dry shell as drying 
proceeds. 


adjacent to the interface will be designated by subscript », that adjacent 
to the exterior surface by subscript 1. The mass of the molecules which 
cross unit area in direction perpendicular to that area is 


Seana 
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where P is the pressure of the vapor, T is the absolute temperature of the 


vapor, and R is the gas constant per gram. 
Thus, for the net rate of flow per unit area through the boundary 


separating the first two spaces we have 


dm _ Ss Pa—Pa (1) 
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Ficure 2. Diagram to illustrate possible paths of water molecules after sublimation from 
ice crystals at the interface between the frozen interior and the dry shell. Solid line indicates 
the physical interface; dashed, that assumed for mathematical computations. 
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where P, is the vapor pressure at the subliming interface, P, is the vapor 
pressure at the exterior surface of the dry shell, and T is the average 
temperature of the shell. 


For a shell of 1 mm. thickness, » ~ 1000; therefore 


dm 1 F; 
dt 1000 (22RT,) }/2° (6) 
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Ficure 3. Diagram of area A in Figure 3 greatly magnified to show relation of water 
vapor, ice, and ice crystal spaces. Lengths of arrows indicate direction and magnitude of 
flow of water vapor. 


Therefore the effect of a thickness of dry shell 1 mm. thick is to reduce the 
rate of drying by a factor smaller than .001 relative to the maximum rate 
of sublimation of ice at the same temperature. 

As interesting as the above argument is, it does not provide a sound 
basis for the design of drying apparatus. One must have some experi- 
mental measure of the fraction of escaping molecules. This can be ob- 
tained by considering the relation between vapor flow and heat flow. If 
we denote the escape probability by f, then the net rate of evaporation is 

UE oye VLA (7) 
dt- _JarT,) 1?" 
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where a is the coefficient of evaporation (Kennard, 1938, p. 69), T. is the 
absolute temperature of the subliming interface, P, is the vapor pressure 
of ice at T,, and R is the gas constant per gram of water vapor. 

Because of the way f appears in (7), we will call it the drying factor or 
the fractional drying rate. The reciprocal of the drying factor is the com- 
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Ficure 4. Sublimation rate of ice as a function of temperature 


bined resistance of the dry shell and the drying apparatus to the flow of 
vapor; 
als 
(2nRT,) 1/2 


is the maximum rate of sublimation of ice at temperature T,. It is plotted 
in Figure 4. 
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It is necessary to obtain the drying factor in terms of the probability 
f, that a water molecule which enters the dry shell reaches its exterior be- 
fore it is recaptured and the probability f, that a molecule which escapes 
from the dry shell is removed from the drying apparatus before it re- 
enters the dry shell. The probability that a water molecule which is 
emitted by an ice crystal at the subliming interface passes through the 
dry shell and emerges from its exterior surface into the drying chamber is 
fs. The probability that this molecule reaches the vapor trap and is re- 
moved from the drying chamber is f.; the probability that it re-enters the 
dry shell is (1 — f.). Therefore a fraction ff. of the molecules emitted by 
the ice crystal will pass through the dry shell and reach the vapor trap 
without re-entering the dry shell, a fraction (1 — f,) will be recaptured by 
the ice crystals without reaching the exterior surface of the dry shell, and 
a fraction f,(1 — f.) will reach the exterior surface, enter the drying 
apparatus, then re-enter the dry shell before they are removed from the 
system. If we assume that the probability that a water molecule entering 
one surface of the dry shell reaches the other surface without being re- 
turned to the surface where it entered is the same in both directions 
through the shell, then of the fraction f,(1 — f.) of molecules which re- 
enter the exterior surface of the dry shell, a fraction f, will be recaptured 
by ice crystals and a fraction (1 — f,.) will be reflected back into the dry- 
ing apparatus. Of the latter a fraction f, will be removed from the system 
and a fraction (1 — f.) will re-enter the dry shell. This process of molecules 
being reflected back and forth across the exterior surface of the dry shell 
can continue indefinitely, with each reflection into the drying apparatus a 
fraction f, being removed from the system. Therefore for the total fraction 
of molecules which eventually reach the vapor trap and are removed from 
the system we obtain 


f= Dssta-s9 asl". (8) 
Since this is a convergent geometric series, we obtain 
peer slew) ee 
that is, 
7 - ats mitt (10) 


This equation is discussed further in Appendix I. 
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Substituting (10) into (7) we obtain 


dm Od alos 1 


of (2s RT ee Soi 
; len as 
The maximum value of (11) is 
dm _ Oia ts 
‘dt (2xRT,) ¥?’ ser 


given by f, = 1. 

In the above discussion the drying of the tissue has been looked at pure- 
ly as a problem in vapor flow. However, heat transfer must be considered. 
In order to sublime one gram of ice, between 600 and 700 calories must be 
supplied, depending on the evaporating and condensing temperatures. 
That is, 


dQ_,dm 
7 Times (13) 


where dQ/di is the net rate at which heat is supplied to unit area of the 
interface between the dry shell and the frozen interior, and Z is the latent 
heat of sublimation per gram of ice. 

In order for heat to be supplied to the interface, its temperature must 
be somewhat below that of the surrounding environment 7... If we assume 
that the rate of heat transfer can be expanded in a power series in (J. — 
T.), we have 


a Rie Ts Se eee 


Neglecting higher order terms in (14) and combining (13), (14), and (11), 
we obtain 


dm K Giers 1 
OO seh (Tig woll amma catia tere (15) 
Pigs PREP EER f, 


where K may be regarded as a general coefficient of heat transfer. 
Combining (15) with the integrated form of the Clausius-Clapeyron 
equation, 
Peeps ai 
log B= R FMT 5 (16) 
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where P, is the saturated vapor pressure of water at temperature T,, gives 
Quel 1 


K 5 
Ee = eet a Xexp | —=Se |. (17) 
1 bigs a } c 8 


To a first approximation 
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which is obtained by substituting 


G ales 
Te T= S aRT) 


in (17). Equation (18) is not very useful for computing dm/dé, particularly 
if (T. — T.) is greater than 1° C. This is better done using (15) and vapor 
pressure tables. It does, however, illustrate the general functional rela- 
tions of the various factors which affect the rate of drying, the resistance 
of the dry shell to vapor flow, the resistance of the apparatus, the condi- 
tions of heat transfer, and the temperature of the environment (drying 
chamber). 


Since 
dx dm 


es ere ey 


where «x is the thickness of the dry shell in centimeters, and w is the ice 
content in grams per cubic centimeter of the frozen tissue, we have for 
the time needed to dry a thickness of shell «x 


t= fy Ce” (F+5 —1)dx. (20) 


Equation (20) is not very useful for computation as it stands, because both 
f. and T, will vary as the drying progresses. However, if f, is known as a 
function of the thickness of the dry shell, K/L and f, as function P of the 
apparatus, shell thickness, and environmental temperature, T, and P, 
can be computed from (15) as a function of the shell thickness and (20) 
can be integrated. Thus, a complete solution of the problem is obtained. 

We have now developed the theory needed for the experimental meas- 
urement of f,. There are several possible ways of doing this—all utilizing 
(15). The one we have used is indirect. From (15) we have 


‘ wit (2nRT,) ? 21 
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and 


17K 
=> —_— — a t: 22 
' of pe T,) dt: (22) 
if «, is thickness of the shell at the end of a drying period, 
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te 


where 0 < #’ <1. 
If we assume K/L is independent of the thickness of the dry shell 
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Thus if the temperature of the interface is measured during drying, the 
average environmental temperature is known, and the final thickness of 
the dry shell is determined, the thickness of shell for intermediate times 
can be computed from (24), the factor f from (25). Then, using (10) anda 
computed or measured value of f., f, can be determined as a function of 
the thickness of the dry shell. 

If, in addition, the weight of the tissue and, accordingly, dm/dé, is 
continuously recorded during drying, additional information can be ob- 
tained. Then 


1 r%dm 


a Bios dt=*(m, me PM (oao 
_ dm (2nRT,) 2 _ 
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Thus K/Z can be determined as a function of the thickness of the dry 
shell. It will be noted that (19) through (29) apply only to a thin exterior 
shell or to slabs of tissue. Correcting factors for other shapes will be dis- 
cussed below. 

Drying rate of guinea pig liver. We have used the above theory to deter- 
mine the drying rate of guinea pig liver. The apparatus is shown in 
Figure 5.* The copper-constantin thermocouple was embedded in ap- 


* To Macleod gauge 


2 


——> To vacuum pump 


Leads to potentiometer 


Rubber stopper 


ir Liquid nitrogen 


f¥ie—_) Thermocouple embedded 
in tissue or ice 


Ficure 5. Experimental apparatus 


proximately cubical pieces of guinea pig liver, weighing between 0.5 and 
1.0 grams. The drying chamber was evacuated to a partial pressure of air 
less than 10-4‘ mm. of Hg. The temperature measured by the thermocouple 
was recorded at short intervals. In a typical experiment, with the tissue 
five centimeters below the cold finger, the temperature of the interior 
would gradually increase from —40° C to —25° C during a three- or four- 
hour period. At the end of this time the tissue would be removed, cut into 
two with a razor, and the thickness of the clearly delineated dry shell 


* For some of the experiments a slightly modified apparatus in which the thermocouple 
wires were led in through the bottom was used. There was no significant difference in results. 
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measured with an eyepiece micrometer. The average environmental 
temperature J. was determined by measuring the equilibrium tempera- 
ture of completely dry pieces of liver. It was found to be 12° C. The re- 
sistance of the apparatus was measured, using a piece of ice instead of 
tissue. In a typical experiment with ice, it was found that during a five 
and one-half hour period a surface thickness of 0.4 cm. sublimed, the 
temperature varying from —70° C to —60° C with an average of —64C. 
This gives an average sublimation rate of .073 gms./cm.? hr. The max1- 
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Rate of drying 
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Ficure 6, Drying rate of guinea pig liver, as a function of the thickness of the dry shell. 
The ordinate is the drying factor f, and so is the fraction that a given thickness of dry shell 
reduces the rate of sublimation relative to ice at the same temperature. As is discussed in the 
text, fs is independent of temperature (see page 424), 


mum rate computed from (1) with the value of a assumed equal to one 
(Kramers and Stemerding, 1951) is 0.295 gm./cm.? hr. The ratio of these 
gives a value of f, of about .25. The value of w was taken to be 0.75 gms./ 
cm.* From these data f, can be computed from equations (10), (24), and 
(25). The composite data from fifteen experiments are shown in Figure 6. 

There are several possible sources of error in the above procedure which 
should be pointed out: 
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(1) The assumption of a constant value of K. Obviously as the shell 
increases in thickness it obstructs the transfer of heat to the subliming 
interface. Computations indicate that under our experimental conditions 
the largest part of the heat transfer is radiative (see pp. 425-26). As soon as 
a dry shell forms, it partly absorbs this radiation, and there will be con- 
ductive transfer through the shell. It will be observed that the general 
effect of assuming K constant is to minimize the variation of f, with in- 
crease in the thickness of the dry shell. The value of K/L as a function of 
the thickness of the shell could be determined to a first approximation by 
continuously recording the weight of the tissue during drying and using 
equation (29). It should be noted that this would also give by (28) a more 
direct measure of f,. The experimental details of this method would be 
more difficult than the one we have used. 

(2) The thickness of the dry shell varies from point to point, some- 
times by as much as a factor of two. This variation is partly due to differ- 
ing amounts of heat received by different parts of the tissue, partly to 
differences in the intrinsic structure of the tissue. 

(3) Error in the value of the average environmental temperature T. By 
differentiating (25) and (26) it can be shown that the values of x/ and f 
are not very sensitive to an error in T,, provided it is not a large fraction 
of T,. — T,. This is because of the way 7. appears in both numerator and 
denominator of these equations. 

(4) The assumption that the temperature of the thermocouple is that 
of the interface. There is some conduction of heat up the wires of the 
thermocouple. This was minimized by wrapping pieces of tissue around 
the lead wires. It was assumed that the temperature of these pieces would 
be approximately that of the tissue whose temperature was being meas- 
ured. The fact that the temperature of the tissue varies slowly also con- 
tributes to this error. This error is not a very large one, as can be shown by 
assuming heat constants for tissue and computing the time necessary for a 
temperature difference through the tissue to be dissipated, or by measur- 
ing the time necessary for a piece of tissue of the size used in these experi- 
ments to warm without evaporation from —100°C to —50° C. This is 
about five minutes, which is a small fraction of the total drying time. An- 
other contribution to this error is from the fact that the surrounding en- 
vironment of the tissue is not uniform, thus one surface of the tissue 
“ooks” at the cold finger and the others at the walls of the drying cham- 
ber. This gives rise to a steady temperature distribution through the block 
of tissue and so the temperature of the thermocouple will depend on pre- 
cisely where it is embedded. 


424 JOHN L. STEPHENSON 


(5) The assumption that the dry shell is completely dry. This assump- 
tion appears to be at least nearly correct. The delineation between the dry 
shell and the undried interior is very definite. As soon as the frozen in- 
terior disappears visually, the tissue as judged by weight is nearly dry 
(weight 27 to 28 per cent original). In addition if the thermocouple is em- 
bedded near the surface, its temperature rises very rapidly as soon as the 
dry shell surrounds it. 

(6) The assumption that f, is independent of T,. This will be justified 
as long as the collisions of water molecules with the dry shell are relatively 
much more numerous than collisions with each other. The vapor pressure 
of ice at —20° C is .78 mm. of Hg. This corresponds to a mean free path of 
about 50 microns, which, as is obvious from tissue microscopy, is large 
compared with distance between collisions within the dry shell. At lower 
temperatures the condition will be even better satisfied. 

(7) Biological variation in the structure of the dry shell. There is obvi- 
ously great variation from tissue to tissue and some variation for the 
same tissue in different biological states. 


Despite the difficulties, discussed above, there is no doubt that a very 
thin shell of dry tissue will reduce the maximum possible drying rate of 
tissue to a very small fraction of that of ice at the same temperature. 

The data in Figure 6, at least that between 0.1 mm. and 1 mm., can be 
fitted by a straight line. The line shown is drawn so that most of the data 
fall below it. This is to give the maximum value of f,, that is, the minimum 
value of the resistance of the dry shell, 1/f,. The line shown in Figure 6 
has the equation 

f= = ; (30) 
where d has the value 1.2 X 10-4 cm.-, x is in centimeters, and 7 equals 
.55. The name “drying coefficient” is suggested for dand “drying exponent” 
for ”. It is to be expected that other tissues will show different values of 
the drying coefficient and the drying exponent. One would expect that 
the drying exponent would have the value of one. The fact that it is less 
may partly result from the assumption of a constant value of the coeffi- 
cient of heat transfer K as discussed above. It may also partly result from 
the marked inhomogeneity of the tissue. Thus as the shell becomes thicker 
there is an increasing chance that a water molecule will find a large crack 
or blood vessel through which it will be very rapidly conducted to the 
exterior. 


For very small values of x, f., as given by (30), becomes arbitrarily 
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large, but by definition the fractional drying rate can never exceed unitv. 
A simple empirical modification of (30) which behaves properly is 


d 


a> fae 


(31) 
For very small values of x, (31) approaches unity and for others approxi- 
mates (30). (For shell thicknesses greater than 1,, this approximation is 
excellent.) From our data for guinea pig liver, we find that the fractional 
drying rate as given by (31) is 0.9 for a shell thickness of 0.2 and 0.1 for 
a shell thickness of 1u. This is reasonable in view of the known dimensions 
of the ice crystal spaces, which vary from 0.14 to 1 near the surface. 
That is, as soon as the subliming water molecules have to pass through 
two or three ice crystal spaces in order to reach the surface there is an 
appreciable reduction in the rate of drying. 

If equation (30) is substituted into (21), the drying time of a slice of 
tissue can be computed if we assume constant temperature of the interface 
and a value of unity for the coefficient of evaporation. We have 


= (27RT,) M2 w yntrl 


t 
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(32) 


Computed drying times are shown in Figure 7. For ice at temperature 7, 


— (2aRT,)” (33) 


t 
P, 


Since the value of (w/d)(x"/n + 1) is about 1000 for a shell thickness of 
1 mm., it is obvious that the tissue will require much longer to dry than 
ice at the same temperature. This does not mean that if a piece of ice and 
a piece of tissue are placed in a drying apparatus the ice will require a 
thousand times as long to dry as the tissue requires to sublime, because 
the temperature drop of the ice relative to the walls of the drying chamber 
is much greater than that of the tissue. For example, in the apparatus of 
Figure 4, we have seen that under similar conditions the temperature drop 
of ice was about 80°C, relative to the average environmental tempera- 
ture, while that of tissue was about 40° C. The drying times will be ap- 
proximately inversely proportional to the temperature drop. 

The temperature drop can be computed from (15) if we know f/-E. 
For the example of sublimation of ice given above we find 9.6 X 10-4 
gms./cm.? deg. hr. for K/L. In a typical experiment with guinea pig 
liver a value of 4.7 X 10-4 gms./cm.2 deg. hr. was obtained. In the case 
of the ice, if one assumes that the heat transfer is radiative and that the 
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walls of the chamber and the ice are black bodies, one obtains a value of 
3.6 X 10-4 gms./cm.2/deg. hr. for K/L. Conduction through the gas filling 
the chamber gives a value 1.0 X 10-4 gm./deg. cm.” hr. In this particular 
experiment #24 gauge thermocouple wire was used, which accounts for 
about 1.0 X 10-* gms./deg. hr. which is distributed over the surface of 
the interface. That these values are all of about the same order of magni- 
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Ficure 7. Drying time of slices of gui + : 
‘ guinea pig liver as a functi 
thickness. PAB on of temperature and 


tude illustrates the difficulty of a precise analysis of the heat transfer 
Nevertheless when the measured values of K/L, f., and fs are used id 
equation (15), results which are in reasonable agreement with experiment 
are obtained. Thus using the value of K/L given above for guinea pig 
liver and a drying chamber temperature of — 30° C and a value of .001 for 
f, one computes a temperature —41°C for T,. This compares with a 
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measured value of —39° C, obtained when the drying chamber was sur- 
rounded with an alcohol bath at — 30° C. As the drying chamber tempera- 
ture is lowered, the temperature difference decreases. At chamber tem- 
perature of —50° C, the predicted temperature is —53° C, the measured 
—51° C. With ice, the value of T. — T, is always much greater, because of 
the larger value of f. For example, using the values of K/L and f, given 
above, one computes that an environmental temperature of about 
60,000° C is required to raise the temperature of ice 7, to 0° C. This is of 
course incorrect, because the value of K/L does not remain even approxi- 
mately constant over such a wide range of temperature, but varies as T°. 
However, taking into account this variation, one still computes a tempera- 
ture of 1200° C. While it was impossible to verify this experimentally, we 
found that if the wall of the drying chamber was heated with a gas flame, 
the ice still remained frozen; whereas the tissue charred. An interesting 
question is to what extent it is possible to increase K by having the tissue 
in good thermal contact with the walls of the drying chamber. This is of 
some practical importance because it is desirable to minimize the tempera- 
ture drop T. — T,. The answer must await further experimental investi- 
gation—particularly of the conductivity and radiation absorption coeffi- 
cient of the dry shell. 

As noted above, the measured values of f, apply to a slab of tissue or to 
a thin exterior shell and hence so do the drying times computed from (32). 
For other shapes an approximate correcting factor can be obtained from 
analogy with the theory of the conduction of heat, if we assume the drying 
exponent 7 has the value one. For a sphere, we have from the theory of 
heat conduction through a spherical shell (Carslaw and Jaeger, 1947, p. 
199) 


dQ_T.—T: ri 
| Say Gay weal (34) 
Tosi T2 


where 7; is the inner radius of the shell, 7. is the outer radius of the shell, 
T, is the temperature at the inner surface of the spherical shell, T, is the » 
temperature at the outer surface, and & is the heat conductivity. In our 
theory d/(2rRT.)” corresponds to k when ” = 1. 
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where a is the outer radius of the dry shell, r is the inner radius, and 


if a 2 
t=w eee | (r-=)ar (36) 
8 0 
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If we define the half-thickness of a piece of tissue as the radius of the 
largest sphere contained entirely within the piece, we see that a sphere 
will dry in one-third the time of a slice of the same thickness. In a similar 
way we find a cylinder dries in one-half the time. For pieces of a given 
thickness, a sphere will dry in the minimum time, a slab in the maximum 
time. Thus the times given in Figure 7 should be reduced by a factor be- 
tween one-third and one, depending on the actual shape of the piece being 
dried; and the actual drying temperature must be computed. 


APPENDIX I 


Using (10) as a recurrence relation, we obtain for a system which is 
divided into » parts 


1 1 
Fie te oes (38) 


where f; is the probability that a water molecule which enters the 7th part 
passes into the (7 + 1)th part before it returns to the (¢ — 1)th. It will be 
noted that always f; < 1, and that when for all 7, f; = 1, f = 1. It is as- 
sumed that f; is the same for both directions of vapor transport through 
the ith part. This assumption is not valid if there are marked differences 
in cross-sectional area of the various parts. It should be noted that in the 
derivation of (10) this assumption is used only for flow through the dry 
shell, that is, for f, and not for f.. 

The term (m — 1) in (38) essentially results from the non-Maxwellian 
distribution of velocities in the water vapor because of net flow. If we as- 
sume an approximately Maxwellian velocity distribution 


f[ P;A; Pi41 Ajai |=", 


(2aRTi) 2 (2eRT ix) 2 (39) 


where P; is the pressure at the junction of the ith and (i — 1)th part, A; 
is the area of the junction of the ith and (i — 1)th part, T;is the tempera- 
ture at the junction of the ith and (¢ — 1)th part, and F is the net flow of 
water vapor, the same through all parts. 
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If we divide equations (39) by f; and sum, we obtain 


Pein a) f 
CRT)! UF sat) 


where it is assumed that P,,; is negligible and P, is replaced Dyaly 
Therefore 


Laer! 41 
vi = ( ) 


Comparing (38) and (41) we see that the assumption of a Maxwellian dis- 
tribution is justified when 1/f>> (m — 1), that is, when the drying factor 
for the system is small. 


LITERATURE 


Carman, P. C. 1948. ‘“‘Molecular Distillation and Sublimation.” Trans. Far. Soc., 44, 529-36. 

Carslaw, H. S. and J. C. Jaeger. 1947. Conduction of Heat in Solids. London: Oxford Univer- 
sity Press. 

Ede, A. J. 1949. ‘‘Physics of the Low-Temperature Vacuum Drying Process.” Jour. Soc. Chem. 
Ind., 68, 340-32. 

. 1949. ‘*The Low-Temperature Vacuum Drying Process as Applied to Green Peas.”’ 
Ibid. 68, 336-40. 

Flosdorf, E. W. 1949. Freeze-Drying. New York: Reinhold Publ. Corp. 

Gersh, I. 1932. ‘“The Altman Technique for Fixation by Drying While Freezing.” Anat. Rec., 
53, 309-37. 

Greaves, R. I. N. 1946. ‘‘The Preservation of Proteins by Drying, with Special Reference to 
the Production of Dried Human Serum and Plasma for Transfusion.” Medical Research 
Council. Special Report Series. #258. 

and M. E. Adair. 1939. ‘‘High-Vacuum Condensation Drying of Proteins from the 
Frozen State.” Jour. Hyg., 39, 413-45. 

Kennard, E. H. 1938. Kinetic Theory of Gases. New York: McGraw-Hill Book Co. 

Kramers, H. and S. Stemerding. 1951. ‘“‘The Sublimation of Ice in Vacuum.” App. Sci. Res., 
3, 73-82. 

Malmstrom, B. G. ‘‘Theoretical Considerations of the Rate of Dehydration by Histological 
Freezing-Drying.” Exp. Cell Res., 2, 688-92. 

Packer, D. M. and G. H. Scott. 1942. ‘“‘A Cryostat of New Design for Low Temperature 
Tissue Dehydration.” Jour. Tech. Methods, 22, 85-96. 

Scott, G. H. and N. L. Hoerr. 1950. ‘‘Drying: Frozen-Dehydration Method for Histologic 
Fixation.” 292-96. Medical Physics. Sec. Ed. Chicago: The Year Book Publishers. 

Sjostrand, F. 1951. ‘‘Freeze-Drying of Tissues for Cell Analysis by Light and Electron Micros- 
copy.” 177-88. Freezing and Drying. (Ed. R. J. C. Harris.) Institute of Biology: Tra- 
vistock House South, London. 

Zamzov, W. H. and W. R. Marshall, Jr. 1952. ‘Freeze Drying with Radiant Energy.” Chem. 
Eng. Progress, 48, 21-32. 


RECEIVED 2-6-53 


lee Mibew, Te (tl are 


naler, ; Kyerta 13 em 
ee 


° - 


Sin eine" a “ed wat ait 


=. - + — 
} 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 15, 1953 


A MATHEMATICAL CONTRIBUTION TO THE STUDY OF 
ORIENTATION OF ORGANISMS 


CLIFFORD S. PATLAK* 


COMMITTEE ON MATHEMATICAL BIoLoGy 
THE UNIVERSITY OF CHICAGO 


Organisms orient themselves to a stimulus by two general methods. One method is by di- 
rected orientation (taxis); the other is by undirected locomotory reaction (kinesis). An equa- 
tion, and the methods for finding the necessary parameters of this equation, is derived for the 
distribution of organisms within a container, with the following limitations: (1) the organisms 
have no accommodation, (2) they are always active, and (3) the stimulus changes slowly with 
position. Necessary modifications of the equation are then derived, so that the last two limita- 
tions may be eliminated. The equation cannot be solved exactly because of its complexity; 
hence an approximation method must be used. This method is discussed, an approximate 
solution is found, and a time constant for equilibrium to be established is derived. Applications 
to various experiments in the literature are then made with fairly satisfactory results. A new 
interpretation of the theory of klino-kinesis with accommodation is found upon application 
of the equations developed to experimental work. Further limitations and uses of these equa- 
tions are then discussed. 


Introduction. Organisms orient themselves to the effect of stimuli (such 
as heat, light, humidity, etc.) in two ways. One is by a directed orientation 
reaction (taxis), in which the direction of motion of the organism is 
influenced by the stimulus. This is a more general definition than that 
commonly used (i.e., that the organism’s motion was directed toward or 
away from the source of the stimulus) in that it also includes the term 
“transverse reaction,’ which means the organism moves at an angle to 
the line joining it to the source of the stimulus. Another method of orien- 
tation is undirected locomotory reaction (kinesis) in which the average 
speed, the average rate of turning, etc. of the organism, but not the direc- 
tion in which it moves, are dependent upon the stimulus. These methods 
of orientation have been thoroughly reviewed and discussed by G. Fraenkel 
and D. L. Gunn (1940). 

For a review of previous mathematical attempts to treat the problem 
of orientation of organisms, the reader is referred to a paper by the author 
(1953) which will hereafter be referred to as I. In that paper equations 
were derived characterizing the movement of particles under the influ- 

* This work was done while the author was Public Health Service Research Fellow of The 
National Institute of Mental Health, Federal Security Agency. 
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ence of an external force and which showed persistence of direction. These 
equations may be applied to the problem of orientation of organisms by 
first showing how the terms of these equations are related to the organ- 
ism’s movement and how their various parameters may be evaluated for 
any specific experimental situation. Certain limitations of these equations’ 
applicability will be pointed out and methods will be developed to remove 
some of these restrictions. Due to the complexity of the equations, only an 
approximate solution may be found. The equations derived will then be 
applied to experimental data found in the literature in order to illustrate 
the methods outlined in this paper, test the validity of the equations, and 
also show how further experiments should be carried out so that the 
equations would be applicable. 

Some modifications of the definitions given in I will be necessary to fit 
the present situation. In this paper the term “persistence of direction” 
implies that when an organism turns it will not turn with equal probability 
in all possible directions, even when there is no taxic component to its 
motion. For example, the organism after turning may travel more often in 
the same direction in which it was traveling before turning than in the 
opposite direction. By external force we mean the “force” of the stimulus 
which causes a taxic component or an avoidance reaction in the organism’s 
motion. The avoidance reaction will be discussed later in this paper. 

In I, a particle (which we consider as representing an organism in this 
paper) was considered to have a net average motion in a given direction 
for two reasons: one was that the organism turned more often in a given 
direction due both to persistence and the external force; the other was that 
the organism traveled further in a given direction than in the opposite 
direction due to an external force. 

The following notations were introduced in I. The parameter e,,(X, #) 
represents the coefficient of external bias in the «,-direction if the organ- 
ism turns at the point X at time ¢, where X denotes the point (x, y, z) and 
%1, %2, “3 represents the coordinates «, y, z respectively. That is, it is a 
measure of the probability that an organism turns in the x,-direction over 
and above the average probability that it turns in any direction. Also, 
¢o(X, #) represents the average speed of the organism; 79(X, #) represents 
the average time between turns of the organism; (ér))(X, #) represents 
the average distance the organism travels between turns, which may also 
be denoted as the mean free path; (@72)o(X, 4) represents the average of 
the squared distance which the organism travels between turns; ¥(X, #) 
represents the coefficient of persistence of direction, which is a measure 
of the probability that the organism will tend to persist in traveling in the 
same direction after turning as it was before turning; and (erj.[X, 8 
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represents the average distance between turns the organism travels in the 
*;-direction over and above the average distance the organism travels in 
any direction—all parameters are for organisms which turn at point X 
at time ¢. 

Thus we see from the discussion above that both e,,(X, #) and (C7) a,(X, 8) 
are measures of the taxic component of the organism’s motion. That is, 
if there is a taxic component to the motion, then the six parameters 
(¢ = 1, 2, 3) above will not all be zero, and conversely. In the derivation 
of the equations in I, the assumptions that (¢r).,(X, !) « (ér)o(X, #) and 
that «.,(X, #) < 1 are made. This implies that the taxic component of 
the motion is small compared to the normal random motion of the organ- 
ism. This is precisely what we wish it to be for this paper as we may 
easily and quickly note a large taxic component in the motion of the organ- 
ism and thus we are only interested in ascertaining the small components 
which escape casual detection. We shall discuss in the next section how 
the values of the various parameters in the equations which we shall use 
may be determined from experimental data. 

The problem of accommodation will not be treated in connection with 
the various types of orientation with which we shall attempt to deal. By 
accommodation (or adaptation) we mean that the organism’s behavior at 
a given point and time not only depends upon the particular point and 
time, but also upon the organism’s history. The equations developed in I, 
which will be used in this paper, were derived under the implicit assump- 
tion of non-accommodation. However, accommodation has been treated 
in I for the case of persistence of direction which depended only upon the 
organism’s direction of travel immediately prior to its turning at the 
given point. The method which was used in the derivation of the co- 
efficient of persistence of direction, ¥(X, ¢), could be used to treat other 
cases of accommodation which involve only the previous step. However, 
such a treatment will not be discussed in this paper, primarily because 
almost all of the cases of accommodation which arise experimentally con- 
cern the organism’s history for more than just one previous step. 

In the derivation of the equations in I, infinite regions were postulated. 
In any actual experimental situation, however, the organisms will be en- 
closed in a finite container of length L. However, we may apply these 
equations, with negligible error, to finite regions if we assume that the 
average distance that the organism travels is much less than the length 
of the container, that is, 

(cr)o(X, 4) KL, (1) 


and that the probability of an organism traveling a distance greater 
than, say, ten times a mean free path is negligible. Then except for 
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those organisms which are very close to the edge, we may ignore the 
effect due to the finite size of the container. Due to (1), the small size of 
the taxic component of the organism’s motions, and the fact that, 
physically, we can expect the number of organisms to vary continuously, 
the number of organisms close to the edge will represent a negligible frac- 
tion of the total number of organisms. We may therefore conceptually 
replace the actual organisms near the edges by hypothetical ones that 
obey the equations and, due to their negligible number, should not ex- 
perience any difficulty. However, it is found experimentally that the 
organisms in certain experiments may tend to travel along the edges to a 
much greater degree than would be expected from the theory. We shall 
not try to treat this effect but shall merely ignore all such organisms in our 
counting of the number of organisms in each portion of the container. This 
point will be treated in more detail in the section on the discussion of 
experiments. 

Although the equations developed in I will apply to any experimental 
situation which obeys the.assumptions listed in the derivation of those 
equations, we shall limit our considerations to the type of experiments 
that have actually been performed. We shall consider organisms that 
either move about in space (three-dimensional case), are confined to 
movement about on a plane (two-dimensional case), or are confined to 
movement along a narrow path (one-dimensional case). We shall further 
limit ourselves to the situation where the stimulus changes along only one 
dimension, which we will call the «-axis. This limitation is made solely to 
simplify mathematics and does not restrict the generality of the method. 
The point X may thus be effectively given by its x-component, and we 
note that 

(er) y(, i) = (er) s(x, t) = g (x, t) =e,(x,t) =0. 
Hence if » denotes the number of dimensions the organism is free to move 
about in and if N(«, #) denotes the number of organisms in a unit region 
about the point x at time ¢, (36) of I may be written as follows 


CNet lle Gat 1+¥(2 1+v(2 22-1) Dae N (cr? ae 
art ot el to slo | 


TO 


~ [24 (u4ve (2) 


¥_(¢r)o 8 [(e*r)0]\ (Er) 0 
+Toppaehe ae =: TO |y : 


ORIENTATION OF ORGANISMS 435 


where all of the terms are functions of position and time. This is the 
principal differential equation to be used in this paper. Due to the com- 
plexity of the above equation, only the steady-state solution can be found, 
i.e., the solution for the case in which the number of organisms in each 
region of the container does not change with time. In order for this situa- 
tion to occur, we assume that the parameters will approach a steady-state 
value. Let one end of the container be taken as the origin, and let 


2{ (er)e(x) + [1+y¥ (*)]e(«) (er) o(«) } [1—W (x) ] 
T 2 a 


(e272) o(x) 


1 


q(x) = 


(3) 


(cr)a(x) d (c? 7?) 9 (x) 
jity (x) [2 Ate? _ 


(77?) o (x) 


r(x)= 


If K is a constant which is found from the boundary condition, then the 
total number of organisms, V7, in the container is given by 


Ne= f N(x) dx, (4) 


and, using the fact that in the steady state the net number of organisms 
crossing a line through any point x is zero, we then have 


ean y cay ew ff ale) +r (alae (5) 


To (x) 


In the alternative chamber type of experiment (Gunn and Kennedy, 
1936) the stimulus is kept constant within each of the two regions of the 
container but differs between the two regions. We will denote the two 
regions by ¢ and y. For this case, the following points should be noted. 
First, since the aspect of the stimulus which affects the organism may not 
be the intensity of the stimulus but its gradient, etc., what we mean by 
a constant stimulus is the constancy of this aspect. In all of the experi- 
ments which have been reported in the literature, this aspect has been 
taken to be the intensity alone. Second, there is usually a sudden change 
of the stimulus in the center of the container, which occurs within a very 
short length. This situation contradicts assumption 5 of I concerning the 
slow change of the functions. However, (2) will be valid for the regions 
inside of the container which are some distance from the center (i.e., at 
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least a few mean free paths away) and hence (5) will hold for these regions. 
For this case, the solution to (2) may be written as 


: : is 
Ne etn eRe: (g(a) +r(x)]del, OS"<5 
To («) . (6) 
Seo tft (g(a) +r(a)ldxt, F<eSL 
AGES AIS L/2 


To (x) 


where equation (4) is one of the boundary conditions necessary to find the 
values of K, and K,. The other necessary boundary condition will be dis- 
cussed in the section on the avoidance reaction. 

Pure ortho-kinesis as defined by Gunn ef al. (1937) solely involves a 
change of the average speed with change in stimulus intensity. However, 
we will omit from their definition any factor involving the inactivity of 
the organisms (cf. Ewer and Bursell, 1950). Thus ortho-kinesis would be 
measured by é)(x, ¢). We shall also adopt Ewer and Bursell’s concept of 
klino-kinesis without accommodation. According to this definition, pure 
klino-kinesis solely involves a change of the average rate of turning under 
different stimuli. Hence klino-kinesis would be measured by 7o(x, #). But 
we see that equation (2) involves such quantities as the mean free path 
(cr) o(m, 4), etc. Thus what is important are the values of 70(x, #), (cr) o(x, 2) 
and (¢2r2)o(, 4), and not the value of éo(x, f) or 7o(x, #) alone. It is not 
that the classification of the organism’s movement into ortho- or klino- . 
kinesis is unimportant, but as far as the equations are concerned this 
classification is not relevant. In this connection it may be pointed out that 
pure ortho-kinesis as defined above implies 


(cr)o(x, t) =700(x, t), 
(c%72) | (x, 4) = r8c2(x, 2), (7) 
to(a, t) =7, 
while pure klino-kinesis implies 
(cr)o(x, t) = cro(x, 2), 
(c8r?) o (a, t) = c72 (x,t), (8) 
to (x, t) = 79(x, t). 


If we inserted either of the above into (5) or (6), we would note that there 
is very little simplification introduced and hence we see that little is 
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gained mathematically by subdividing the motion into pure ortho- or 
klino-kinesis. 

Evaluation of parameters. The equations which have been derived in I 
contain various parameters, such as (¢7)o(X, 2), (c7).,(X, 2), €z,(X, #), etc., 
which must be known to the experimenter so that he may solve the equa- 
tions for the number of organisms in each region. In order to obtain data 
from which the above parameters may be determined, observations may 
be made by either of two general methods. The first of these is the “cross- 
section method” (Bentley e¢ al., 1941) in which the observer makes his 
observations at certain prescribed time intervals, each observation being 
instantaneous. Average values of the various parameters may be found by 
this method, as desired. Its major advantage is that many experiments 
may be carried out simultaneously, while some of its disadvantages are 
that the experiments must last for a longer period of time than those done 
under the “continuous observation method” described below, parameters 
such as 70(X, #), W(X, #), etc. cannot be found by this method, and the 
precise behavior of the individual organism cannot be found. This last 
disadvantage is important only if the experimenter is interested in seeing 
if the organism obeys the assumptions leading to the derivation of the 
equations. If this point has once been settled, then this disadvantage is 
obviously of no consequence. If the organism changes its characteristics 
very much (i.e., adapts) during the course of the experiment, this method 
is not applicable. The length of time between observations should, for 
maximum efficiency, be large enough so that there is little or no correlation 
between any two successive observations. The total length of time 
needed for the experiment will be determined both by the characteristics 
of the organism and the length of time between observations. 

The second method, which may be called the ‘continuous observation 
method,” involves continuous observation of the organism. In this meth- 
od, only one experiment may be performed at a time. In order to find how 
often the organisms turn, by how much they do turn, and their precise 
behavior, this method must be used. If the only information required is 
such facts as the average speed, the amount of inactivity, etc., then either 
of the two methods described above may be used. Also, the results of the 
two methods may be compared with each other in order to check the in- 
ternal consistency of the theory. 

In the derivation of (2), the assumption was made that there is no 
interaction between the organisms. This may be easily achieved by allow- 
ing only one organism in the container in which the experiment takes 
place during each experiment. The data may then be examined either for 
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each separate organism or lumped together for many experiments, each 
with a different organism. The results of both methods may be treated by 
equation (2). The continuous observation method may perhaps be more 
useful since results are found for the average organism of the species and 
there is more data, permitting results that are more valid statistically. 

The organism may adjust to the total environment of the experiment in 
such a way that the various parameters at each point may change with 
time until the organism reaches a steady state. This is not the same as ac- 
commodation, discussed previously, since the organism adapts to the 
total environment and not to any particular portion of it. If we are inter- 
ested in the value of the various parameters at different times, then we 
have to find this average value by considering the average over the differ- 
ent organisms in the various experiments for the same time. This con- 
sideration need not apply in the steady state, because we may then find 
the average value of the parameter fora single organism, since this aver- 
age value is the same for all times. Therefore, for the steady state, we may 
average over the values of the parameter found whenever the organism 
appears at the desired position. 

Certain of these parameters are extremely simple to ascertain. For ex- 
ample, to find (ér)o(X, #) the experimenter need only observe the average 
length of the organism’s path, regardless of direction, if the organism last 
turned at the point X at time ¢. The value of 7)(X, #) may be found by 
observing the average length of time a particle travels between turns, while 
(C2r2)9(X, t) is found by averaging the square of the lengths of the paths, 
regardless of direction, if the organism last turned at the point X at time ?. 
For (¢r)z,(X, ¢) the experimenter need only note one-half of the net average 
extra distance which the particle travels in the positive «,-direction over 
that of the negative x,-direction. 

The method of determining e,,(X, #) and W(X, #) is not as obvious as 
the methods listed above. In order to find these parameters, the following 
approach, which is illustrated for the two-dimensional case, may be used. 
In the notation of I, if an organism enters and turns in a region about the 
point Y at time ¢ after traveling in the direction given by the angle 4, 
the probability that it will then leave the region by traveling in the direc- 
tion given by the angle @ is 

E(62| VY, t) I(62| 6:, Y, t) d. 


JE Co Y, t) I(03| 0, Y, t)d 0, (9) 


Since we are considering the situation wherein ¢,(Y, #) is 0, J (02|61, Y, 2) 
depends only upon the angle a between the incoming and outgoing paths 
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(a = 0; = 62), I'(a| ie t) = I'(-a| ¥e is and E(6» | 4 t) > 47[1 
+ €.(Y, ¢) cos 6:]; therefore (9) becomes 


[1+e.(Y, t) cos 6] I’(a| Y, #) d 6, 
Fe eh a bee (ls t) cos 6; ; 


Since we are only interested in the coefficient of persistence of direction, 
¥(Y, ¢), and in e,(Y, é), we may consider certain averages. Let a wavy line 
across the top of a function denote the average value of the function for a 
given value of 6;. For an organism which enters the region about the point 
Y at time # from the direction given by the angle 6,, the average value of the 
cosine of the angle which the organism turns would be 


bie V4) cos 05} (a 7,2) 
Ley, the. (Y, t) cos 6, 


2x 
CS a(Y,) = f COS a 
0 


Qr (10) 
CYS?) +eCri ap ws af cos? al’ (al Y, t) da 
— 0 . 


1 -WieY a he, CY t) cos 6, 


If we similarly find sin a(Y, #), cos 62(Y, #), and sin @,(Y, #), then solve 
each of these equations for 


2a 
fi cos? al’(a| VY, t)da 
0 


and equate the solutions, we find 


[1+we, cos 6;] cosa —y 
cos 6, 
_e,sin 6:— [1+ We, cos 6;] sin a 
- sin A; 


(11) 


[1+ We, cos 0:] cos 02—y cos 0;—e, sin” 41 
cos 24; 


aw , sin 20 
[1+ We, cos 6;] sin 0.—y sin Pap oror aa 


sin 2 0, 


where all of the functions, except those involving 4, are functions of posi- 
tion and time. From the derivation of the above equations, it can be seen 
that they are valid only if e,(Y, ¢) or none of the denominators of the above 
fractions is zero. If e,(Y, ¢) is zero, then in place of (11), all the numerators 
of the four fractions should be set equal to zero and the resultant equations 
may then be easily solved [see eq. (13)]. If any of the denominators 1s zero, 
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due to a particular choice of 6, then the numerators of these fractions 
should be set equal to zero and the resultant equations solved. 

In expression (11) there are three independent equations for two un- 
knowns, and hence the method can be checked. Since we have assumed 
in I that (Y, A)e,(Y, t) «0, the relations in (11) simplify to 


cosa(Y,t)—w(Y,t) _e(Y, t)sin 0:—sin a (Y, 2) 
cos 6; sin 6, 


_ cos 62(Y, t) —W(Y, t) cos 6,—e,(Y, t) sin? 0; 
=: cos 2 6, 


(12) 


Ss in 26 
sin 6.(Y, t) me ed t) sin Ae te = 


J 


sin 26; 


which may be handled much more easily than (11). The statement follow- 
ing (11) also applies to the above expression. 

If, by the nature of the experiment, we can physically be assured that 
¥(Y, ¢) and e,(Y, #) are constant within a certain sector of the container, 
we may then simplify our task by considering all of the organisms that 
enter parts of this sector from the direction given by the angle 0, and then 
turn in these parts. 

If «,(Y, ¢) is 0, then, as stated before, if the numerators of the fractions 
of (12) are set equal to zero and the resultant equations solved, we have 


cos 62.(Y,?#) _ sin 42(Y, #) 


Y,i)= es = a = 
¥ POO Ek a cos 4; sin 6; 


(13) 
0=sina(Y,?). 


If the organism does not possess any persistence of direction, then J’(a| Y, 2) 
= 1/27 and y(Y, #) = 0. Hence we may evaluate 


Qr 
J ‘cost al’ (a|Y,t)da 
0 


and find the values of e.(Y, ¢) directly from (10) and equations analogous 
to it which are not shown in this paper. We then find 


ay('¥, ¥) aero SOR og Sieh. 


cos. 63° «= Bah Oa a 5 2 cos #2(Y, #) , 


(14) 
0 =sin 0.(Y, ar 


It may be too difficult experimentally to find a large number of organ- 
isms entering the region from a direction given by the angle 6;, and there- 
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fore we may consider the average of the various functions for all 6,’s. 
We shall illustrate the general method for the special case of eV i=, 
Since cos a(Y, #) and sin a(Y, é) are independent of 6,, and if we average 
these over all ;, letting the double wavy line indicate the average value 
of a function for all 6;’s, we have from (13) 


¥V(Y,t) =cosa(Y,?#), 


(ee 
0=sina(Y, 2). 


For the two remaining relations of (13) we may average over the 6;’s by 
means of the “least squares” method. Thus we have 


> [WC Y, t) cos 6;—cos 62] 2= minimum , 

Z[¥ (CY, ¢)sin 6;—sin 6.] 2= minimum , 
where the summation extends over all of the experimental data. If we 
then differentiate the left side of the above equation with respect to 


V(Y, t) and set these values equal to zero (since we are looking for a 
minimum), we have, after simple rearrangements, 


cos 6; cos 4 
D cos? 6; 


ren eee 
(16) 


eae sin 6, sin A 


2 sin? 6, 

From our knowledge of the physical basis of the experiment, from data 
on the relations between the averages, and from our measurements on 
(cr) a(X, #) we may easily ascertain whether there is a taxic component 
in the organism’s movement. Another check of this is to perform the 
experiment in a container long enough so that there will be a crowding up 
of the organisms at one end [see eq.(5)]. In the next section another meth- 
od will be outlined for determining if the organism has a taxic component 
to its motion. 

Avoidance reaction. In many experimental situations it has been found 
that when an organism reaches the dividing line between two different 
regions it may suddenly react to the new region as if it was trying to avoid 
it. For example, if Paramecia are swimming in a container divided into 
a cool and a hot region, whenever one encounters the hot region it will 
respond by suddenly backing up and turning. It will continue doing this 
until it eventually turns into a direction along which it may travel with- 
out encountering the hot region (Jennings, 1904). Excellent reviews of this 
type of behavior are to be found in S. O. Mast (1911), G. Fraenkel and 
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D. L. Gunn (1940), and D. W. Ewer and E. Bursell (1950). In the present 
paper the term “avoidance reaction” (Jennings, 1906) will be used to 
describe this type of response. P. Ullyott (1936) studied the reaction of 
the organism Planaria to various light stimuli and came to the conclusion 
that the avoidance reaction could be explained as simple klino-kinesis 
with accommodation, that is, the organism upon encountering the region 
of rapid change increases its rate of turning tremendously, so that its 
T)(X, #) for this narrow range is very small. However, as Ewer and Bursell 
(1950) point out, this explanation will not suffice for many of the situa- 
tions which have been studied experimentally in which the organism stops 
at the boundary and then, after waving its antenna of forelegs about, 
turns around and walks away. Another experiment in which it was found 
that this explanation is not valid is that of V. B. Wigglesworth (1941), 
about which he says: 

... when the insect approaches the boundary obliquely. It is true that it may 
occasionally turn into the adverse zone—as would be expected if the direction of 
turning were entirely random as in “‘klinokinesis’—but in the vast majority of 


instances the louse inclines away from the boundary and back into the favour- 
able side. 


Thus, in this paper, we shall consider the avoidance reaction as distinct 
from klino-kinesis and attempt to derive certain relations for it for the 
steady state. 

Mathematical considerations would also lead to this conclusion. Avoid- 
ance reaction, by its very nature, is usually concerned with changes of 
stimuli which occur within a space of less than one mean free path of the 
organism. Thus assumption 5 of I, concerning the slow change of the 
functions, is completely violated and the mathematics which have been 
derived are not applicable to this case. That is, with the formulae which 
have been derived we cannot tell if the organism “avoids” the region by 
means of klino-kinetic behavior or by some other method. To utilize the 
mathematics available, the concept of avoidance reaction will be used for 
which further formulae may be developed; we will not assume that these 
reactions may be explained by klino-kinetic behavior. Of course, these 
reactions may actually be due to a klino-kinetic behavior, but this fact may 
be determined by other methods, such as actual observation of the organ- 
ism’s responses at the boundary. 

We do not have to assume complete avoidance for the purposes of this 
paper. The organism may tend to shy away from the region which it is 
avoiding, but it may still enter it in a certain fraction of its encounters 
with the boundary. This point is most clearly brought out by Wiggles- 
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worth (Joc. cit.) in the statement quoted above. However it should be 
noted that the term avoidance reaction, as it will be used in this article, 
implies an asymmetry. That is, if the organism in going from one region 
into the other avoids the new region by the same amount regardless of 
whether it is going from the first into the second region, or vice versa, 
then the steady-state distribution of organisms will be the same as if there 
was no avoidance reaction at all, providing there is not complete avoid- 
ance. An analogy with a semi-permeable membrane in a solution may 
provide an example of the above. If the membrane is equally imperme- 
able in both directions to the solute, then the equilibrium, i.e., steady 
state, distribution of the solute particles will be the same as that which 
would have existed if the membrane were completely permeable to the 
solute particles. 

Several types of taxic responses can be classified under the term avoid- 
ance reaction (as used here). Klino-taxis [analogous to the older term of 
Phobotaxis (Fraenkel, 1931)], which is the comparison of the stimuli 
which are successive in time, and Tropotaxis with anteroposteriorly 
placed receptors (Ewer and Bursell, 1950), which is the comparison of the 
stimuli which are successive in space, will be treated by the methods of 
the avoidance reaction. This will be done not only because we wish to treat 
all avoidance reactions by one method, but also because in this type of 
behavior the organism’s movement at the boundary would depend on its 
previous position; hence the organism would be exhibiting accommoda- 
tion. This might be treated by use of the persistence of direction method. 
However a simpler and more general method is described below. 

In order to derive the relations which are desired, it is first necessary to 
define certain functions and point out how they are related to the problem 
at hand. Consider a two-dimensional container divided into two regions, 
¢ and y (generalization to the one- and three-dimensional cases, as can be 
seen, is obvious). As before, we shall assume that the functions vary only 
in the x-direction, and hence all points may be effectively given by their 
x-coordinates. Let 7(&, &) represent the probability that an organism 
which crosses the line through the point &; (all lines are taken to be per- 
pendicular to the x-axis) traveling in the direction toward the line through 
the point £ will cross the line through £ before it again recrosses the line 
through §&. Similarly for 7(&, &) where the points & and & are inter- 
changed in the definition above. Let v(&, &), the avoidance function, be 
defined as 


_ 19 (E1, &2) 17) 
v (£1, &2) Bey Gees en Geucen , ( 
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Let us choose the point & to be in region ¢ and &, to be in region y, both 
points being many mean free paths distance from the dividing line between 
the two regions. If there is a positive avoidance reaction, i.e., the organism 
tends to avoid traveling from region ¢ into region y, then it is physically 
obvious that 7(&, £,) will be less than the value which it would have had 
if there were no avoidance reaction, while n(é,, &) would be greater. 
Therefore, »(&, £,) will be less than its value for no avoidance reaction. 
Similarly, if there is a negative avoidance reaction, then »(£;, £,) will be 
larger than its value for no avoidance reaction. Hence if some method is 
found to calculate »(&, &,) for the case of no avoidance reaction, then by 
simply observing 7(&, &,) and n(é,, &) and thus finding »(&,, £,) experi- 
mentally, it is possible to tell if there is an avoidance reaction, and, if so, 
how much. With this goal in mind, some further definitions and deriva- 
tions will now be made in order to obtain certain general results and then 
derive a formula for a specific example to illustrate the methods involved. 

Let (&, &, +) be the number of organisms that cross the line through 
per unit time traveling in the direction toward the line through & and 
then, before crossing the line through &, again in the opposite direction, 
cross the line through &. In an analogous manner, let (&, &, —) be the 
number of organisms that cross the line through & per unit time traveling 
in the direction toward the line through &, and then cross the line through 
£; again in the opposite direction before crossing the line through 
£. The functions m(&, &, +) and m(&, &, —) are defined in a manner 
similar to the functions above, with & and & interchanged. Let the num- 
ber of organisms that cross the line through & per unit time while travel- 
ing in the direction of increasing x be N.(&) and in the direction of de- 
creasing « by V_(&,). If we choose & > &, then an organism which crosses 
the line through & while traveling in the direction of increasing x must 
either cross the line through & before it again crosses the line through &, 
or it must cross the line through & before it crosses the line through £. 
An organism which crosses the line through &, traveling in the direction 
of decreasing «, must either have crossed the line through & while travel- 
ing in the direction of increasing « and then crossed this line again before 
crossing the line through &, or it must have crossed the line through & 
while traveling in the direction of decreasing x and then before crossing 


this line again must have crossed the line through &. Thus in the steady 
state 


N+ (&1) =n(&, f+) +n(&, Exim, 
\ (18) 


N- (1) =n(&, £o,— ) +n (és, Ear a 
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Since, in the steady state, 


N+ (§1) = N_ (£1), (19) 
therefore, substituting the above equation into (18), we find 
(£1, £2, +) =n (ko, &1,+). (20) 
From the definitions it is obvious that 
ra _ (hy, te) 
n (£1, &2) ous ae 8%) tek (21) 
and 
bey Sein) 
n (2, &1 aes te (22) 
Then from (17), (20), (21), and (22) we have 
TWEE) 
v (£1, &2) NG)" (23) 


To find V,(£), a method completely analogous to that used in I may be 
used. Thus we have 


N+ (£1) Mpftity ais nany Woaeem ees Seite 


eee Sete 


(24) 


Using the same expansions and approximations used in I, equation (1), 
and the fact that since we are dealing with a steady state, the net number 
of organisms that cross a line through the point x per unit time is zero, we 
find that (24) yields 


Neacesiae > N (é,) obey (25) 
To (&1) 


Inserting (19) and (25) into (23), since & is arbitrary, we find 


N (&2) to( £1) (er) 0 (és) (26) 
N (£1) ro(&2) (er) 0(&1) ; 


If the two regions ¢ and y are identical, then from the above discussion 
we see that v(&, £,) should be equal to one. However, if the two regions 
are different, then even though the organism does not actually exhibit 
an avoidance reaction at the boundary between the regions, it may still 
possess an avoidance function that is different from one. In order to de- 
termine whether the organism does or does not have an avoidance reac- 
tion, the theoretical value of »(&, £,) which the organism should possess 
under the experimental situation, providing it does not have an avoid- 
ance reaction, must be calculated. We shall let the subscript , represent 


v (£1, E>) — 
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the theoretical value computed under the assumption that there is no 
avoidance reaction. 

A specific example will be considered in order to illustrate the general 
methods which may be used. As stated previously, one of the reasons for 
the introduction of the avoidance reaction concept was the difficulty of 
dealing with sudden changes in the stimulus, i.e., those taking place 
within a region of less than a few mean free paths. Therefore we shall con- 
sider the stimulus as changing gradually over many mean free paths. That 
is, even though the stimulus may not change gradually, we shall consider 
that its effect will be the same as if it did and any experimental devia- 
tions from the theoretical results may be attributed to an avoidance re- 
action. Then, since it has been assumed that the container is much greater 
that a mean free path [eq. (1)] we may solve the equations by ignoring any 
of the gradual stimulus changes and treating it as if there was a sudden 
jump in the stimulus. The reasoning for the above is completely analogous 
to that used in the discussion on ignoring the number of organisms near 
the edge of the container. 

We will assume that the taxic component of the motion, the coefficient 
of persistence of direction, and the mean free paths are constant within 
each region (although they may be different in the two regions) and that 
the regions are equal in size. From the assumption of the constancies of 
the functions, it can be seen that (3) reduces to 


g(x) = qe, 0s0<5, 
= gy, F<2SL, 
for(e)dz=0, 052<2, 
0 2 
Ve 
= am < 
Ks 7 SSL, 


where q;, gy, and s are obviously constants. The values of g and g, may 

be easily found, while finding the value of s may be more difficult. This 

will be discussed later in this section. Thus if & is taken to be in region 

§(0 < x < L/2) and , is chosen to lie in region y(L/2 < x < L), then 

ie inserting the values of (27) into (5) and these results into (26), we 
ave 


Ven (Es, fy) 


_ (ér)0(&y) (¢?r?) 0 (Er) 


— ee eE) CN ey OP LoL] + ofS ee] est. 
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If we now choose two points fi, and £, so that they both lie in region ¢ 


and (f — fr) = (Z/2 — &), then by the same argument as above we 
would have 


16, 
rn (Erts E22) =exp | ax[5— Er] t. (29) 
Analogously, if any two other points &, and &, are chosen so that they 
both lie within region y and (&, — &,) = (é, — L/2), then we would have 
rn (Em Em) exp fa —Z] (30) 
Hence, using the notation 


Vy, (Er) =v, (E1r, for), 
(31) 


¥en, (Ey) Sy, (Er Ear), 


where £1;, £, &,, and £, are as defined above, we have, after substituting 
(29), (30), and (31) into (28), 


ee on (edo (E,) et. (32) 


(cr) 0 (&y) 


Vin (Es = ae 


If the value of s can be found, then only the values of the mean free paths 
and of vi_(&-) and »,(é,) are needed to determine the value of vi, (&, &). 
If the assumptions are correct, the functions v,(é) and »,(é,) should 
equal their values found experimentally, as there should be no avoidance 
reaction within each region. If it is found that the experimental value of 
v(&, &,) does not agree with the theoretical value, then an avoidance re- 
action may be assumed to be present. We shall let the subscript ex denote 
the value found experimentally. In order to obtain a quantitative measure 
of what may be called the true experimental avoidance reaction we will 
use the function 2,(&, &), which is defined as 


* _Yex (£1; E>) 33 
Pox ( £1 2) ~ Ven (E1, Ee) * ah 


In some experiments, the avoidance function cannot be determined 
directly, i.e., by the use of (17), due to a lack of proper data, but must be 
found by the use of (26) where the values of the parameters on the right- 
hand side of the equation are determined experimentally. For this case, 
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if we let the subscript a) denote the calculated value, we may then use 
the function v%,)(&, &), which is defined as 


_ Yea (E1, £2) 34 
Vent (£1, £9) ack bis Es\-4 ( ) 
Inserting (26) into the above equation, we find that 
& (£2) | 
* ee N (£1) ex 35 
Veal (£1; E>) a pea ] : ( ) 
N (£1) Jin 


As is obvious, if (26) is valid, then 
Veot (E1, £2) =Vex (£1, £2), 


an expression which may also be used to test the validity of equation (26) 
and hence of the method itself. 

It should be emphasized again that (32) holds only under the assump- 
tions mentioned. These assumptions should not be construed to be gen- 
erally valid, since they were chosen merely for their mathematical sim- 
plicity and illustrative uses. In fact, they may only be true for very special 
cases, such as for no taxic component to the movement. 

Another example will be given to illustrate how the value of s may 
be found. For this situation, we shall assume that the coefficient of per- 
sistence of direction has the same value in both regions and that the prob- 
ability distributions of the c’s and r’s are narrow so that (c2r2)o(x) 
= (¢7)?(x). For this case (3) becomes 


ay s53 a 
(cr)o(x 
Pea a (er) 0 (x) Taek (6 “| opm (Tee Dil. 


Therefore, inserting the above into (27), we find 


doth vilken atte) 
va tte os lal tela (Sa) 


Inserting this value of s into (32) and using the assumptions which were 
stated for the derivation of equation (36), we find 


cr (IW) /a+¥) 
| PACA (37) 


an equation which might be used for the very simplest of cases. 
The avoidance function which we have been using to test for the pres- 
ence or absence of an avoidance reaction is also useful for other purposes. 


Vin (Eo ee, =| 
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One of these concerns itself with finding the steady-state solution of the 
number of organisms at each point within the container. As mentioned 
previously, one of the major difficulties in solving the steady-state cqua- 
tion has been the sudden change of the stimulus within less than a few 
mean free paths. A suggestion has been offered by W. H. Thorpe e¢ ai. 
(1947) on how to get around this difficulty. Their method is to use the 
fact that in the steady state the number of organisms crossing the bound- 
ary is the same for both directions. Then, from an equation analogous to 
our (25), with (cr)o(*) = é(«)79(«), they find, if L~/2 represents a point 
immediately to the left of the boundary and L+/2 represents a point 
immediately to the right of the boundary 


ne). 
= —. 
sG) G) 
Using the above result and (4), we then have the two necessary boundary 
conditions to evaluate the two constants of (6). However the above meth- 
od will not be used in this paper for three reasons. One is that it is not 
applicable if there is an avoidance reaction at the boundary. The second, 
and more basic, reason is that (6) does not actually hold a¢ the boundary, 
but only a few mean free paths from the boundary (as mentioned in the 
Introduction). Thus, since N(Z~/2) and N(Lt+/2) in the above equation 
represent the actual values of the number of organisms at the boundary, 
and not the values which are found from (6), (38) cannot be used as a 
boundary condition. The third reason is that the above equation was de- 
rived for the case of a boundary of zero thickness. In any experiment the 
actual boundary between two regions may be as wide as one or two mean 
free paths. These difficulties may be eliminated by the use of the avoid- 
ance function. From (26) we see that the number of organisms at two 
points on either side of the boundary are related through Yel Ewe) 
these points being chosen to be many mean free paths from the boundary. 
Hence we are now able to solve the differential equation completely and 
find the value of W(X) at all points. This value of V(X) may be compared 
to experimentally derived values of V(X), and if the two values disagree, 
further inquiries may be made into the cause of disagreement. 

Another use for (£;, &) may be found in quickly and easily ascertaining 
whether there is a taxic component to the movement of an organism for 
the case where the mean free path is constant within each region, though 
differing between regions. This may be seen if we note that for the condi- 


(38) 
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tion of a constant mean free path and for no taxic component g(x) 
= r(x) = 0 in each region. Thus from (5) we see that N(X) is constant 
within each region, as would be expected from physical considerations. 
Then, from (26), we see that v(é, &) equals one if £ and & are both in the 
same region, which is also physically obvious. However, if the mean free 
path is constant in each region, while the taxic compound is not zero at 
all points within a region, we see, by reasoning similar to that used above, 
that if £ and & both lie in one region and the taxic component is not zero 
for the sub-region lying between & and &, v(&, &) does not equal one 
except for some exceedingly unusual combinations of £1, &, and the taxic 
components. In general, we may feel fairly confident in saying that if 
v(£, ) within a region is experimentally measured and found to be one, 
there is no taxic component to the organism’s movement in the sub- 
region between & and &. On the contrary, if v(&, &) does not equal one 
then we may say that there is a taxic component to the organism’s motion 
in the sub-region between & and &. The values of e,(X) and (c7).(X) 
might be calculated by use of equations analogous to (29) and (30), but 
except for some very unusual cases only the function g(a) may be found 
and only very rarely is it possible to separate out the values of the e.(X) 
and (cr).(X). However, this method does provide an independent check 
on the values of those parameters which are calculated by the methods 
outlined in the previous section. 

Inactive organisms. If some of the organisms in an experiment are in- 
active, then the fourth assumption of I, that the time spent in turning is 
negligible, is violated and a different approach must be used. It may be 
noted that inactivity only signifies that the organism’s speed is zero. This 
is obviously true if the organisms are at rest, and it is also effectively true 
if the organisms move about in small circles, i.e., exhibit ‘“circus-move- 
ments” or exhibit any other movement which has no translocatory effect. 
This has been called “virtual inactivity” by Gunn and Pielou (1940) and 
generalized by Bentley et al. (1941). In agreement with these authors, we 
shall consider organisms which exhibit virtual inactivity as being inactive. 
Inactivity is completely analogous to what, in a study of diffusion of 
particles in a solution, may be called reversible binding of the solute 
particles to a fixed matrix. 

Only the steady state will be considered here. Since in the steady state 
there is an equilibrium between the number of active and inactive organ- 
isms in each region, the ratio of the number of organisms which are active 
to those which are inactive at each point is a constant with respect to time. 
It will be assumed here that the organism enters the inactive state only 


ORIENTATION OF ORGANISMS 451 


at the time when it turns, i.e., the probability distribution of 7’s is the 
same for those organisms entering the inactive state as it is for all the 
active organisms. Then, in the steady state, each time an organisin turns, 
the probability that it enters the inactive state is equal to the probability 
that another organism at the point of turning leaves the inactive state. If 
it is further assumed that the probability that an organism enters the 
inactive state is independent of its past history, then the number of organ- 
isms which leave a given region per unit time is determined by the same 
methods as used in I, provided the number of organisms in the region is 
equated with the number of active organisms in the region. Thus we may 
consider the situation as if there were two populations: an active one in 
which equilibrium is established in the manner derived previously, and an 
inactive population which is in equilibrium with the active population at 
each point. 

Let the subscripts , and ; refer to the active organisms and the inactive 
organisms respectively, and let o;(X) and o.(X) represent the fraction of 
inactive and active organisms respectively at the point X. Then 


o;=1—4,(X) 


N:(X) (39) 
= Oe O.eP: 


If, as before, we let V(X) represent the total number of organisms at the 
point X, then 
N(X) = Ni(X) + Na(X) 


Re NS, (40) 
~og(X) 1—4;(X) 


If, as before, we assume that the stimulus varies only in the x-direction 
and is constant in the y- and z-directions, then (4), (5), and (26) may be 
used to find NV,(X). In the derivation of equation (26), only active 
organisms were considered. Hence, if we include inactivity, (26) becomes 


(é é ) oT, Na (&y) 0(&:) (cr) 0(éy) 
rostt 8) WN, (Et) tol £7). Er) 0€s) As 
a N (&y) to (és) (cr) o(éy) da (Ey) ) 
~ N (€r) 7° (Ey) (7) 0 (Et) oa (Es) 


We therefore see that (35) remains valid even for the case in which some 


of the organisms are inactive. ; 
To consider some further characteristics of o,(X), the following rela- 
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tion is derived. If g(X) is the probability that an active organism enters 
the inactive state each time it turns at the point X, then by an argument 
similar to that used in I, the average number of organisms that go from the 
active to the inactive state per unit time at the point X is 
ee. ¢ X)= 
Na (X) £ QO =e = X)° 
Thus if «(X) is the average time that an organism remains inactive, then 
the total number of inactive organisms at the point X is 


Ni(X) = Na (X) g(X)= Ty = ; u(X). 
Inserting the above into (39), we find 
1 
g(X)u(X)=— yy 
Gey) = abt 3 (42) 


a g(X) 4 (=e 
If either u(X) or g(X) is zero, then o;(X) is zero; if either u(X) or g(X) is 
infinite, then o;(X) is one. These results are what we would expect from 
physical considerations. However, if #(X) is infinite and g(X) is zero, 
then o;(X) is indeterminate. That is, physically we have a situation in 
which a certain number of the organisms are “permanently”’ immobile 
and their positions are determined by chance. For this case, the problem 
may be treated by ignoring the permanently inactive organisms and con- 
sidering the active organisms as if they were the whole population. 

Approximation method solution. In previous sections of this paper it has 
been shown how the equations derived in I can be applied to the problem 
of the orientation of organisms and these equations have then been solved 
for the steady-state distribution of the organisms. Since all parameters in 
the equation are arbitrary functions of the coordinates and of time, it is 
not possible to find the general solution of the equations for the non- 
steady state. However, it would be useful to obtain some idea of what the 
actual solution would be like. To do this the approximation method as 
outlined by N. Rashevsky (1948) will be used. In this method, instead of 
attempting to find the exact solution of the equations for all points in the 
container at all times, we will focus our attention on the average values of 
the concentrations in the two regions of the container, and attempt to 
find how these average values change with time. 

We shall consider a container divided into two equal regions, which we 
may denote as regions ¢ and ¥ as before. Since, as stated above, we are 
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not interested in exact solutions but are only interested in obtaining cer- 
tain crude but, perhaps, useful results, we shall take the average values of 
the various parameters in each region and treat these parameters as if 
they were constant within each region. Furthermore, we make the follow- 
ing assumptions: There is no variation in the y- and z-directions, the total 
number of organisms in the container remains constant, there is no taxic 
component to the organism’s movement [e,(X, #) = (¢r)2(X, #) = 0], the 
coefficient of persistence of direction is constant throughout the container, 
there is no avoidance reaction, (€2r2)o(X, #) = (cr)?(X, 2), and all of the 
values of the parameters are independent of time. 

From the assumptions above, it can be seen that in the steady state the 
distribution of organisms in each region will be uniform. We will not at- 
tempt to find the distribution of organisms within each region, but will 
only be interested in the average number within the regions. Therefore, 
we will treat the problem as if there is always a uniform distribution of 
organisms within each region. As in I, let I',(x, ¢) represent the net num- 
ber of organisms that, per unit of time about the time /, cross a plane 
of unit area perpendicular to the x-axis at the point X. We may then write 
(2) as 


¥ 
—¥ 


Ne, t) t) 
To (x) 


1(1+y) @ (er) 0(*) 
a ieee oe Sq| Ne OS =p ad ie 


Px (2,1) =24>¥ (er) 0(2) (er) 0(#) 1 


which, upon rearrangement, yields 


fexie (eNCl 4y)0 


Ceri te) 
Prieta so = Aerie Aa) | 


sal ¥ 4) T9(X) 


Letting H be one for the one-dimensional case, the width of the container 
for the two-dimensional case, or the cross-sectional area of the container 
for the three-dimensional case, we have approximately 


dN (¥, 1?) H(1+y¥) PTa intone ieee) 
eat in) a3 2. + Im(1-y) 2 
(cr ib 2/ (14) (y) (63 ee orb 2/(14y) (¢) (43) 


x 


iis 
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We shall assume further that the organisms move only under the effect 
of pure ortho- or pure klino-kinesis. If we let (@7)o represent the average 
value of (€272)9(X)/70(X), i-e., 


Af eo tee 


COT) e— Flea rae 


and if | (cr)o(v) — (cr)o(S)| is less than, say, .2 (cr)o(f), we then have 
to a first approximation that the solution of (43) is 


N (y,; #) 
4(1+y) (cr)o (44) 
n(i-y) FL 


= [N(y,0)-N(y @)] exp }— t+ N(y, ©), 
where NV(y, 0) is determined from the initial experimental distribution and 


Na Gere” (2) 


HL to (£) 
N (yy, ©) == =e ‘ (45) 
70(f) 70 (¥) 


It should be noted that the value for V(y, ~) is the same as that which 
is found from the exact steady-state solution of (2). If equations analo- 
gous to (44) and (45) for V(¢, #) had been derived, we would find com- 
pletely similar results to those above except for the interchange of ¢ and y 
in all of the parameters of (44) and (45). If we then take the ratio of 
N ({, ~) to N(y, ©) we find it to be 


N(f,@) _ [ (cr) o(y) | ede 3 (46) 
N (y, ©) (er)o(¥) Toy) 


Since the distribution of organisms in each region is taken to be uniform, 
the above equation also gives the ratio of the total number of organisms in 
each region. 

It can be seen that for the above conditions the expression for the dis- 
tribution of organisms in each region contains an exponentially decaying 
transient term, the time constant, ¢*, for the decay being 


ee n(1—y)L? 
4(1+y) (c?r)o 


(47) 


For any experimental situation it is highly improbable that all of the 
assumptions will be fulfilled, and hence equation (44) will not hold for 
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these conditions. However, as stated before, it may be assumed that (44) 
will hold very approximately for most experimental situations, i.e., the 
distribution will change in a manner which may at first glance be con- 
sidered to be exponential. Also, it is to be expected that the form of the 
time constant will not be very different from that given in (47) even 
though the assumptions leading to this expression are not fulfilled. Thus 
we may expect it will take a longer time for the distribution of organisms 
to reach equilibrium in a long container than in a short one, etc. Also, ex- 
pression (47) may give the experimenter some idea of the order of magni- 
tude of the time constant for a particular experiment, and thus may 
indicate very roughly how long he will have to wait before an equilibrium 
distribution is established. 

Discussion of experiments. To illustrate some applications of the equa- 
tions which have been developed in this paper, and also to test their 
validity and applicability, some examples of animal orientation taken 
from the experimental literature will be analyzed. However, only a small 
segment of the existing literature will be treated, due to the fact that there 
are only a few experiments with a sufficient amount of data so that the 
equations may be intelligently applied. 

1. Experiment of Uliyott (1936). In the work of Ullyott (loc. cit.) it was 
clearly shown that organisms may congregate in a preferred region by 
means of non-directed movement. In this case Ullyott showed that the 
organism’s ‘‘rate of change of direction” (r.c.d.) changed when it entered 
a region of increased illumination and that the r.c.d. accommodates to a 
basic level as the animal becomes adapted to the new environment. This 
type of behavior was given the name “klino-kinesis” by Gunn et al. (1937) 
and later the name was changed to “klino-kinesis with accommodation”’ 
by Ewer and Bursell (1950). Ullyott further showed in his paper that this 
type of behavior should yield congregation for steep gradients of illumina- 
tion and not for mild gradients. His demonstration of this theory was 
based on the assumption that a change of r.c.d. implies a change in the 
rate of turning, i.e., a change of 7o(X, #) in the notation of this paper. 
However, this assumption, which seems to be quite plausible and obvious, 
is actually false, as will be demonstrated below. Although this may in- 
validate the argument used by Ullyott, it need not invalidate his con- 
clusions. | 

To see what information may be obtained from Ullyott’s paper, it is 
necessary to find values for the parameters of the equations given in this 
paper and also tosee if the assumptions involved in deriving the equations 
are valid. With this in mind, Figure 6 of Ullyott’s paper, the only piece 


456 CLIFFORD S. PATLAK 


of data available (the track of the organism in a container in which the 
illumination was suddenly increased), was used. While many more tracks 
than this should be examined before a definite conclusion can be stated, 
it may be assumed that this track is representative, especially since the 
results derived from it are so striking and self-consistent. 

The tracks were examined under a hand lens. The number of turns 
were counted for each five-minute period (the total track was for a period 
of a little more than 25 minutes) and the angle each separate segment of 
the track made with a fixed direction was measured to the nearest five 
degrees. From this type of measurement it is obvious that, in the nomen- 
clature of I, the angle so measured would be the angle @ of the segment 
and would, similarly, be the angle 6, for the next segment. By taking the 
difference between consecutive angles, the angle of turn (a) could be found. 
Since, for this example, the light is uniform throughout the container, we 
may assume that there is no taxic component to the motion and hence 
€x(X,t) = e,(X,t) = (cr)2(X, t) = (7), (X, #) = 0. It is obvious from the 
experimental set-up that a turn of zero degrees (which for all practical 
purposes we may take to mean a turn of less than two degrees) cannot be 
observed. 

The following results were found. The total number of turns for the 
25-minute period was 208. An examination of the angles made with the 
fixed direction (6; or #2) shows that there is a fairly uniform distribution of 
angles from 0° to 360° and that the average angle of turn is 191°. For a uni- 
form distribution of angles between 0° and 360°, the standard deviation of 
the angles should be 104° (the actual standard deviation for the distribu- 
tion of the angle @ is found to be 115°) and thus the standard deviation of 
the mean would be 104°/ W208 which is 7.2°. Thus the value of 191 is less 
than two standard deviations of the mean from the expected value of 180° 
for a uniform distribution. A similar analysis was carried out for the values 
in each of the five-minute periods and the values were found to agree with 
those expected from a uniform distribution in a way analogous to those 
above. 

A numerical analysis of the angle a, as well as an inspection of a graph 
of its values, shows that the values of a are symmetrical about the origin, 
which indicates that the organism has no right- or left-hand preference 
in its turnings. Thus, as may be expected, an evaluation of sin a for each 
of the five-minute periods indicates that this average differs from 0 by 
less than 15 standard deviations of the mean, which would also be in agree- 


ment with the assumption of no taxic component to the motion, as indi- 
cated in (15). 
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The number of turns in each five-minute period varied between 40 and 
42, thus showing that 7)(X, ¢) was constant during the change in the 
r.c.d. In order to evaluate the coefficient of persistence of direction, both 
€oSa [eq. (15)] and ¥ cos 6; cos 62/E cos? # [eq. (16)] were computed. Com- 
puting the average values plus or minus one standard deviation of the 
mean, we find for the €6sa for each five-minute period the values .39 + 
09, .59 + .06, .61 + .05, .60 + .07, and .70 + .04, with an average value 
of .58 + .03. For = cos 4, cos 62/= cos? 6; we have for each five-minute 
period the values .37, .61, .59, .44, and .79, with an average value of .56. 
We thus see that the two methods yield results which agree as well as 
could be expected. 

There is no track given to illustrate how the organism moves under 
the influence of weak stimuli, but it can be inferred that it is not appreci- 
ably different from that of the last five-minute period of the track dis- 
cussed above. If this is so, then we see that the organism upon entering an 
area of increased light intensity will not change its rate of turning but will 
commence making larger turns. It will then accommodate by decreasing the 
size of its turns, i.e., the coefficient of persistence of direction increases. 

Let us now see what results the equations which have been developed 
in this paper would yield for this case. We may assume that there is no 
taxic component to the organism’s motion; that, from the above discus- 
sion and from Ullyott’s statements on the speed, the parameters (Cr) 0(X, #), 
(c272)9(X, #), and 7,(X, #) are constant; that, as Ullyott states, there is 
no avoidance reaction; and that only the coefficient of persistence of 
direction varies with position. If we further assume that there is no ac- 
commodation, then (5) reduces to V(X) = Constant, which is not what 
is experimentally observed. Since the above assumptions seem valid 
except for the one involving non-accommodation, and since the other 
assumptions leading to (5) are fulfilled as shown above, we infer, in 
agreement with Ullyott, that the assumption of accommodation is 
essential for the congregation of the organisms in this case. 

It may be of interest to point out here that even if the organisms orient 
themselves only by means of pure klino-kinesis with no accommodation, 
i.e., To(X) is a function of position and the coefficient of persistence of 
direction is a constant, it is possible to have the animals congregating in 
different parts of the container for a given 7o(X) distribution, depending 
on certain other factors. This point may best be brought out by an illus- 
trative example. 

Consider an organism moving about under pure klino-kinesis with no 
accommodation or taxic component to its motion. Then (8) would apply. 
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Let us further assume that 72(X) = /73(X), where / is a constant. Then 
for this case (5) reduces to 


N(X) = Krphi@D/+4¥ 21 (X) . 


Thus if we assume that y > 0, then since y < 1, the organisms will con- 
gregate in the region of lowest 7o(X), i.e., in the region in which the organ- 
isms turn most often. 

However, let us now attempt to take into account experimental condi- 
tions of a steep gradient of the stimulus which may cause a rapid change 
of To(X) along the container. This may be done (if we do not wish to treat 
the problem by use of the avoidance coefficient) by the use of an artifact. 
This assumes that the organism’s 7 is changed at any moment from the 7 
it would have possessed if the stimulus was constant over its walk between 
turns [this change being in a manner analogous to the change of T9(X)]. 
Since the organism will turn after it has traveled for a length of time equal 
to the 7 the organism possesses when it turns, then we may crudely take 
for 7-(X) the value of (é7)o(X)d70(X)/dx instead of the value of zero as 
we did before. For this case, we then find that (5) reduces to 


N(X) = Kry@-)0-)/+¥.2-4)] (XX), 


Hence if we assume that y > 0 and that # < 2, then the organisms will 
congregate in the region of highest 7)(X), i.e., in the region where the 
organisms turn least often, contrary to the former case. If it is assumed 
that the r.c.d. of Ullyott is proportional to 1/t9(X) (which, as seen above, 
is actually false), then we see that the second of the two above equations 
would yield qualitatively the same result as that found experimentally. 

It should be mentioned at this point that in the discussion of the effect 
of the stimulus upon the organism, a slightly different viewpoint than 
that of Fraenkel and Gunn (1940) may be adopted. If we conceive that 
there is something, which need not be specified further, being produced 
in the organism which causes it to turn when a threshold value of this 
“unknown” is reached, then Fraenkel and Gunn (Joc. cit., pp. 51-53) would 
assume that the threshold is fixed while the rate of production of the 
“unknown” is governed by the intensity of the stimulus. However, it 
may be assumed that the rate of production is constant while the threshold 
is determined by the stimulus. To illustrate this distinction, consider an 
organism orienting under pure klino-kinesis with no accommodation, 
moving at constant speed, and traveling in a container where the stimulus 
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increases linearly with distance. Let us also consider two points, & and 
£2,such that if the organism turns at point & and travels toward point 
£, it will on the average next turn at point &. Then, if the organism 
turns at point & and travels toward point &, it will on the average next 
turn at point £, according to the hypothesis of Fraenkel and Gunn; 
whereas it will not on the average next turn at point &, according to the 
other hypotheses. The equations of the previous two paragraphs illus- 
trate this point, the latter equations being an illustration of Fraenkel 
and Gunn’s hypothesis, the former of the other hypotheses. More experi- 
mental work must be done to decide which, if either, of the two theories 
is correct. 

2. Experiment of N. Waloff (1941). Fairly complete data suitable for 
the equations which have been developed here can be found in Waloff’s 
paper, part of which (experiments on the organism Porcellio Scaber) will 
be used here. In the experiments the container used was basically the same 
as the ring container of Gunn and Kennedy (1936) except that it was oval 
in shape and a total length of 583 inches. The stimulus used was humidity. 
Two types of experiments were performed, those in which the humidity 
was kept constant throughout the container, and those in which two 
different humidities were established at opposite ends of the container and 
equilibrium values of humidity were allowed to be established throughout. 
This type of container was used in place of the usual alternative chamber 
because it was observed that the organism exhibited a large amount of 
“thigmokinesis,” i.e., the animals walked with one antenna pressed 
against the side of a wall for over 80% of the time. Hence by use of the 
ring type of container, Waloff took advantage of this behavior pattern 
and the behavior reduces itself to that of one-dimensional orientation. As 
we require, only one organism at a time was used in each experiment. 

First it may be noted that the values of various parameters such as 
speed, activity, etc. for the case of the organism in a container with uni- 
form stimulus is mot the same as for the organism in the container with 
non-uniform stimulus. For example, if the humidity in the ring type con- 
tainer varies from 75 to 95% (these being the extreme values), then it is 
found that the organisms are always active on the dryer half, while only 
50% are active on the wetter portion. However, if the humidity is made 
uniform, then it is found that the percentage of active time would vary 
from 90% for the relative humidity (R.H.) of 75% to 62% for the R.H. of 
95%, which is different than the observed values in the former case. For 
the case in which the humidity varies from 75 to 95% the average humid- 
ity in the dryer half is 80%, and, as stated above, all of the organisms are 
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active. However, for the case in which the humidity varies from 65 to 85% 
the average humidity in the wetter half is also 80%, but in this case only 
63% of the organisms are active. The most probable explanation of these 
discrepancies is that the organism is adapting to the total situation. If it 
accommodates to the entire situation within the first ten minutes (the 
time allotted by Waloff for equilibrium to be reached) then the equations 
of this paper may be used, providing that values of the parameters are 
computed from the data found by observing the organisms in the experi- 
mental situation. But if the organism adapts to each humidity in each 
separate portion, then the equations of this paper will not be valid. 
From graph 2 of Waloff’s paper we may infer that accommodation is 
very slow, taking about one hour to complete. Due to Waloff’s lack of 
comment on possible accommodation to the new situation each time the 
organism enters a new portion of the container, we will assume that the 
organism adapts to the entire situation, but perhaps slowly. If this adapta- 
tion is slow enough, we may assume that we have a series of equilibrium 
values and we are therefore taking the average values of the equilibrium 
distributions. This point will be discussed further in the next section. 
The data as given in the article are close to that required for applica- 
tion to the appropriate equations of this paper, but are not exactly what 
is wanted. Also, in some cases it is not precisely clear what measurements 
the values are for. Hence certain assumptions and simplifications must be 
introduced in order to use these data. The first point which must be tested 
is the validity of the values of the “average number of minutes between 
turnings.”’ The question that arises here is whether this value is 79(X) or 
T9(X)/oa(X), ie., whether Waloff measured her 7 for all organisms or just 
for those organisms which were always active. To determine this we will 
consider the ratio \(X)/éo(X)r(X), where A(X) is the average distance 
between turnings, ¢(X) is the average speed, and 7(X) is the average 
number of minutes between turnings. For the ten sets of values we find 
that the ratio varies from .82 to 1.13 with no correlation found between 
these values and the experimental values of the “percentage of time at 
rest” [o;(X)]. We thus see that 7(X) is the same as To(X), i.e., the meas- 
urements were presumably made for organisms which were always active. 
The expected value of the above ratio should be one. The discrepancy of 
as much as 20% may be explained as due to experimental errors (which 
should only constitute a small fraction of the error), to the fact that 
(¢7)0(X) probably does not equal é(X)ro(X ), to the fact that accommo- 
dation may have been present so that é(X) and 7(X) are constantly 
changing during the course of the experiment, and to the fact that we are 
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taking the average values of the parameters for a range of humidity in 
each half of the container. 

Another point that should be noted is that o;(X) is measured by count- 
ing only those organisms which are inactive for a full 30 seconds. Thus if 
an organism is inactive for only 29 seconds, we may infer that the author 
would have counted this organism as active for the full 30 second period. 
This probably did not introduce much of an error since in most cases if 
the organism rests it will remain inactive for a long period of time. How- 
ever, a more serious error is introduced due to Walofi’s counting organ- 
isms in “‘virtual inactivity” as active organisms. Both of these factors 
show that the values for o;(X) which are given in her paper are less than 
or equal to the true values of o;(X) for the organism. 

In the experiment discussed above, the organism was observed either to 
make a complete 360° turn, a 180° turn, or to turn 90° and then turn back 
to the direction in which it was traveling before it turned, this being 
equivalent to a 0° turn. The number of each type of turn is not given in 
the data, but the total amount of turning, where the 0° turn is counted as 
1, the 180° turn as 2, and the 360° turn as 4 is recorded for the total num- 
ber of organisms which were involved in each experiment. We may also 
establish the total number of times an organism turns since we know how 
long each experiment lasted (50 minutes), the average time between 
turns, the fraction of time that the organism spends in each half of the 
container, and the fraction of time that the organism is active. Thus if x 
represents a particular half of the container and if («) is the number of 
0° turns, (x) is the number of 180° turns, and 73(x) is the number of 
360° turns that an organism makes in the x-half during the time of the 
experiment, then for any experiment 


(“total no. of turns” in x) ‘2 
(no. of organisms in the container) 


M(x) + 2n2(x%) +43 (x) = 
(48) 
50ca(«) (fraction of total time spent in Nite - 


my (*) + m2 (x) + 3(%) = 1 (x) 


We therefore have only two equations for three unknowns. Because of. 
this, a further assumption must be made. This will be that n3(x) =.0. This 
assumption is valid for the higher humidities (over 85% R.H.) but not. 
for the lower ones; however, it is made so that some solution may be found 
for the above equation. Since the coefficient of persistence of direction was 
defined in the one-dimensional case so that the probability of an organism 
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continuing in the same direction was 3[1 + ¥(x)], which in the above case 
is equal to m(x)/[m(x) + mo(x)], then 


n(x) — Ne (x) 


VNt): = ea eae ae 


Hence inserting the values of m(x) and m(x) found from equation (48) 
into the above equation [with 13;(x) = 0] we find 


7 (x) (total no. of turns in x) 


~ 250, (x) (fraction of total time spent in’ 
x) + (no. of organisms in the container) 


v(x) =3 (49) 


It may be noted that the number of organisms in each half must be used 
to calculate y(x) and that this value of ¥(x) will then be used in turn to 
find the number of organisms in each half. However the value of ¥() is 
not too critical in computing this number of organisms [see eq. (46)] and 
in an experiment in which m(), m2(a), and 23(*) would be given directly 
as part of the data the number of organisms in each half would not be 
needed to compute any of the parameters of the equation. Upon evaluat- 


TABLE I 

PERCENTAGE OF ORGANISMS IN 
RELATIVE THE WeT HALF 
Hoummity Experimental Theoretical 
75-95% 87.2 79 
65-85% 73.8 80 
55-75% 70.8 66 
20-73% 74.5 60 
42-60% 67.7 58 


ing ¥(x) for 10 different cases, i.e., for the two halves of the container for 
five experiments, it is found that in three of the cases for the damper half 
of the container the value of ¥/(x) is slightly greater than one. However, 
¥(«) by its definition cannot be greater than one. This error may be due to 
reasons similar to those mentioned previously in the consideration of the 
ratio A(X)/é)(X)7(X). For these cases ¥(x) is taken to be equal to one. 
If we now take the values of the various parameters which we have com- 
puted or taken directly from the data and insert them into (46) and then 
from the derived ratio of the numbers in each half of the container find 
the fraction in the damp half we have Table I. 

The following interesting point about the results may be noted. In one 
case the experimental value is smaller than the theoretical one, whereas in 
the remaining four cases it is larger. This, then, would tend to indicate 
that there is no avoidance reaction present, or if there is it is very slight. 
This might be expected from the fact that there is a maximum of 50% 
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R.H. difference, this being over a length of 29 inches, i.e., there is a very 
small gradient of R.H. The most likely situation in which there might be 
an avoidance reaction is in the case of the 20-73% range, which would fit 
in quite well with the fact that for this case there is the largest discrep- 
ancy between theory and experiment. There are two factors which tend 
to minimize the importance of the deviations between theory and experi- 
ment. One of these is the 20% error in the data, mentioned previously, 
and the fact that the above differences are less than or about 20% in size. 
The other is the fact that since the humidities are ot uniform in the two 
halves of the container, (46) is only approximate. Considering all of the 
factors mentioned previously, however, agreement between experiment 
and theory is satisfactory and does seem to indicate that there is no taxic 
component or avoidance reaction in the orientation of the organism. 

Another point in Waloff’s paper which may be discussed concerns ¢*. 
In order to obtain some general idea of its order of magnitude the fol- 
lowing values (representative of those in her paper) of the parameters are 
used: Y = .6,7) = 1.3 minutes, LZ = 29 inches, and é = 6 inches/minute. 
We take LZ to be 29 inches instead of the full length of the chamber (58 
inches) because in the derivation of the time constant a simple alternative 
chamber was used, implying reflection at the ends of the chamber. How- 
ever, in the ring container there are no ends, but the point at which the 
stimulus is a maximum or minimum for each half (i.e., at the center point 
of each half) may be considered to be equivalent to the end in the case of 
the alternative chamber. That is, the two “arms” of each half are con- 
sidered to be equivalent to each other. Thus Z is actually one-half of the 
true total length of the container. If we further assume that the probability 
distributions of the c’s and r’s are narrow, we then find that (@7)o = Gro. 
Inserting the above values into (47) with = 1, since we are dealing 
with a one-dimensional case, we see that ¢*~ 1 minute. Thus the value 
of 10 minutes which Waloff used in allowing her organisms to reach 
equilibrium seems to be sufficient time. However, since the value of;t* 
above is only very approximate, and since at least 3¢* are required to 
reach equilibrium, then it may be advisable that a somewhat longer time 
should have been allowed for the organisms to reach the steady state. 
Also, the time it takes for the organisms to accommodate to the experi- 
mental situation should be taken into account. 

3. Experiment of J. S. Kennedy (1937). We shall now discuss some ex- 
periments reported in the literature in which, contrary to Waloff’s work, 
there is evidence of an avoidance reaction. The first paper which we will 
consider is that of J.S. Kennedy. Figure 4 on page 194 of his article will be 
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analyzed. In this figure the path of the organism (the African Migratory 
Locust), which is in an alternative chamber, is shown. We will consider 
only that portion of the path which is not along the edge of the container 
and will draw lines parallel to and one-half inch from the central line be- 
tween the dry and moist sides. We will then count the number of paths 
which cross each of these lines heading toward the center and follow these 
paths so as to find how many of them cross the other line before crossing 
the first line again. Doing this we find we have four paths crossing the 
line in the dry side of the chamber and traveling toward the moist side, 
only one of which crosses the line in the moist side before crossing the dry 
side line again. For the line on the moist side we find three paths which 
cross this line traveling toward the dry side, all three of which cross the 
dry side line before crossing the moist side line again. If d represents the 
dry side and m the moist side, we then have n(d, m) = } and n(m, d) = 1. 
Thus from (17) we have »,,(d, m) = }. If we assume that these data are 
representative of the organism’s path, we may say that the organism 
avoids the moist side, which is in agreement with what may be observed 
by further examination of the path and of the organism’s movements. 
However it should be noted that although nothing was found by this 
analysis that was not known before, the method yields a quantitative 
statement of how strong the avoidance reaction is. 

4, Experiment of Gunn (1937). The next paper which will be examined 
is that of Gunn on the behavior of the organism Porcellio Scaber, the 
stimulus being humidity. The experiments were conducted in an alterna- 
tive chamber, which may explain why the results are different from those 
of Waloff (loc. cit.), who worked with the same organism. Since no data 
are given for the behavior of the organism in the alternative chamber, we 
are therefore forced to use the assumption that the behavior of the organ- 
ism in each half of the alternative chamber is the same as its behavior in 
the same humidity in the constant humidity chamber. Gunn states that 
congregation of organisms is solely due to the different degrees of inactiv- 
ity of the organism in the different humidities. From Figure 4 of his paper 
and the discussion in the text we find 


o;(%) =.19 0% <R.H. < 30% 
=.0102(R.H.)—.115 30%<RH.<70% 
= .60 70% <R.H. < 100%. 


To find the avoidance function the following simple relation will be used. 
Since the only difference in the organism’s behavior between the two 


ORIENTATION OF ORGANISMS 465 


halves of the chamber is due to the difference of activities, then the theo- 
retical value of N.(X) for the dry and the moist halves should be the same. 
Hence using (40) we see that 


N (m) © oF) 
Poe ae 28) 
Inserting (50) into (35) it can then be seen that 
* oe N (m) 1 —o;(m) 
via (d, m) = | ‘elie ay I some 


Using the above equation we find that in this experiment, for the mean 
R.H. of less than 85%, the avoidance function is always greater than one 
(ranging from 1.25 to 5.15), i.e., the organism always avoids the dryer 
side, which agrees with Gunn’s statement on pages 179-80. However, the 
value of the avoidance function various tremendously, with no obvious 
connection to either the mean R.H. or the range of R.H. used in the ex- 
periment. This may be due to the fact that there are wide experimental 
fluctuations and also large differences between organisms, as shown in his 
Figure 4. It should also be noted that many of the statements made by 
Gunn are at variance with those of Waloff. For the mean R.H. greater 
than 85% the avoidance function varies from .75 to 1.25. Due to the rea- 
sons mentioned above, we may assume that the avoidance function for 
these R.H. ranges actually should be one. Hence we see that in agreement 
with the statement of Gunn the organisms do not distinguish very well be- 
tween two humidities, if both are high. However, contrary to his conclu- 
sions, it appears that the organism can distinguish between two low 
humidities. 

5. Experiment of Bursell and Ewer (1950). Another paper which we 
shall examine is that of Bursell and Ewer. This deals with the behavior 
of the Wood-Mason in an alternative chamber, humidity being the stimu- 
lus. Since the exact data which would be needed to find the parameters of 
the equations are not given, rather arbitrary and simplifying assumptions 
will be made here. There are 1) that the distributions of the c’s and 7’s 
are such that (c272)o(X) = @2(X)73?(X); 2) that ¥(X) is constant through- 
out the chamber (we will choose y to be .6); 3) that the authors did not 
include “virtual inactive” organisms among the number of active organ- 
isms, so that Figure 4 of their paper yields o,(X) directly; 4) that the 
humidity varies linearly along the container so that the average value of 
the humidity in each half of the container is found by averaging the 
extreme value of the humidity in that half with the average value of the 
humidity in the container; and 5) that the behavior of the organism in 
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each half of the alternative chamber is the same as its behavior in the 
same humidity in the constant humidity chamber. Another assumption 
made is that there is no adaptation. This fact is experimentally confirmed 
by the authors and therefore may be considered valid. 

Since we assumed y to be a constant, therefore (€7)9(X) is proportional 
to 1/r.c.d., measured in units of degrees/cm. Since we are finding ratios 
of concentrations it will not be necessary to find the value of the constant 
of proportionality for this relation. Then, from (40) and (46), after some 
simple rearrangements, we have 

[Aw] -2@ co (d) Aisa rama oes 
N(d) Jin Ga(m) co(m) L (rec.d.) (d) 


Inserting the above equation into (35) we obtain 


ria (d, m) =| 


| oa(m) Co(m) [ (r.c.d.) (d) (itachi (53) 


N (d) Jex Ga(d) co(d) L (r.c.d.) (m) 


Inserting the values of the various parameters found from the data given 
in the text into the above equation, we find that some interesting con- 
clusions may be drawn. The first is that 


[pases (d) pean? 
(r.c.d.) (m) 


varies only from 1.00 to .93. Thus this factor is almost negligible in its 
effect on the final results. Next for the three cases in which the mean R.H. 
was greater than 94%, the avoidance function was less than one, a result 
which may not be expected However, due to the large number of as- 
sumptions made and the fact that we are dealing with a situation in 
which the R.H. constantly changes with position in the chamber so that 
the equations used are only approximate, this result may not be correct. 
The final conclusion is that for all remaining cases in this experiment, the 
avoidance function is greater than 1 (ranging from 1.53 to ©). Upon com- 
paring the v2,(d, m) with the difference between the R.H. extremes, it is 
noted that, in general, if the two values for the average R.H. less than 
40% are ignored, the v2(d, m) increases slightly as this difference in- 
creases until a difference of about 15% is reached, at which point the 
vea(d, m) increases very rapidly until at a difference of 20% the 
vii(d, m) is greater than 15. These results are in qualitative agreement 
with the experimental results found by the authors (see Table 9 of their 
article). The two low R.H. cases have v%,,(d, m)’s of 1.6 and 3.5 respec- 
tively, which, although significant avoidance functions, are still very small 
compared to what would be expected from the discussion above. This may 
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indicate that the organisms do not distinguish the differences between 
R.H. very well when both R.H.’s are low, which is also in agreement with 
the predictions of the authors. 

Another point which should be noted in regard to Bursell and Ewer’s 
paper is that concerned with the time course for equilibrium to be estab- 
lished. In their paper the terms W and D are the same as the terms 
N(m)HL/2 and N(d)HL/2 used here. Thus using (44) and assuming that 
initially there is the same number of organisms in each half of the cham- 
ber, we see that the intensity of the reaction (I.R.) is given by 


LR. = 100 epee raul = eee gieth 
—[FF- Na, @) Jew wa, 0) | AE (54) 
= roof Vim eV NG) Ba — e-t/t*) 


= (final IL.R.) (1 — e-‘/**), 


If we then plot 100(1 — e~’’*) vs. ¢ with ¢* taken as 2 minutes, we find 
that the points of Figure 1 of their paper, in which experimental values of 
the percentage of final I.R. are plotted vs. #, fit the theoretical curve 
extremely well. 

In order to see if the value of ¢* which was found from experimental 
data agrees with that predicted from (47), the following assumptions 
would have to be made. From Table 1 of Bursell and Ewer’s paper we see 
that @& ~ 16 cm./minute and that the angular deviation is 8 degree/cm. 
If it is assumed that the organism makes on the average a 60° turn, then 
we may say that y is approximately equal to cos 60° which is .5. Since 
the angular deviation equals 8 degree/cm., which from the value of é 
above is equal to 128 degree/minute, therefore r = (60 degree/turn)/ 
(128 degree/minute) = .5 minute/turn. If we also assume that the prob- 
ability distributions of both the c’s and the 7’s are narrow, then (@7)o 
= 7, = (16 cm./minute)? X (.5 minute/turn) ~ 130 cm./minute turn. 
We will also assume that the chamber is 22 cm. long, as in the experi- 
ments of Gunn and Kennedy (loc. cit.), although Bursell and Ewer do not 
give any value for the chamber size in their paper. Therefore, inserting the 
above values into (47) with m = 2 since we have a two-dimensional case, 
we find * ~ .6 minutes. Although this is only 3 of the value found experi- 
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mentally (see above), it is at least of the same order of magnitude. The 
difference between theory and experiment is not necessarily as poor as it 
may appear since the actual values which should be used for the organ- 
isms in the alternative chamber were not known and many of these values 
(such as y, L, and 7) are, at best, estimates. 

6. Experiment of Gunn and Pielou (1940). The next paper which will 
be examined is that of Gunn and Pielou on the Mealworm Beetle in an 
alternative chamber, humidity being the stimulus. As seen from their 
Figure 3, the organisms accommodate to their environment. However, we 
shall first attempt to see what the equations developed in this paper will 
yield if we assume that there is no accommodation. The authors state 
that there is no ortho- or klino-kinesis in the organism’s motion. Hence the 
difference in time spent in the two halves of the chamber is due to the 
difference in activity of the organisms in the two halves, and, possibly, to 
an avoidance reaction. If there is no avoidance reaction, the ratio of the 
number of organisms in each half would be given by (46). 

Some assumption must first be made concerning the time spent by the 
organism close to the edge of the chamber. It is obvious that if an organ- 
ism travels around the container by moving along the edge, it is not mov- 
ing about randomly and hence the equations of this paper do not apply, 
i.e., the organism is effectively moving about as if it were in a one-dimen- 
sional situation while we are interested in its two-dimensional motion. For 
an organism of this type one may either perform an experiment similar to 
the one done by Waloff (Joc. cit.) or consider that when an organism is 
moving along the edge it is not part of the experiment, which is what will 
be done here. However, we must be careful with the latter approach for 
the following reasons: 1) if an organism moves along only a very small 
portion of the ledge, this should be counted as part of the experiment 
since the theory assumes that the organism will be reflected from the wall 
(which is what a slight movement along the wall may well be); 2) an 
organism may move in a path parallel to the wall even if the edge was not 
present; and 3) if an organism does move along the edge, it will not be at 
the same place or even necessarily in the same half as it was before it 
started to move when it leaves the edge to resume its motion in the 
chamber. This factor may not be important if the points at which the 
organism starts walking next to the edge and the points at which it leaves 
the edge are both random. i 

Using the above assumption of not counting the organism when it is 
moving along the edge, and counting the time spent in “virtual inactiv-: 
ity” as part of the inactive time, we find that [V(d)/N(m)]p, is equal to 
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16.7 in contrast to [V(d)/N(m)].., which is equal to 54.2. It should be 
noted that the data for o;(~) were recorded for the organism in actual ex- 
perimental conditions, which is what is desired in order to apply the re- 
sults of this paper. Using (35) we see that v%j(m, d) is equal to 3.24, 
indicating that the organism avoids the moist side, in agreement with the 
conclusions of the authors on the basis of a detailed analysis of the organ- 
ism’s path. If we ignore the above discussion with regard to organisms at 
the edge of the container and count the time spent close to the edge as 
part of the active time, we would find that [V(d)/N(m)],, is equal to 7.38 
and [NV(d)/N(m)]-x is equal to 23.7. These values would yield a v*,(m, d) 
of 3.21. However, if the time spent close to the edge had been considered 
as part of the inactive time then we would have found that [V(d)/N(m)|u, 
was equal to 7.35 and [V(d)/N(m)].. was equal to 23.7. These values 
would yield a v¢,;(m, d) of 3.24. The fact that all three values of v*,(m, d) 
are almost identical seems to be entirely fortuitous. 

As stated above, the previous computations are based on the premise 
that the organism does not accommodate to its surroundings. However, 
it should be noted that the experiments leading to Figure 3 of Gunn and 
Pielou’s paper, which indicated adaptation, were not carried out in a 
manner which would allow us to draw conclusions for the alternative 
chamber situation. The experiments were done by placing the organism 
into a constant humidity chamber and then activating it by pushing it 
with a pipe cleaner. This is obviously not the situation with which the 
animal is confronted in its motion about the alternative chamber. For the 
purpose of this paper the experiment should be done by placing the organ- 
ism in the experimental situation after it has become adapted to some 
specific R.H., placing it in so that the organism is not disturbed (which 
should not be too difficult). This experiment would then show whether the 
organism accommodates to a new R.H., and, if it does, by how much, in 
what time, and how much the accommodation depends on the initial and 
final R.H. The experiment may be modified by allowing the organism to 
walk from the initial R.H. to which it had been adapted into the test 
RH. In this method the organisms are moving initially instead of being 
initially inactive as in the former case. 

7. Experiment of Pielou and Gunn (1940). In these authors’ work on 
the same organism, results are graphed in their Figure 4 which shows the 
time course leading to the equilibrium situation. If we divide all the values 
in the graph by the maximum percentage found, 92, we then have a 
graph analogous to that discussed previously in Bursell and Ewer; hence 
(54) should apply. However, upon closer examination of the curve it can 
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be seen that its agreement with (54) is only very approximate, the best 
agreement being for a ¢* of approximately 5 minutes. We shall attempt to 
evaluate * by the use of (47), will take y = .6, 7 = 2, L = 22 cm, 
(c7)o = 2 cm., & = 50 cm./minute, and assume a narrow distribution of 
the c’s and 7’s so that (@r)o = (é7)oép. The value of (cr)) was found by 
examining the paths of the organism in Gunn and Pielou’s Figure 4 
(loc. cit.), while the value of é was found by personal observation of 
several Wood Beetles of a different species. This latter value may be in 
error due to species differences, but since we are only interested in orders 
of magnitude this should not be too serious. Inserting the above values 
into (47) we find that ¢* ~ .6 minutes, which is in poor agreement with 
the ¢* found above. 

A major reason for the lack of agreement between theory and experi- 
ment is probably that of the organism’s accommodation. However, the 
difference between theory and experiment may also be due to such facts 
as the large difference in the values of the organism’s parameters between 
the two halves of the chamber, whereas the theory requires small differ- 
ences; the inactivity of the organism part of the time, a factor which was 
not included in the theory; and the approximation method by which the 
equations were derived. The fact that Bursell and Ewer’s data (discussed 
above) are in such good agreement with theory may therefore only be an 
accident in which the different factors tend to cancel each other. 

8. Discussion of avoidance reaction and experiment of Bentley (1944). 
In the three papers which have been discussed, namely, Kennedy, Bur- 
sell and Ewer, and Gunn and Pielou, there was strong evidence of an 
avoidance reaction contrary to the results of Waloff’s paper. The primary 
reason for this is probably that in the former papers there was a very 
steep gradient of stimulus, whereas in Waloff’s experiment there was a 
very slight gradient of stimulus. This assumption of a very steep gradient 
inducing an avoidance reaction is in agreement with the observations of 
various authors. For example, Wigglesworth (1941) states: 

It is obvious in some experiments, however, that the response at the boundary, 


where there is a steep gradient in the intensity of stimulation, is a directed 
response. 


Lees (1948) gives an explanation of this as follows: 


The reason is clearly that in the linear apparatus the gradient of about 1° C/cm. 
is so slight that the tick cannot perceive without further exploration whether it 
is rising or falling. When the gradient is steeper, this information is provided 


immediately by the testing movements of the forelegs, and the tick orients 
accordingly. 
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Another possible reason for the inducement of an avoidance reaction by 
a steep gradient, and for the relation between avoidance reaction and 
stimulus intensity, was also put forward by Lees (1943) when he stated: 


At low humidities the larvae are always intensely active, and as a result, those 
which are crossing the arena are rapidly carried over the steepest part of the 
humidity gradient into the lower, though relatively uniform, humidity in the 
opposite side of the chamber. This has the effect of virtually eliminating the 
avoiding mechanism which is so well shown at higher humidities. 


A different explanation has been introduced by Fraenkel and Gunn (1940) 
on page 18 of their book, where they say: 


Shock reactions may be observable when the change of intensity is rapid 
enough, and there may be no reaction at all if sensory adaptation keeps pace 
with the change. 


There are, however, many cases of an avoidance reaction to a slight 
gradient of stimulus. One of these can be found in the experiment of Lees 
(1948). Another is found in the experiment of Bentley (1944) on Ptinus 
Tectus. In this experiment the animals were first placed in a ring con- 
tainer, where it was found that by the second day the organisms were con- 
centrated at the point of lowest humidity. Due to the length of time of the 
experiment and from the author’s discussion, it can be assumed that the 
orientation involved is due to differences of activity at different humidities 
and to an avoidance reaction. 

It should be noted that although Bentley finds that o,(X) changes over 
a period of 18 days for the organisms in a constant R.H. chamber, this 
does not mean that the organisms accommodate in this experiment, since 
we may consider the total experiment to bea series of shorter experiments, 
each of which lasts for about one day. During each day o.(X) remains 
fairly constant. In this case accommodation means adapting to a new 
stimulus, the adaptation occurring within period of time small compared 
to one day. The specific length of time depends on the gradient size and 
the speed of the organism, which will be discussed further in the next sec- 
tion. Bentley gives no evidence of whether or not the organism adapts al- 
most immediately or in a short time on being presented with a new stimu- 
lus. However, we shall neglect the effect of accommodation. 

If the assumption of no accommodation is correct, then, using the data 
shown in Figure 2 of Bentley’s paper together with (50), we see that the 
resultant values do not yield a congregation of organisms as intense as 
that found or an LR. as large as that shown in his Figure 4. It may thus 
be inferred that there is an avoidance reaction even though there is a 
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slight gradient. On the other hand, the congregation of organisms may bé 
due to accommodation or to the animals’ preference to remain near each 
other. The latter factor would obviously accentuate congregation and 
could possibly account for the results found. This could easily be tested 
by having only one organism in the ring gradient at a time, as was dis- 
cussed in the section on the evaluation of parameters. 

In the experiment of Gunn and Walshe (1942) on Ptinus Tectus Boie, 
the authors report that 
With Ptinus in a temperature gradient, however, it is the steep gradient which 
fails to produce a reaction against cold and the gentle gradient which succeeds. 


This result seems opposed to all expectations, but the existing data are 
too meager for any precise statements. It is hoped that further experi- 
mental work will be done on this organism so that this confusing situation 
may be clarified. 

It should also be pointed out that, contrary to the results of Waloff 
(1941), the orientation of lower organisms almost always includes the 
avoidance reaction. In fact, as Mast (1938) states: 


There are two fundamentally different types of photic response in animals and 
both are found in nearly every animal that responds to light. The one is quanti- 
tative. It consists primarily of increase or decrease in the rate of activity, corre- 
lated with the amount of light absorbed by the photosensitive substance. 
It may be designated a kinetic response. The other is largely qualitative. It 
consists primarily of change in the direction of movement, correlated with the 
rate of change in the amount of light absorbed by the photosensitive tissue. 
It may be designated a shock-reaction. 


Discussion. Some amplification of the assumptions involved in the 
derivation of the equations in I and in this paper may prove of value in 
view of the applications of those equations made here. The six assumptions 
listed in I will be considered first. Assumption 1, the negligible interac- 
tion between organisms, has been thoroughly discussed in the section on 
evaluation of parameters. Assumption 2, the non-persistence of ¢ and 7, 
may not be completely valid for all actual situations. However, this error 
is not too important for the majority of experiments, since the standard 
deviation of the usual distribution of speeds is very much less than the 
average speed of the organism. The same argument applies to the dis- 
tribution of r’s. Assumption 3, the negligible time spent in turning, should 
cause no difficulty. If this assumption is found to be false in an experi- 
ment, then the time the organism spends in turning may be considered as 
inactive time and the problem may be treated by the methods outlined in 
the section under inactive organisms. Assumption 4, the slow change of 
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the system, should present no problem because we are usually only inter- 
ested in the steady state, which is by definition a state wherein the system 
does not change at all. For the non-steady state, the lack of validity of this 
assumption might prove troublesome. However, since the equations are 
solved by means of the approximation method, the error involved might 
be less than the error involved in the derivation of the approximation 
method equation. Assumptions 5 and 6 concern the slow change and 
continuity, respectively, of the functions in I. If these assumptions are 
not valid in a given case, the problem may be treated by means of the 
avoidance reaction method. 

An implicit assumption in I was that the organism traveled in a straight 
line between turns. This is found to be valid for most experiments. How- 
ever, if there is doubt as to whether or not an organism made a slight turn 
the experimenter will have to decide the issue. Then the values of 7o(X, é) 
and y(X, #) will be modified so that if one increases the other decreases 
and vice versa. Yet it may easily be seen from (5) that these changes will 
tend to cancel each other out, so that the value for V(X) willremain about 
the same regardless of the decision as to whether or not the organism 
turned. 

Thus it can be seen that all of the assumptions in I discussed above are 
either valid for experimental situations or, if not, the situation may be 
treated by use of the methods introduced in this paper. Those assumptions 
mentioned in the derivation of the coefficient of persistence of direction 
will not be discussed here due to the lack of experimental data on this sub- 
ject. However, as in our analysis of the experiment of Ullyott (1936), 
any data would have to be analyzed to see if these assumptions are valid 
before the equations developed in this paper may be applied. 

The assumption of (1) in this paper, concerning the fact that the mean 
free path of the organism is much less than the length of the container, is 
not generally valid for most of the experiments which have been carried 
out. Some examples discussed in the previous section are those of Waloff 
(loc. cit.), in which (cr) = 8 inches, while L = 29 inches; Bursell and Ewer 
(loc. cit.), in which (¢r)) = 8 cm., while L = 22 cm.; and Gunn and Pielou 
(loc. cit.), in which (ér)9 = 2 cm., while L = 22 cm. If by “much less than” 
is meant less than 1/20, then it can be seen that (1) is satisfied in none of 
these examples, and thus the results derived for all of the above examples 
cannot be considered completely valid. However, to be sure that the equa- 
tions of this paper apply, the experimental conditions must be such that 
the assumptions are fulfilled, even if this is found to be difficult. 

As stated before, if the organism exhibits accommodation, the equa- - 
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tions of this paper are not applicable. However, for the situation in which 
the organism accommodates either very slowly or very quickly, informa- 
tion may still be gained for the steady state. For the case of slow accom- 
modation to each new stimulus, if the organism enters and leaves a region 
within a time much less than the time of accommodation, we may con- 
sider that the organism does not have any accommodation at all. If the 
organism accommodates slowly to the entire situation, the experiment 
may then be divided into many sub-experiments, each lasting sufficiently 
long for enough data to be gathered for valid conclusions, but short 
enough so that the organism does not accommodate by much. The value 
of the parameters would not change by much, and therefore the average 
value of the parameters for each sub-experiment may be used, and their 
change with time subsequently noted. 

For the case of rapid accommodation, the equations of this paper may 
still be used. If the organism accommodates quickly to each new stimulus, 
then each time it crosses the boundary between the two regions of the 
alternative chamber it will adapt to the new region. If it adapts fast 
enough so that by the time it makes a few turns in a region it has adapted 
completely to the region, then the experiment may be treated by means of 
the avoidance reaction. The region a few mean free paths about the divid- 
ing line in which the organisms will adapt completely may be considered 
the region in which the organism possesses an avoidance reaction. The 
equations of this paper may then be solved with the use of the experi- 
mental avoidance function, but the value of the true experimental func- 
tion cannot be found for this case, since, as a result of adaptation, the 
ratio of the number of organisms in the two regions cannot be calculated 
for the case where there is no avoidance reaction. If the organism accom- 
modates quickly to the entire situation, then it is only necessary to wait 
a short time until it adapts completely and then perform the experiment. 

If the time of accommodation is between the two extremes mentioned 
above, the equations of this paper cannot be applied. It should be noted 
that in any experiment to test the organism’s accommodation to a new 
stimulus the total history of the organism is important. For example, an 
illustration of the effect of culture temperature on the organism’s adapta- 
tion to a new environmental temperature is found in Gunn and Hopf 
(1942), while another example, the change in the rate of accommodation 
as an organism is repeatedly presented to the test stimulus, is to be found 
in Lees (1948). 

In conclusion, it may be stated that we have attempted to derive equa- 
tions in this paper so that if a prescribed experiment is carried out and 
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various measurements made, computation of certain parameters and their 
insertion into the equations will give values for the number of organisms 
at specific points in the container of the experiment. If these values then 
agree with the values found by direct measurement, we may infer that 
our equations are valid for this situation. If the values do mot agree with 
experimentally found values, then it may be inferred that our equations 
are not valid for this case, and therefore further experiments must be 
planned to attempt to find the cause of the discrepancy. 

It is hoped that this paper will provide a systematic method for experi- 
mental investigation of the problem of the orientation of organisms. It is 
also hoped that by use of the knowledge of the various parameters not 
only may the behavior of the organism in its natural environment be 
predicted but also the search for the underlying physiological mechanisms 
for this behavior may be facilitated—for, as Mast (Joc. cit.) has said: 


These and many similar facts make it evident without argument, that classifica- 
tion of organisms in reference to the process of orientation must necessarily be 
so indefinite that it has but little value. It would therefore seem to be much 
wiser to attempt to ascertain precisely what occurs in the organisms in various 
physiological states, during the process of orientation under various accurately 
controlled conditions, than to attempt to ascertain whether the process of 
orientation in this or that organism should be designated thus or so, e.g. photo- 
taxis, tropotaxis, telotaxis, or menotaxis in the terminology of Kuhn. 


The author wishes to express his appreciation to the Committee on 
Mathematical Biology for their generous aid and assistance in the prepa- 
ration of this paper, and, in particular, to Mr. Robert Macey and Dr. 
H. D. Landahl for their many stimulating discussions and encouraging 
remarks. 
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The equations governing the time course of the exchange of substances between the blood 
in the capillaries and the extracellular space are solved for the case of substances which do 
not penetrate the cells. The equations given relate the time course of the exchange process to 
the various tissue and circulation parameters such as the specific capillary wall area, the pore 
area, the inter-capillary distance, the size of the extra-vascular, extra-cellular space, the dif- 
fusion coefficient in this space, and the velocity of blood in the capillaries. 

Some experimental work on capillary exchange is discussed in relation to the theory and 
estimates are made of the relative importance of the various tissue and circulation parameters 
in the exchange of substances in different tissues. 


Introduction. The problem of the exchange of substances between the 
blood in the capillaries and the intercellular space, and between this 
space and the cells has been considered previously by many investigators. 
A brief review of some of the work has been given by the present writer. 
(Schmidt, 1952, hereinafter referred to as I.) 

Although the problem is a complicated one and quantitative descrip- 
tions of the capillary anatomy and of the velocity of blood flow in various 
tissues are lacking for the most part, some data are available. Such data 
may be used as guides in constructing an idealized tissue model, and the 
equations governing the exchange relations in such a model may be formu- 
lated. Such a model has been constructed and the equations deduced (1). 
The aim of the present paper is to review briefly the capillary anatomy; 
the tissue model, and the pertinent equations, to extend the previous 
treatment slightly and to present the solutions of the equations governing 
the time course of the exchange of substances which do not enter the cells, 
or which enter so slowly that the cellular concentration of the substance 
may be taken to be constant. The equations developed give the concen- 
tration in the venous blood, the mean capillary concentration, the mean 
extracellular concentration, and the transient osmotic pressure, all in 
terms of the concentration in the arterial blood and other measurable 


parameters. 
The capillary anatomy. The capillaries form an maastomnoeiig net-work, 
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with central channels and side-branches in some tissues. Arterio-venous 
shunts, which directly connect the smallest arterioles with the smallest 
venules, are present in most tissues. Two types of capillaries are distin- 
guished, the arteriolo-venular capillaries, and the true capillaries. The 
arteriolo-venular or a-v capillaries normally have a relatively rapid, con- 
tinuous flow of blood through them, whereas the true capillaries are 
marked by a slower and intermittent flow, with flow ceasing completely in 
various different true capillaries from time to time. The ratio of the num- 
ber of true capillaries to the number of a-v capillaries varies from one 
tissue to another and constitutes a tissue characteristic. The relative 
number of true capillaries through which blood is flowing at any time, in 
any given tissue, is dependent upon the metabolic activity of the tissue and 
is subject to large variations. Furthermore, the various tissues of an intact 
animal vary greatly in their specific blood perfusion rates (Jones, 1951). 

For example, the thyroid in man is perfused with 3.0 to 10.0 volumes of 
blood per unit volume of tissue per minute. This is about 150 to 500 times 
the perfusion rate of resting muscle; about 0.3 to 2 times the rate for the 
kidney; about 2.5 to 9 times the rate for the liver; about 6 to 20 times the 
rate for the brain; and about 20 to 60 times the rate for the marrow. 

In view of these data attention will be confined to a single homogeneous 
tissue, and the assumption made that the arterial and venous blood, enter- 
ing and leaving this tissue, is sampled. Furthermore, the distinction be- 
tween a-v and the true capillaries will be neglected. 

But, even in a homogeneous tissue, transport delays occur in the vascu- 
lar system. Thus we may sample the arterial blood at time #, but the blood 
entering any specific capillary, say the 7th, has a concentration equal to 
the concentration in the artery at the site of sampling at time ¢ — &i, 
where /; is the arterial transport delay, i.e., the time required for the blood 
to flow from the arterial sampling site to the 7th capillary. Similarly, the 
concentration of a sample taken from the venous blood leaving the tissue 
at time ¢ is a weighted mean of the concentrations of blood leaving the 
various capillaries at times from ¢t — 4 to ¢ — fn, where times hh, fo, . . . , tn 
are the venous transport delays of the blood in flowing from the various 
capillaries to the venous sampling site. 

Any theory which purports to describe capillary exchange ought then 
to take into account the two kinds of capillaries, in tissues in which the 
distinction can be made, the shifts of flow from one to another of the vari- 
ous true capillaries and the fraction of the blood which flows through the 
arterio-venous shunt, and hence does not go through any capillary. Further- 
more, the arterial and venous transport delays ought to be considered. 
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A theory taking into account the two types of capillaries has been con- 
structed (I). Only minor corrections need be made in this theory to take 
the arterio-venous shunts into account. However, the theory is of limited 
utility, since in order to use it, data are needed which are not at hand. 

The total transport delays in the hind limbs of cats may be roughly 
estimated from data supplied by E. M. Renkin (1952). No separation can 
be made of the arterial, capillary, and venous delay times, and only crude 
estimates may be made of the minimum, maximum, and mean total de- 
lays. These are, respectively, a few seconds, six to eight minutes, and one 
and one-half to four minutes, for a perfusion rate of about 6.5 X 10-4 ml. 
per sec. per gm. of limb. 

If the capillary exchange rate is sufficiently slow relative to the mean 
transport time through the tissue, we may simplify the problem by 
neglecting the exact distribution of transport times, and by assuming that 
the transport time through all the capillaries is equal to the mean time. 
This, of course, means that our treatment of the problem will be re- 
stricted either to extremely high perfusion rates, or, more usefully, to 
normal perfusion rates, but to large molecules with slow exchange rates. 

The tissue model and its exchange equations. The actual capillary anat- 
omy will be represented in the model tissue by an arrangement of capil- 
laries, all in parallel, each of length Z, and of radius r. Each capillary in the 
model runs through a cylinder of tissue of radius Ry and of length L. Five 
different concentrations are considered; namely, the concentrations of the 
arterial and venous blood, the mean capillary concentration, the mean 
concentration at the wall of the capillary in the interstitial space, and the 
mean concentration in the interstitial space. The cellular concentration is 
assumed to be zero at all times, i.e., the substance does not enter the cells. 
The exchange of the substance, across the capillary wall, is assumed to be 
governed by diffusion, as has been experimentally confirmed. (References 
in I.) 

The equation governing the capillary plasma concentration K(z, ¢), 
anywhere along its length, which runs from z = 0 to z = L, is 


LGD ee Se ee gs 
Ss = IK (e ) -E( 1. A) 


Here v is the velocity of the blood along the capillary, & is the permeability 
coefficient of the capillary wall, r is the capillary radius, f is the fraction of 
plasma in the blood, and £’(z, ¢) is the concentration near the capillary 
wall in the interstitial space. The equation is similar to equation (11) in I, 
except that the factor f was inadvertently omitted from the latter, v and f 
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are here immediately assumed constant along the capillary, and E' (gz, t) is 
here stated to be a function of z. 

At the wall of the capillary we have 

Rel Beit) — 2 ays =e [E’, (t) —E(t)). (2) 

Here D is the free diffusion coefficient of the substance in the interstitial 
fluid, gis the ratio of the cell-free area through which diffusion may occur 
to the total area normal to the gradient of the concentration in the inter- 
stitial space, / is the ratio of the short distance between two points in the 
tissue to the distance which must be traveled between two points without 
crossing any cell walls, K»,(¢) is the mean plasma concentration, £,,(é) is 
the mean concentration along the wall of the capillary in the interstitial 
space, and E(¢) is the mean concentration in the interstitial space. The 
equation is approximate in that the exact gradient at the capillary wall is 
replaced by an approximation, and the speed of diffusion through the 
interstitial space is assumed to be very great relative to the speed of trans- 
port through the cell walls. 

In the interstitial space we have 


2 
ee [K, (t) -E‘()]. (3) 
In this equation a is the fraction of the extracellular, extravascular vol- 
ume, not occupied by cells and available for the distribution of the sub- 
stance. 

The final equation is one giving approximately the osmotic pressure 
acting across the capillary wall due to the substance under consideration. 
This is 

AP=RT[K, (i) —E’ (i)]. (4) 


Here AP is the osmotic pressure, R is the gas constant, and T is the abso- 
lute temperature. 

The solution of the equations. In order to solve the set of equations, (1) to 
(4), we use the technique of the Laplace transformation. We multiply each 
equation by e~?' where the real part of # is positive, and integrate with 
respect to ¢ from 0 to ». Let the Laplace transforms of the various func- 
tions be given as follows: 


K (2, t) ~-w(2, p); 
Ky (t) wm (Pp); 
E' (2, t) +x! (2, p)5- 
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El (t) +x! (p); 
E(t) +x (p); 
APU Cas 


Let us also introduce the following notation 


B=p+uh, 

6 =1— e~(4/) (e+) = 1 — e—8, 

3% 
ee 

Meet Sa it 

_2D¢l (5) 
Pav iby 

2rk 

oo Rr)? 
T=-. 


v 


Then we obtain the transforms of (1) to (4) as follows 


d 
pro(s, p) = — 9 WEEP) san(s, p) + ux"(2, 6), (6) 
w, (P) — #1, (p) = plz! (p) — 2 (9)], (7) 
px(p) =¢ lw, (p) —2!,(P)], (8) 
and 
u (p) =RT (w,, (p) — 2%, (p) 1. (9) 


Here the initial concentrations are assumed to be zero. 
We are not interested in £’(z, #) itself, but only in its influence on the 
other functions. Thus the solution of (6) is 


w (2, p) = woe Jw (0, p) +2 fox’ (2, p) ewmedst, (10) 
but we may rewrite this as 
o(z,)) = e198, (p) +2", (o) 5 (1 — ¢-(e/s)8) t (11) 


where x,,(p) is the mean value of x’(z, p) from ¢ = 0 to z = L and wo(p) 


= w(0, ). 
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From (11) we obtain the equations for wn(p) and w(L, p) = w,(f). 
These are 


and 


Introducing (13) into (12) we have 


w,+rbw,— Tux, =w, - (14) 


Equations (6), (7), (13), and (14) become 


w,, vx + px =O, (15) 

— ¢w,, tox), +px=0, (16) 
bw, — pox =8(1—f)»w,, (17) 
w,+Thw,— Tux, =, - (18) 


Let us designate the determinant of the coefficient of the unknown 
functions in this set of equations by A; then we have 


0 il —v p 

Sic $ p 
A= 

8 ts 0 

1 7B —Th 0 


or 
A= -wjp(o+u)[pt2o+u [+ oteot. (19) 


Routine calculation gives the following 


wop (P+u)[p+2(6+u) | 
M0 Wi a ae pee pea eT eee (20) 
p(p+u)[pt+2ot+u) |+#£(p+¢0) 
u=RT wopp (btu) £ 
p ue P (21) 
vw |p(o+u) [p+2o+n) +2490) } 


and 
‘Sent wodp (pt nu) & 


—_ (22) 
wv jp(o+u)[pt+2o+u) +X (p+) 
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pale (p+¢p) (23) 


™plo+u)[pt2o+u) |] 


then we may rewrite (20), (21), and (22) as 


Wot 
aad a pa ae (24) 
aye ered res (25) 
| p+2(6+u) | +90 
and 


vp p+ 2 (6-+u) | +n) 


If the real part of # is sufficiently large, |n¢| < 1 and we may write 


1 = n,n 
uy sag os amie hc (27) 


ape ews 97 
ae: ee epee eg ieee yas ae ay 
Letting 
y 
Fay eae age chp ge car ae ce Se cea poy LD 
p(p+u)[p+2o+u) |+£eteo) 142 °7" 
Ty n 
we have 
n—1 
ep eth. 40) 


The inversion of (30) is readily obtained in terms of partial fractions. 
Using the results in (28) and this in (24), (25), and (26), we may evaluate 
exactly the functions of interest. However, the results are so complex as 
to be of little use. Besides the complexity of the expressions another con- 
sideration lessens the interest in them. This is the complication introduced 
by the time delays between the arterial sampling site and the capillaries, 
and between the capillaries and the venous sampling site. The equations as 
developed describe the process at one capillary. In order to describe the 
whole process of capillary exchange even in a homogeneous tissue, where 
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all the capillaries are similar in length, spacing, permeability, etc., one 
would need to know the distribution of transport delay times, from the 
arterial sampling site to the capillaries, in order to introduce the correct 
concentrations at the entrance to the capillaries. Then one would need the 
venous transport delays in order to know how to combine properly the 
expressions for the concentrations in the blood leaving the various capil- 
laries. 

The needed data are not available. But, even if they were, they would 
not help much, for the initial details of the exchange involve some fast 
processes, with time constants of the same order of magnitude as the trans- 
port delays. For the rate constants of the exchange processes are given by 
the roots of the equation A = 0, where A is given by (19), and the smaller 
roots are less than —(p/v)(@ + yu) which is likely to be of the order of mag- 
nitude of 10-? to 10-* per second. But the slow processes are so slow, at 
east for relatively large, lipid soluble molecules, that the transport delays 
are negligible. This point will be discussed later at greater length. 

We will confine our attention then to equations (19), (20), (21), and 
(22). Combining (20) and (21) and inverting the results gives 


dAP 
aRE dt Ont t (btu) on= Ko— Kx, (31) 
or 
PRES RT ( Ko— Kx 
: + (+n) =A (A). (32) 


In order to compare this with a somewhat similar expression given by 
J. R. Pappenheimer e¢ al. (1951) we need an interpretation of the con- 
stants p, v, ¢, uw, and 7 in more familiar terms. Let us assume that a unit 
volume of tissue has WV capillaries, each of length L. Then if each capillary 
supplies a cylinder of tissue of radius Ro we have R2LN = 1. If the radius 
of each capillary is 7, then the capillary wall area per unit volume of tissue 
is A = 2nrLN. If the average linear velocity of the blood in each capillary 
is v, then the total perfusion rate, per unit volume of tissue is V = mr2vNV. 
If fis the fraction of plasma in the blood, then the plasma perfusion rate is 
Q = fV. We have neglected the volume occupied by the blood vessels, and 
any arterio-venous shunts. Let us also assume that the permeability of the 
capillary walls is due to pores (cf. I and Pappenheimer, loc. cit.). Then we 
may write 

_D Ao 

uohte Aa’ 
where D is the diffusion coefficient of the substance exchanging and 
Ao/Ax is the virtual pore area per unit path length through the pores in 


(33) 
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the capillary walls. Thus Ao/A is the ratio of the virtual pore area to the 
total capillary wall area and Az is the effective length of the pores. 
We can now write the following equations 


TP ara a f" 


ua 20. LReae 


_2Dgl_ 4gir 
peek Aw peck okt (34) 
Ax 9 


Using the expressions for the various constants in (29) we obtain, approxi- 
mately, 


Ao RB 
e(, 4m ) 2naP , de ORT (Ke) (35) 
DR Argl di eo D AP ‘ 
Neglecting for the moment the first term on the left, we see that this is the 
same expression used by Pappenheimer ef al. (loc. cit.) to evaluate 
A/Ax, except that the values of the concentration at the entrance and 
exit to the capillaries are involved, rather than the arterial and venous 
concentrations. After a sufficient time this difference is very small. An 
inspection of all the osmotic transient data presented by Pappenheimer 
also indicates that the osmotic transients are well represented by a single 
decreasing exponential after about five minutes from the start of the ex- 
periments. Using this data and reasonable values for the constants in the 
left-hand term of (32), we find that this term may be neglected. Hence the 
theory here presented may be regarded as having some experimental con- 
firmation since the equation (35) involves all of the assumptions of the 
whole theory. 

Returning now to equations (20), (21), and (22), and assuming for the 
moment that Ky is a constant, we find that for large ¢ all three functions, 
Ko — K,, AP, and E vary as does 


( 
+ aCe 


ppl (36) 
et 


es) 
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For ru X 1, ie., for a very great blood velocity or small permeability of the 
capillary wall, and also neglecting ¢ with respect to » we have 


expj =P | (37) 


If tu>> 1, ie., for small blood velocities, or large permeability of the 
capillary, we have 


exp] 1+ 7 py . 

The expression (37) was given in I, as were graphs of the time constant 
as a function of Dgl, and , for various values of Ro. The expression (38) is 
probably more widely applicable. In terms of the constants defined by 
(34) we can write (37) and (38) as follows: 


exp {| — ——+—_— ru<1 (39) 
aAx aRo te 
DA. 4rDgl : 
and 
= ; ru>>1 (40) 
a aR}, 
Q' 4rDgi Se 
If 
S igh 
mast 
we have for (36) 
exp (41) 
a aR} F 
— (1 — e—(B/Q)(Ag/ Az)— 
ie a Saepel4 


Pappenheimer et al. (loc. cit.) have determined the value of the factor 
A,/Ax in the perfused hind legs of cats for various substances. The values 
given in their paper in units of cm. per gm. of muscle at 36° C. follows: 
NaCl, 1100; glucose, 750; sucrose, 580; raffinose, 420; and inulin, 180. 
The free diffusion coefficients of these substances in water at 37° C. are in 
units of cm.? sec. X 105: NaCl, 2.03; glucose, 0.90; raffinose, 0.64; sucrose, 
0.75; and inulin, 0.21. The plasma flows used were such that the values of 
DA)/QAx are greater than 4 for all substances except inulin and one 
experiment with raffinose. 

In general the expression (39) may be applicable to the exchange of 
large lipid insoluble molecules, M.W. > 1000, in tissue with large perfu- 
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sion rates, Q > 10-* ml./sec./gm. of tissue, while (40) is applicable to the 
exchange of small lipid-insoluble molecules, M.W. < 500, in tissues with 
moderate or small perfusion rates, Q < 10-3 ml./sec./gm. of tissue. Ex- 
pression (41) may be applicable to most tissues, for all lipid-insoluble mole- 
cules, while the general expression (32) ought to be applicable to all tissues 
for all lipid-insoluble molecules which do not penetrate the cells. 

Conclusion. The theory developed here may be useful in estimating the 
relative importance of blood flow rate, capillary permeability, capillary 
spacing, and diffusion coefficients in tissue in the exchange process. It is 
evident that for a sufficiently small lipid-insoluble molecule the exchange 
rate always depends on the blood flow rate [cf. (40)]. The influence of the 
flow rate relative to the capillary spacing and the diffusion coefficient in 
the tissue will of course vary. 

The factor aRo/4rDgi may be written as 


aR) _ anRoL ce (42) 
4rDgl 2nrk \2Dgi 
Then a7R;L is approximately the tissue volume per capillary available for 
the distribution of the substance, 27rL is the capillary wall area of each 
capillary, and 2D/Rp is a “tissue permeability.” 

The factor aAx/DA, is approximately 


iia (43) 


Here & is the capillary permeability. 

Thus in (39) the time constant is the product of the ratio of the avail- 
able tissue volume to the capillary wall area by the sum of the reciprocals 
of the capillary and tissue permeabilities. 

The factor a/Q is approximately 


a7RoL ‘ (44) 


Thus a/Q is the ratio of the available tissue volume per capillary to the 
volume of plasma that flows through each capillary per second. The time 
constant of (40) is then the product of the available tissue volume per cap- 
illary by the sum of the reciprocal of the volume of the plasma flowing 
through the capillary per second, and the reciprocal of the product of the 
capillary wall area by the tissue permeability. 
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The factor ru which appears in (32) and is equal to the factor DA»/QAx 
in (41) can also be written as follows: 


sc DAo _2xrLk 
TE OAn ar? fv’ 


(45) 


Thus rp is the ratio of the product of the wall area per capillary by the 
permeability coefficient of the wall to the plasma flow through the capil- 
lary per second. 

The analysis given here may be extended to substances which penetrate 
the cells and which are consumed or produced in the cells. The results are 
complex and of limited utility without a better knowledge of the parame- 
ters used in this paper. However, if experiments are done which yield the 
parameters used here, then a further development of the theory, dealing 
with substances which do penetrate the cells may be useful. 

The author is indebted to Dr. H. G. Landau for a critical reading of this 


paper. 
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A mathematical theory is developed which permits the determination of certain parameters 
of an inhomogeneous tissue, such as a nerve trunk without its epineurium. The parameters 
are the permeability coefficients for entrance into an exit of a substance from the nerve fibers, 
and the diffusion coefficient of the interstitial material. The experimental data required are 
the dimensions of the cross-section, the average diameter of the fibers, and the ratio of the 
cross-sectional area of the fibers to the total cross-section, as well as the time course of the 
decrease of the fraction of the substance left in the nerve trunk, when the trunk is immersed 
in a bathing solution containing none of it. 


Introduction. The determination of coefficients which quantitatively de- 
scribe the exchange of substances between a biological system and its 
environment seems to be dependent upon the development of theories 
which take into account the complex geometry of the system, or the isola- 
tion of the cells of the system from the interstitial matrix without damage. 
The latter procedure is at present impossible with most tissues. Hence 
theories of the exchange process which are realistic in considering the 
geometrical details of the tissue dealt with would be useful. 

The present paper is concerned with the exchange processes between a 
system of circular cross-section, very long relative to its radius. The sys- 
tem is composed of longitudinal fibers, separated from each other by an 
interstitial matrix. A nerve trunk with its epineurium removed may be 
approximately so described. It is hoped that the exchange process in a 
nerve trunk with its epineurium intact will form the subject of a later 
paper. 

Description of the tissue model. We consider then a system, which for 
simplicity will be called a nerve trunk, which is of circular cross-section 
with radius R, and which is so long that we may assume that any exchange 
between the trunk and its environment takes place in a direction normal to 
the axis of the trunk. The nerve trunk is made up of round fibers, with 
their axes parallel to the trunk axis, separated from each other by an inter- 


stitial matrix. 
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The interstitial matrix is characterized by a diffusion coefficient D, 
which must be modified by two factors. The first factor must take into 
account the reduction in the area, normal to the flow path, across which 
material may flow without entering the fibers. We define the factor g as 
the ratio, to the total area, of the area normal to the concentration 
gradient in the interstitial material through which diffusion may occur 
without entering the fibers. The second factor / takes into account the 
greater distance which a molecule must move in traveling from one point 
to another without passing through the fibers than if it passed through the 
fibers. We define / as the ratio of the shortest distance between two points 
in the trunk to the average distance between the two points which must 
be traveled while not passing through any fiber. Both g and / must be com- 
puted over distances equal to several fiber diameters, and in a direction 
normal to the trunk axis. With these definitions the effective diffusion 
coefficient in the interstitial material of the trunk is Dgl. We also consider 
that the interstitial material may not all be available for the distribution 
of the substance and define 6 as the fraction of the interstitial material 
which is available. 

Distributed throughout the trunk are fibers of average radius p. This 
assumption is the weakest point in the theory, for the fiber radii in a nerve 
trunk are certainly not all equal, and the distribution of radii is probably 
bimodal or trimodal. However, the introduction into the theory of more 
than one fiber size, while possible, would yield such complicated mathe- 
matical results as to be of doubtful utility at the present time. 

Each fiber is characterized by two “permeability” coefficients, &; and 
ko, such that under equilibrium or steady-state conditions the concentra- 
tion inside the fiber, C, is related to the concentration in the interstitial 
space adjacent to the fiber, HE, by the equation k\E = k,C. These co- 
efficients may be passive permeability coefficients modified by a partition 
coefficient, or may be related to the rate of an active transport mechanism. 

It is assumed that of the intra-fiber volume, a fraction, a, is available for 
the distribution of the substance. The total cross-sectional area of the 
trunk is divided into two fractions, the fraction composed of fiber area, fa, 
and the fraction composed of interstitial area, fs; we have f, + fr = 1. 

The experimental conditions assumed. The experimental conditions as- 
sumed are as follows: The nerve trunk, with its epineurium removed, con- 
tains the substance whose transport is to be studied in equilibrium con- 
centration, at time zero. This condition is attained by immersing the trunk 
in a solution of the substance until each of the two compartments has 
attained equilibrium everywhere in the trunk. Then if E is the uniform 
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interstitial concentration, and C is the uniform intra-fiber concentration, 
we have ki = k.C. The experiment with the epineural sheath removed is 
convenient in avoiding the very great mathematical complications arising 
from the presence of the sheath, but is itself complicated by the swelling 
which occurs when the sheath is removed. However, if the sheath is re- 
moved at the beginning of the immersion process, the preparation becomes 
approximately stable dimensionally before the soaking-out process is be- 
gun (Trunat, 1953). But it must be remembered that the swelling changes 
the characteristics of the interstitial material. 

At time zero the nerve trunk is immersed in a soaking-out solution, 
whose concentration in the substance considered is zero, and remains ap- 
proximately zero. This may be accomplished by using a well-stirred solu- 
tion of large volume, or by using smaller volumes which are changed fre- 
quently. In either case the experimental data should give the fraction of 
the initial amount remaining at any time after the beginning of the soaking- 
out of the substance. 

The exchange equations. The fundamental idea of the method to be used 
is similar to that used by I. Opatowski and G. W. Schmidt (1952), suitably 
modified for the circular cross-section. 

We divide the trunk, of radius R, into N concentric circular rings each 
of width 7, so that rV = R. We calculate the exchange occurring through 
the interstitial space between each ring and its neighbors, between the 
interstitial space of that ring, and the fibers of each ring, and between the 
interstitial space of the outer ring and the solution in which the trunk is 
immersed. 

Let the average concentration of the material in the interstitial space of 
the nth ring be £,(#) and that in the fibers of the ring be C,,(¢). Then, ap- 
proximately, the amount entering or leaving a one centimeter length of the 
fibers of the th ring in time df is given by 


(hE, — BC.) “fort? (20 — 1) dt. (1) 


If we assume that the core has M fibers per square centimeter of cross 
section and pis the fiber radius, then zp is the total fiber area per square 
centimeter which is equal to f,. We also have 27pM as the total perimeter 
of the fibers per square centimeter of core. Hence 

m=, (2) 


ap?’ 


and 2f,/p is the fiber perimeter per square centimeter of core, 
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The amount of material accumulating in or leaving the fibers of the th 


ring in time dé is approximately 
dC 


ar?(2n—1) faa ry din Gen) 
Equating (1) to (3) and simplifying we have 


On...” Rp hia ee (4) 
Ot ap 


The amount of material flowing from the (m — 1)st ring to the mth ring 
in a one centimeter length in time df is approximately 


(FAS) Deltas (n— 1) at; a= 2, to Nl 


Here the first term is the approximate gradient between the two rings, 
Dgl is the diffusion coefficient, modified by the geometrical factors arising 
from the fiber spacing, and 2rr(m — 1) is the circumference of the bound- 
ary between the two rings. Similarly the amount flowing from the th to 
the (n + 1)st ring is 


(“*) Dglderndi ; (6) 


The amount of material entering or leaving the fibers of the mth ring is 
given by (1). The amount accumulating in the wth ring in time di is 


Pr eA Bus (7) 


Equating (7) to the sum of (5), the negative of (6), and (1) and simplifying 
we have 


dE, 2Dgl 
2 + Ses _—" 
(2n — 1) ry 


= Fag (AEA (20 = 1) Eg + Ena] 
(8) 


(ner) 2a eet eee ae Dees 
pbfy 


Special considerations are necessary for the 1st and Nth rings. 
The amount of material flowing across the boundary of the 1st ring to 
the 2nd ring in time d¢ is 


E, —E, | 
(Ao : )Del2mrdt. (9) 
The amount entering the fibers of the 1st ring is 


2 
(RE; — kC)) pietriat : (10) 
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The amount accumulating in the interstitial space of the first ring in time 
dt is 

OF 
Es 


Equating (11) to the negative of the sum of (9) and (10) we have, after 
simplifying, 


wr fyb dl. (Gla 


dE, + 2Dgi 


aS 
CU ELYA mehr ey perce o tg 


The amount of material entering the Nth ring from the (NV — 1)st ring 
in time df is 


Ey1—E 
(2 *) alter (WA) at. (13) 


The amount leaving the Vth ring into the soaking-out solution is 


ie. (14) 
2 
The amount entering the fibers of the Vth ring is approximately 


(hiEy — FxCx) “fon? (2 = 1) dt. Gis 


The amount accumulating in the interstitial spaces in the Nth ring is 


approximately 


Qn (2N—1) bfy 22 dt. 2-116) 


The material balance requirement gives, using (13), (14), (15), and (16) 
and simplifying, 


7 gees eee 
(2N e ot r-bf, 


aye 1) ole hyBy— byl 


[(N —1)Ey-1— (3N — 1) Ey] | 
(17) 


The initial conditions are, at ¢ = 0, 
E, =E, a constant form =1,2,..., N; 
C,=C, a constant form =1,2,...,N; (18) 
hE, = ke, « 


The solution of the equations. In order to solve the set of equations we will 
use the Laplace transform, i.e., we multiply each equation by e~?‘dt and 
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integrate from zero to infinity. Here p is a number whose real part is posi- 
tive. We symbolize the Laplace transform of any function F(t) by L[F(@)] 


(Doetsch, 1950). 


Let us then set L[C,(é)] = xn(p) and L[EZ,(é)] = y»(p). It will be con- 
venient to adopt the following notation for the combinations of constants: 


Then the transform of (4) is 


$5) Ca he ee 
ap 


Solving for x, and using (19) we have 


Di 
pane eek? oh jot Et 
; p+ 2 b+ d2 


The transform of (8) is, using (19) and (21) and simplifying, 


= 2 p 
(n= 1) ya — 2n—1) [1 +8 PU4+se. Jon tyne 


=— (2n— NE (1+5H); . 
2 


The transform of (12), using (19) and (21), is 


= {Eee eg 


The transform of (17), using (19) and (21), is 


(N= 1 ywa-[w+ @v—-1) f1+2 (1459 Vy 


xia 2st} FUt525 


Using the notations 


a=1+2(145%.), a= 7 (14585 ; 


Ee 


1 
+ G2)” 


(19) 


(20) 


(21) 


(22) 


(23) 


(24) 


(25) 


NERVE TRUNKS 495 


we may rewrite (22), (23), and (24) in a more convenient form as follows: 


(a — 1) Yn—1 — (2% — 1) arya t+ nyn41 
(26) 
= — (2n—1)a2., Renan s+ Nl) 
— aiyit yo= — az ; (27) 


(N—1) yn~-1— [N+ (2N—1) ail yw= —(2N—1) a2. (28) 


The complete solution of equation (26) is given by 
Yn= AsPys(a1) + AsQmi(an) +2, = 2,3,-..,N—1; (29) 


where P,,1(a;) and Q,-1(a1) are the Legendre polynomials of the first and 
second kinds, respectively; A; and A2 may be functions of p, but do not 
depend on n. 

In order for (29) to satisfy (27), A2 must be zero. Substituting (29), with 
A, = 0 into (28), we find that 


2E 1 
ees Ue HES cone oe 


Hence the solution of the set of equations (26), (27), and (28) is 


a a 2P,—1 (a1) +e 
m= (1 CHRIS ICH PAS ey lad Oe dete 6 


Putting (31) into (21) we find 


>. ki E de 2 (a1) wall 
»=25[1-GSe) me aac mes Rane” 


Calculation of the fractional amount remaining at any time. Let us now 
consider the amount of the diffusing substance, Q(t), remaining at any 
time in a one centimeter length of the nerve trunk. At time? = 0, we have 


Q(0) =7R°bf,E+ mR’af.C ; 
or 

Q(0) = aR /jB(1+2). (33) 
Let X,,() be the amount of the substance in the fibers of the mth ring at 


time ¢, and Y,,(#) be the amount in the interstitial spaces. Then 


N 
OW) = >> 1xX,0 FY WT: (34) 
n=1 
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Since the area of the mth ring is 7r?(2n — 1), we have 


X,@) =ar? (2n—1)a7C, 
and 


Y,(t) =ar?(2n—1) bfrE, (t) - 
If we let L[Q(4)] = g(p), we have 


(0) = d3 On— 1) arf bfeye tb afetal (35) 
n=1 


Using (31) and (32) in (35) we obtain 


2p ae 
q (Pp) = 16D) (2n— nF a (1 ers mel: (36) 
Since SP (2n — 1) = N2 and rN = R, we find that (36) becomes 


n=1 


q(p) = nR2bfoE (+2 4h ; 
(37) 
Qa 


~ Np [Py +Pw—-al 


—1 (a1) . 


a 


Using an identity in the theory of sued polynomials (Hobson, 
1931) we can perform the summation in (37). 
Since 


(2k + TF, (a1) = Piss (a1) — P,-1 (a1) ’ 


where the prime denotes differentiation with respect to the argument, we 
have 


> (2m —1)Py-i (a1) = 2 (2k +1) Pr (a1) 


(38) 


N—-1 


=D) (Pei (ar) — Pra (a1) ] +Po (a1) =P'y (a1) +P ys (a1) - 
k=1 


Defining G(¢) as the fraction of the substance remaining at time ¢, and 
letting L[G(t)] = g(p), we have, after substituting (38) in (37) and divid- 
ing the result by (33), 


¢1 \ Pw (a1) +Py—1 (a3) 
8 (P) =5- mote (14595 Py (a1) +Py-1 (a) ° ed 


The only factor of g(p) whose inverse transform is not obvious is that 
involving the Legendre polynomials. Here we need the following theorem: 
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if a function F(z) has a simple pole at = 2, andif u = u,isa simple root 
of the equation z(w) = z,, where z(z) is an analytic function of » at 
“u = up, then u = u, is also a simple pole of F[z(w)] considered as a func- 
tion of u. The proof of this theorem is given by Opatowski and Schmidt 
(loc. cit.). 

In order to apply the theorem we need to prove that the roots of 
Py(z) + Py_,(z) = 0 are distinct and thus are not roots of P (8) 
+Py_,(z) = 0. Let us assume that Py(z) + Py_,(z) = 0 has a multiple 
root. Then Px(z) + Py_,(z) = 0 for this root. For all z we have (Hobson, 
loc. cit.) 


(z2—1)Py(s) = N [Py (2) —Py-1(2)], 


(2?— 1) Py-i(s) = (NW—1) [#Py-1(%) —Py—a(2)1]. 
Adding these expressions, and using the identity VPy — (2N — 1)zPy_, 
+ (NV — 1)Py_, = 0, we find 
(2?— 1) [Py () +Py-1(2)] = (+1) N [Py (2) —Py-1(2)]. 


We conclude that the only possible common roots occur at z = — 1. But 
Px( — 1) + Py_,(-— 1) = N(— 1)*", while Py(— 1) + Py_,(- 1) = 0; 
hence there are no common roots. 


We have finally 
_,§Pw (a1) +Pw-1(a1)) _ ¥ Fb) me (40) 
de aay era leas ae 
where p, and p, are the roots of 
P+ Pf (ditdet+é[1— ww] ) +g (1—y) =9, (41) 
and y, are the roots of 
Py (2) +Py-1(2) =0. (42) 


The roots y, are obviously all real and lie in the range ae Varma ls 
Then the real parts of p, and p, are all negative, and the imaginary parts 


are zero if 
(ditdeté[1— Yn] )?> 4&¢2(1— Yn)» 


or 
eee ere Cy.) Ss 265 (43) 
We have 
2Dgl _2Dgl iv 
E(L—w) = Fagg, (1-0) = Rap, Nw) (44) 
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Using the relation (Hobson, Joc. cit.) 
2 
lim P,( 1-525) =Jo(2), 


we see that 2N%(1 — y,) > (2.405)”. Here 2.405 is the smallest root of 
J,(z) = 0. Upon analysis, we find that (43) is always true, and that p° 
and p, are real, negative numbers, for all positive ¢1, ¢2, and &. 

In order to express dp/da; conveniently we put 


tr =Patder:, 2 =Pr toe. (45) 
Then 
dp E 
=|} + =—-——_.. (46) 
dai! Pn 12% 


Then, inverting the expression for g(p) we find 


y 2 grt? 
ese eS ges = = 


Fi (47) 
grt? eet =e rs 
eas > Gea) ( rs ): 
Since rr = — digo and re +r, = di — do — &(1 — y,) we find 
N +2 
Tei wl é é 
2) Geetre) 3 2) at 
(1-4) n(i-Z 
Yn T» (48) 
. 3 E ie 
> he aaa li 
We also find 
N 
E rs 
SAE da! p* 7 = Sree r+ —r- 
(49) 
é fr Nv 1 
Te — $2 area 1—-¥y7, 
We may also prove that 
MoGQ 1 
2 ™ n=1 Se bag 
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We consider a polynomial of degree N with real coefficients, f(x), with 
real, distinct roots 1, 72,..., fn. Then identically in a, for all x ¥ n, 
fone, fh, we have 


Take f(x) = Py(x) + Py_,(x). Since P,(1) = 1 for all n, we have 


Peel) Py (1) (51) 


This is true since 2Py(1) = N(N + 1)Py(1). 
We may now write the final expression for G(#), using (46), (48), and 
(50) in (47). We have 
2 Edo = 1 
N? (1+ $2) 4 (p57 — pF) 


Baad) Pei (1) iter N? 
= ee 


G(t) = 
(52) 
Pa téoitds ot! pp téditds 97 
x 1 rr ay = it be i 5 
If we arrange the roots of Py(z) + Py_,(z) =0 in the order 71 
> ¥2 >... > Yy, then we find that the roots of (41) are in the order pr 
>p, >...>py>p, >p, >... > py. Furthermore p* > —¢», while 
p, < —¢2. Consequently for large ¢ the terms of G(¢) containing exp (9, #) 
are negligible compared to those with exp (p,1). 
If we assume that ¢: ~ ¢: and that £(1 — y,) > du, then p, ~ —¢n, 
and 


2d: Pn toitde 21 (53) 


N?(6:+¢2) (b, —Pp) bp N®0—%) (br 92) * 
For V = 5, ” = 1, this becomes 
0.700¢1. 
git do ’ 
while for V — ©, = 1, we have, instead of 0.700, the value 0.692. 
Therefore, for large t, ¢1 ~ ¢2, and &(1 — 2) > ¢1 


0.69261 set 54 
Aoi g1+ 2 Set nt OR? 
For small ¢ we have 
G(t) = eee ee eee , (55) 
N (¢i1+ $2) 


The coefficients ¢: and ¢2 appear here only because the distribution of the 
substance in the two compartments depends on these constants. 
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The total amount of substance inside the fiber and in the interstitial 
spaces can easily be calculated from (31) and (32), as can the amount in 
the fibers and interstitial spaces of each layer. 

The intra-fiber concentration of the first ring, that is, at the axis of the 
nerve, is given by 

a gin! 


bh ES sper ? 
p> Py iq) Op = Pad | 28) 


For large ¢, V = 5, we have 
ky 


Ci) = 1.59 7 Ee Feng (57) 
For NW ranging from 6 to 10 the numerical constant in (57) is 1.59 with an 
inaccuracy of less than one per cent. 

Discussion of the theory. The theory here developed, while approximate 
and somewhat complicated, may be of use in estimating some biologically 
important parameters in an objective manner. Used in conjunction with 
anatomical data on nerve trunks, a rapid and convenient way is given for 
the estimation of k; and k,. With somewhat more labor, an estimate of 
Dgl may be obtained. However, as mentioned previously, the nerve trunks 
without sheaths swell, and the Dg/ obtained may not be applicable to 
intact nerves. 

The use of the theory in connection with experiments in which various 
electrical characteristics are measured during the soaking-out process 
would seem to be particularly interesting. Here correlations between con- 
centrations of substances within the fibers at various distances from the 
nerve axis and the changes in electrical properties may be accurately esti- 
mated. 

The author wishes to thank Dr. I. Opatowski for inspiring this work and 
Dr. H. D. Landahl for a critical reading and discussion of the paper. 
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Certain parameters are defined which roughly characterize the internal structure of net- 
works. A given network structure uniquely determines the values of the parameters, but the 
reverse is not true. The parameters therefore define certain classes of networks. 

One of the parameters, the dispersion D(S) gives an indication of the “‘compactness” of the 
internal structure. 

Addition theorems and inequalities are derived relating the dispersions of sub-systems to the 
dispersion of the complete structure. 


Introduction. The terminology and definitions of this paper are consist- 
ent with those used in a previous one (Shimbel, 1951). It is assumed that 
the reader is familiar with the contents of that paper. 

Structural parameters. If we choose any site z in an adequate network 
structure, then, by definition, it is possible to reach any other site 7 by 
passing through a chain of directed links which connect them. 

In general there may be many different pathways from 7 to 7 of equal 
or unequal length. (Here length is taken to mean the number of links.) 

One or more of the pathways from 7 to/ will be the shortest, that is, will 
contain the least number of links. We shall call this number the distance 
from z to 7 and represent it by the symbol /(7, 7). 

Notice that since the links between the sites are asymmetric, /(z,7) may 
be different from /(7, 2). Also we will take /(7,7) to be zero for alli equal to/. 

The function /(i, 7), as we would expect of a metric, satisfies the tri- 
angular inequality. Thus 


L(i, j) SIG, k) +1, 9). (1) 


In discussing networks we would like to use descriptive terms which 
would apply to large classes of nets no matter how different they may ap- 
pear in their structural details. One such classifying property is the aver- 
age number of links per site. Another is the solution time T.. 

Consider the sum of all the /(i, 7) taken over 7 andj for a given network 
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S. If this sum is small, the network is “compact” in its connections. If this 
sum is large, the network is highly dispersed. bi: 
We shall call this sum D(S), the dispersion of S. By definition, then, 


n 


DS) =) ae (2) 
i=1 7=1 

Notice that D(S) isa sum ona sum. The first summation on /(z, 7) taken 

over j is also of interest. If this sum is large, then we know that on the 

average the various sites of the network are far removed from site 7. If it 

is small, then the various sites of S are readily reached from 7. We shall call 

this sum the accessibility of S to i and label it A(z, S). By definition, then, 


AG, 5) = 1G, 9) (3) 
and 
Dis) = AGS). (4) 
i=1 


On the other hand, if we sum /(i, 7) over 7, then we have a number 
which tells us how accessible 7 is to S. We shall label this sum A~(j, S). 
Note that 


D(S) = 34,5) = AG, S). (5) 


So far we have defined four parameters: 

1. 1(z, 7), the distance from 7 to j, 

2. A(i, S), the accessibility of S to 7, 

3. A(z, S), the accessibility of i to S, and 
4. D(S), the dispersion of S. 


The first of these is a relation between any two sites. The second and 
third are relations between any given site and the rest of the structure. 
The last is an over-all property of the network. 

Although every structure has a unique D(S), the converse is not true. 
In Figure 1 two different structures are shown, but both have a dispersion 
of 18. 

On the other hand, the totality of all the /(7,7) does uniquely determine 
the structure of the network. In fact, if we would arrange them in a square 
array, placing /(a, 6) in the ath row and bth column, if we then replaced all 
of the /(7,7) # 1 by zeros and put ones along the major diagonal, we would 
have a structure matrix of the network. The square array D(S) (we shall 
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call it the dispersion array) is therefore easily modified to give the struc- 
ture matrix. The reverse is somewhat more involved. 

In order to obtain the dispersion array D(S) from the structure matrix 
S, we proceed as follows. First we must find all of the status matrices Lo, 
th, ... Ir, where T is the solution time of the network. But from a the- 
orem previously derived (Shimbel, 1951), we know that these successive 
status matrices are simply the corresponding powers of the structure 
matrix S°, $1, S?,..., S57. The element in the ith row and jth column of 
the dispersion array is formed by adding all of the zeros in the ith row and 
jth column of all of the matrices S° through S?. 


Wz @ 


@ @ 
@ ©) 


@ ® 


Ficure 1 


That this procedure is valid follows from the fact that /(z, 7) is the short- 
est distance from z to 7 and the number of zeros in the 7th row and/th col- 
umn of the status matrices is exactly equal to the time it takes for a mes- 
sage to go from 7 to 7. These two numbers are of course equal. 

It might be noted that the sum of all of the elements of the dispersion 
array D(S) is equal to the dispersion of S, namely, D(S). 

Connection between networks. Suppose we are given two adequate net- 
works M and N, of order m and » respectively, and that we must join 
them by a single symmetric link so as to form a single network S of order 
(m +n). Let us suppose that one end of the symmetric link which joins 
the networks M and N is a* (located in M), and the other end is joined to 
b* (located in NV). Let the sites of M be denoted by a;, those of N by 8;. 

Let us first derive the quantity A(a:, S), that is, the accessibility of the 
combined network to a typical member of M. Now since A(ai, S) is 
simply the sum of all distances /(a;, x), with x ranging over the combined 
network, A (a;, S) contains A(a;, M). But in addition to A(a:, M ) we must 
find the sum of the path lengths from a, to the sites in V. However any 
path from a; to a typical site in N (say, 6) must pass through a* and b* 
since these are the only two sites linking the component networks M 


and NV. 
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Since every site in V can be reached by a site in M only by passing 
through b*, A(a;,.S) must contain A(6*, NV). For each such path we must 
pass first from a; to a* so that A(a,, 5) also contains ni(a;, a*). Finally 
we must add the integer to the total to account for the fact that the 
symmetric link must be crossed exactly # times in order to account for all 
of the paths from a; to V. This leads us to our first expression, namely, 


A (a;,S) = A(a;,M) +A (b*, N) +n1(a;, a*) +n. (6) 


Now if we sum the expression (6) over all a; we obtain 


S A (a; S) = D(M) +mA (b*, N) +nA-(a*,M) +mn. (7) 


t=1 


Without further argument we can write an expression completely sym- 
metric to (7), namely, 


>) 4 (bj, S) = D(N) +n (a*, M) + mA-'(b*, N) +mn. (8) 
j=1 


It can be seen from the definition that D(S) is simply the sum of the 
left-hand members of expressions (7) and (8). 
Thus we obtain 


DiS) =D (MY -+- DCA) me (Ae, Nye es ee) 
(9) 
+n [A (a*, M) + A~*(a*, M)] +2mn. 

Expression (9) relates the dispersion of a network made by joining two 
networks with a symmetric link to the dispersions of the constituent net- 
works and certain properties of the joining sites. 

Mean path length. If we define 1(M) to be the average /(i, 7) taken over 
all pairs (i, 7), then, by definition, 

m?|(M) = D(M). (10) 


This of course assumes that /(i, i) = 0 for ali 7. 
Furthermore, if we define /(i, M) as the average path from i to any other 


site in M and Hg. M) as the average path length from any site in M to 1, 
we have 


nl (i, M) = A (i, M) (11) 
and 


nl! (i, M) = A-1(i, M). (12) 
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Now by using expressions (10), (11), and (12) and substituting in equa- 
tion (9), we obtain 


(m+n) *1(S) = m2] (M) + n°(N) + mn [i(b*, N) +J-1(b*, N)] 


; a (13) 
+ mn [1(a*, M) +1-1(a*, M)] +2mn, 
from which it follows that 
i(s) = m* . cy 4 8 mn - 
1(S) Tm tnye WD) = (mia)? Tete) t {1 (a*, M) 
(14) 
+7(b*, N) +7 (a*, M) +1-1(0* _2mn 
( ) +i" (a*, M) +17 (6 SR sre eters? 


Expression (14) gives the mean path in the combined network as a func- 
tion of the mean paths in the constituent nets, the size of the constituent 
nets and the particular sites used to join them. 


FIGURE 2 


If short communication paths are considered desirable (and there are 
many reasons why they might be so considered), then the best way to join 
two networks by a single symmetric link would be to choose the joining 
sites in such a way as to minimize the right-hand member of (14). 

Note, however, that if the constituent networks M and WN are given, 
then, according to expression (14), the mean path length of the combined 
network can be minimized only by choosing the sites a* and 6* in such a 
way that the swm of their mean outgoing and mean incoming paths is a 
minimum. 

The central station. Suppose that a group of & adequate networks are 
joined together by symmetric links to a central station as indicated in 
Figure 2. Let the ith network be identified by the symbol NV; and let the 
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joining site of the ith network be labeled x. Also let 1; be the number of 
sites in N,. Finally, let the newly formed total network be called C. 

By reasoning similar to that which led to expression (9) we can arrive at 
an analogous equation for the combined network C, namely, 


k 
D(C) = >) D(M™) 
i=1 
kk 
+. a >; O50: [A (zs N;) = A-1(3}, N;)] 
j=1 i=1 
k k ok 
+2S°2:+4>5 >» 6 :jN iN; 
i=l ci eof 
k 


+ >* [A (at, Ni) +A (x7, N)1- 


t=1 


(15) 


Now it follows from expression (14) that 
k 2 


ic) = > — 1) 


“E, 


k 
+ >> > 6,;—** Ind (FN) +20d7 (xf, ND] 
i=1 j=1 
(2 (16) 
k = k —gk Aone 
+2 Gs ni) BD na) p> Du Bastian 
i= a= i=1 i= 


k ceooak e 
+ (Sma) SO tat, we) tude, 1. 
a=1 i=1 


ieee expressions (15) and (16) the symbol 6,; is the familiar Kronecker 
elta. 

Expression (16) gives the mean path length of the combined network C 
as a function of the mean path lengths of the constituent networks, their 
sizes, and the connecting points. 

One point of interest is that the minimization of the mean path length of 
C can be achieved only by the same process as that used in minimizing the 
mean path length of S. That is, for all of the connecting sites «*, the 
quantities " 


al U(x, Ni) +17 (aF, Ns) 
must be minimum. 


COMMUNICATION NETWORKS 507 


The triangle inequality. Returning now to expression (1), note that if we 
sum over all the sites 7 of the network S, we obtain the expression 


A*(S, 7) SA*(S, k) +nl(k, j). (Ti) 
Now if we sum expression (17) over the sites 7 we obtain the inequality 


D(S) Sn[A1(k, S) + A (k, S)]. (18) 


An alternative form of (17) is 


1(S) SI-'(k, S) +1(k, S). (19) 


The inequalities (18) and (19) can be easily substituted in equations 
(9), (14), (15), and (16), thus obtaining inequalities which contain only the 
dispersions of the constituent networks and their sizes or only mean path 
lengths. 

These inequalities would put lower bounds on the dispersions and mean 
path lengths of the connected systems. These estimates, however, can 
probably be greatly improved. 

Other parameters. The parameters used so far say little or nothing about 
the amount of “‘work’”’ done by each site in a communication net. 

Suppose that in order for site z to contact site, site k must be used as an 
intermediate station. Site & in such a network has a certain “responsibil- 
ity” to sites z and 7. 

If we count ail of the minimum paths which pass through site k, then 
we have a measure of the “‘stress”’ which site & must undergo during the 
activity of the network. A vector giving this number for each number of 
the network would give us a good idea of stress conditions throughout the 
system. 

It would be interesting to investigate this vector for a combined net- 
work in a manner similar to the analysis of mean path length and disper- 
sion carried out in the foregoing. 

This “stress” vector can be derived directly from the structure matrix 
but unfortunately the process is not simple. 

This investigation is part of the work done under Contract No. AF 
19(122)-161 between the U.S. Air Force Cambridge Research Laboratories 
and the University of Chicago. 


LITERATURE 


Shimbel, A. 1951. ‘‘Applications of Matrix Algebra to Communication Nets.” Bull. Math. 


Biophysics, 13, 165-78. 
iophysics, RECEIVED 2-17-52 


m5) 
ari 
is 
4 
4 


ite: 7 

in 

.* i) 

ful 
He? 
we 
Lan 
nf 


7 
au 
iv 
4 = : | 
| | 
nin 
« 
> 
ri 
" 
sonata? 
a) 


ra < b 
a 4 
arty 
—? 
: as 
or) i 
a 
4 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 15, 1953 


ENZYME LOCALIZATION AS A MECHANISM OF APPARENT 
ACTIVE TRANSPORT 
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THE UNTVERSITY OF CHICAGO 


A model based on enzyme localization is developed which gives rise to an apparent active 
transport of a metabolite into or out of cells. The model is applied to three simple situations, 
using Fick’s equation and the Rashevsky approximation. It is shown that the apparent ef- 
ficiency can be made as large as desired if, for constant reaction, the outer cell region is made 
sufficiently small, or, for autocatalytic reaction, if the metabolite concentration in the outer 
region is sufficiently small. The physical limitations imposed by this mechanism are developed 
for all three situations. 


It is now generally accepted that concentration differences of metabo- 
lites between cell interior and exterior cannot be accounted for solely in 
terms of permeability properties of the membrane but frequently involve 
active transport against the concentration gradient on the part of the 
whole living cell (Steinbach, 1951; Teorell, 1949; Sacks, 1948). Such active 
transport involves, of course, the continuous expenditure of free energy by 
the cell; the transport mechanism must therefore be coupled to the cellular 
metabolic system (Rosenberg, 1948; Spiegelman and Reiner, 1942). A 
model of such a coupling, utilizing osmotic pressure, has been described by 
J. Franck and J. E. Mayer (1947). 

Some of the probable characteristics of such a transport system have 
been described in the literature (Ussing, 1949). They include: 


1. A carrier capable of forming a complex. 
2. A system whereby the metabolite is set free from the complex. 
3. Spatial separation of the components of this system. 


In this paper we shall set up a model whose main feature is the localiza- 
tion of enzymes within the cell. Such localization actually occurs in cells; 
for example, the well-known localization of the oxidative enzymes in the 
mitochondria or of phosphatases in the cell membrane region. Such locali- 
zation, under the proper conditions, will give rise to a concentration dis- 
tribution within the cell which will show all the characteristics of “active 
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transport,” i.e., metabolites will seemingly flow against the concentration 
gradient as long as we compare the average intracellular concentration 
with the outside concentration. Of course, in biological experiments it is 
exactly this cell average which is determined, and on the basis of which 
data concerning ‘‘active transport” are calculated. 

The model. We postulate a cell consisting of two concentric regions, in 
each of which a different enzyme is localized. A metabolite diffuses into 
and out of the cell; it undergoes opposite reactions in these two regions, 
e.g., it is broken down in the outer region and produced in the inner one. 
We shall assume that these two regions are covered by semi-permeable 
membranes with permeability coefficients / and /2. We shall also assume 
an infinite medium, such that the concentration of the metabolite is a 
constant (Co) at infinity. 

Such a system can very easily give rise to apparent active transport. 
Let us assume for example that the metabolite is produced in the outer 
region and consumed in the inner one. If the steady state concentration 
gradient at the outer membrane is positive, the metabolite will diffuse out 
into the medium. If the consumption of the metabolite in the inner region 
is so large that the average concentration over both regions is smaller than 
the outside concentration, we have a situation where a substance seem- 
ingly flows out of a region of lower concentration into a region of higher 
concentration. Actually, of course, flow here is always along the concentra- 
tion gradient, but there is such a distribution of sources and sinks that the 
effective flow over the average of the cell is against the gradient. The preser- 
vation of some of these sources and sinks is dependent upon a continual 
supply of free energy by the cell, and is therefore directly coupled to the 
metabolic system. Figure 1 gives a simple schematic picture of this situa- 
tion, where C, and C; represent the average concentrations in the two 
regions, Cy) the average concentration in the medium, and the dotted line 
represents the average concentration over the whole cell. As can be seen, 
this average is lower than the external concentration although there will be 
a continuous flow of metabolite from the higher concentration in region JJ 
to the lower one in the medium. 

An exact treatment of constant reaction, We shall first treat the simplest 
possible case, that of constant reaction rate. The physical picture is as fol- 
lows. We shall assume a spherical cell of radius 72 consisting of two con- 
centric shells. In the inner shell of radius 7» an enzyme is localized catalyz- 
ing the reaction A = B; the forward reaction rate is K 1, a constant; the 
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reverse rate is so small that it can be ignored. In the outer shell another 
enzyme, for example, a phosphatase, is localized, catalyzing an energy- 
liberating reaction. The energy liberated here is then coupled to the re- 
verse reaction B = A at a constant rate Ke. Therefore the total result of 
this enzyme localization is that in the inner region, region J, A goes to B 
at the rate Kj; in the outer region, 77, B goes to A at the rate Ky. The 
preservation of this cycle is dependent upon the continual supply of free 
energy in region JJ by the cell, and this can be viewed as the “‘cost”’ of this 
process. 

Mathematically the model is extremely simple, especially if we omit 
permeability considerations. We shall also assume that B has a constant 
concentration Cy at infinity and A a concentration zero. 


FIcureE 1 


We shall use the following notation: Cy is the concentration of A in 
region J, Ci is its concentration in region JJ, and Cj; its concentration in 
the medium; Cx, Cx, and C23 are the respective concentrations of B. If we 
now define X as rC, the diffusion equations then become (Rashevsky, 
1948) 


2 a?X 
Dien tr Ki 0, ieee Lh 0, 
d2X d? Xo = 
- oa 0: Dear RF (1) 
d?X 99 D aXr3_ 
Recrergar bam at a 
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The boundary conditions are as follows 


Cu (0) Fo, Co (0) Fo, 

iim Cra. lim Co3 =Co ) 

Cu (ro) =Cr (ro); Cio (72) =Ci3 (12) 

52) a. a 52) = a) ( 2 ) 
OT Ji ONaE  ee dr /,, dr /,,’ 

Cai (70) =C22 (10) Coo (r2) =Cos (172); 


G2). -($*) s2) = a2) 
“ dr /;, dr /,, 


The solutions of equations (1) are 


ie 


Ci = me ee ae 
K 
Cu=t+ 02, Cu = Bi + 72 (3) 


K 2 
: #25 Cu=4m. 


An evaluation of the various coefficients of (3) by means of (2) finally 
gives us the following expression 


Kz 2 ro Ky 
Cu = apy (a 7 opr 
2 3 
re Ko 2 fo 1 
Seely Fie MG Mates AGIA Doe 


3 
Ke. 3 To 


Hf 
Cu=—| 53 CE ee me 
3D 3D 
r : 1 (4) 


K 
Cs — 3p," ate 
r 
RO pa tet iah ee ats Ks 
2 
1 Ky : To 
Co3 = Cy — — aa 
23 0 = Lak? eas, 
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Since we must require that the concentrations of both A and B remain 
non-negative for all 7, the following inequalities are imposed 


3 3 
Ki< Ke (re ro) 


ce : 
K,< F262 = 10) (5) 
r? : 
0 
Ko cee Kiro 
Pee aay, (re ro) wre 


We shall now assume that 4 is some kind of a complex containing B. 
Therefore in the calculation of the concentration of B we must also include 
the concentration of A. The flux across the cell membrane is therefore 
given by the sum of the gradients of Cy and Cy at r,. From (4) this can 
be easily calculated. It is 


ic dCn\ _ 
T=D:G*) +Di G2) =0. (6) 


Hence we find that the flux is zero, a result which was to be expected, con- 
sidering that we are dealing with a cell in the steady state. 

To calculate the average cellular concentration C, we simply integrate 
the various concentrations of A and B over the respective regions and 
divide by the total cell volume. The final value for ¢ is then 


C=C,— 


D2—D,{ Ker3rs | Kirsro 2 5 Kiro eg 
eA bert Sa ee eT 
ADD, | PAE Sie ace AE elon ripen ee 


As can be seen from (7), for ro = 0, € equals 


275Ke(D2— Dy) | 
Cot 5D, Dz, ? 


for ro = 1, C equals 


Hence for all ro > r*, 0 < r* < ro we preserve a smaller intracellular con- 
centration, while for all 7. < r*, 0 < r* < 7 we preserve a larger intracel- 
lular concentration than in the outside. The specific value of r* is deter- 
mined by the parameters of the cell, i.e., Ki, Ke, Di, and De. 

It is of interest to calculate the efficiency of this apparent active trans- 
port. We know what the cost of the process is: It is simply the amount of 
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free energy expended by the cell to keep the reaction B > A proceeding at 
the constant rate K. in region JJ. The problem arises as to what shall be 
defined to be the work derived from this process. Rather arbitrarily we 
shall assume here that this work (or, rather, apparent work) is the differ- 
ence between free energy preserved between the inside of the cell and the 
medium, and we shall take € as the inside concentration and C» as the 
outside concentration. The amount of work done is, therefore, simply 


' a 
TT? [woaCo — HiaC] 


The efficiency of the process is then the following 


_73{ [Hoa ERT In Co] Com (oa +R In C)C} (8) 
iia Ke ls oe 


An inspection of (8) immediately discloses that 7 goes to infinity as ro 
approaches 72. Hence we get the seemingly impossible result that the 
efficiency can be made as large as we please as long as the inner shell is 
made sufficiently large. Actually there is nothing wrong with this conclu- 
sion, as long as we remember that the efficiency is only apparent here, in 
the same sense as the active transport here exemplified is only apparent 
active transport. The chemical reaction does not actually drive substance 
B against its concentration gradient; as a matter of fact, as pointed out, 
all substances flow along their gradients. But what we do find is that asa 
result of a certain distribution of sources and sinks and a certain expendi- 
ture of free energy—even if only a very small one (i.e., ro very close to 72)— 
we can get an average concentration inside the cell which is different from 
that in the medium. And insofar as that suffices for the requirements of the 
cell, we then have here an extremely cheap mechanism whereby the cell 
can accomplish such differences. 

It is interesting to remember that phosphatases which seem to be in- 
volved in active transport are localized in the membrane region, i.e., our 
region JI. According to our mechanism such a localization would then give 
rise to a very large efficiency. Of course it must be pointed out that we 
have not been considering the cost of preservation of the two regions, i.e., 
the preservation of the localization of the enzymes. Such a cost cannot be 
assessed here. 

Approximation equations for two concentric regions. An attempt to apply 
this model to somewhat more complicated cases gives rise to mathematical 
difficulties if we insist on an exact solution. Instead we shall use Rashev- 
sky’s approximation method. It is assumed here that the reader is familiar 


ENZYME LOCALIZATION 515 


with the derivation of the approximation equation (Rashevsky, 1948, pp. 
16-23). Applying those principles to the case of two concentric regions, we 
find the approximation equation for the inner region to be 


ee 3D: Ii(Ci-C2) _ 6 Din (C1 —Cr) (9) 
dt ** "ry (2DiFiars) 11 (2Di-F Iara) | 


In the outer region it is 


aes 3 ri Dik Ci — Cs) 6 r2Dirihy (Ci — Cs) 
dt jah t7 5) 2) 7h) WP eine Ae (2D --- her) 


i. 


= 37D; D,he (C2 —Co) BOO WEDD 
et. tert) PD Ait. 2%.) +2 DADs oh] 


es 67374.D:0. 8: (Co —G) pot 
eee er) rr, 2 DD oh) | 


where g; and gp are the respective reaction rates, C, and C; the respective 
concentrations, 7; and 7 the dimensions of the inner region, and 73 and rg 
those of the whole region (for detail, see Fig. 2). 

To solve for the steady state, we set (9) and (10) equal to zero. We then 
solve for C; and C2 separately, and after some tedious calculations find the 
following 


2 
FES g2 1112 qi 11 
C2 te Ser: ene OB 
and 
= Ar® (a:+ a2) + rire (Bi + Bo) q2 
ie Te Oe ae aS Ie Co 
Gel anatase ae a eee 
where 
a 377D;D.he Sr 
1 Ar? [Dehe (rs — 2) + 2D; (D. + 5h2) 1’ 
a 67374D;Dhe 
*2 °° Br [Delta (14 = 71) + 2D: (De + he) 1’ (12) 
B = 3D hy 
"re [2Di+ rshi) ’ 
6 Diks Ar® = rata _ or; . 


B= (2D: ral)’ 
Equations (11) and (12) provide us with the equations for the case of two 
concentric regions. 
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Since we have already treated the situation of constant reaction rate by 
using Fick’s equation, we shall apply the above approximation expressions 
to two cases: 1) autocatalytic reaction and 2) a second-order reaction. 

Autocatalytic reaction. We assume here a similar model as in the previous 
case, except that E is broken down to Bin region J at a rate proportional to 
its own concentration, and re-formed in region JJ also at a rate propor- 
tional to its concentration. To be specific, g: equals — aCy and g@ = aC. 
For B we have gq = “Cu and gq = —a@Cy. We also assume that the 
medium concentrations of E and B are Cy and By respectively. 

To keep the model somewhat general, we shall postulate that Ci: is kept 
at the constant concentration CyA, where the sign and magnitude of A 


Ficure 2. The dimensions of the cell regions as shown are denoted by ny, 72; 73, rs; Ci is 
the average concentration in J, C2 the average in JJ, Cio the average in the medium, h; the 
permeability coefficient between J and IJ, he between JJ and the medium, D; the diffusion 
coefficient in J and JT, D, the diffusion coefficient in the medium, and 6 is a length of approxi- 
mately the dimensions of one of the 7’s. 


are not prescribed. In other words, we begin by assuming a constant flux 
of unspecified direction. According to (11) it follows that 


Cy = Sot A (aa = 01 = a2) rie 


ya ’ sore oe (13) 


But Cy is also [from the second equation of (11)] 
ae Y qe 
cf Freteryisesenn na mulaiaecia 


Therefore, solving for a2 between these two equations gives us 


Be a [y (Bi + Bz) (A+C)) + A (ai+a2)] +A (ai+ as) (Bi + Bo) 
RO (FS TE EES hw 


If we substitute (14) into (13) we finally find the following expression 
for Cy 


as 


_ (Bi +82) (A +0) 
Coe aro 645) 
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We can now apply (15) and our postulated value for Cis to calculate the 
values for Cx; and Cz. Again using (11) we find after some calculation 


Jeo eat oe 
Co. = Bo — A e 


es a; (Co — AB, — A Bo) 
(16) 


If we assume that B forms a complex with £, we must then include B’s 
concentration in the calculation of the average C. Under those conditions, 
the average is simply 


1 


2 
137% 


C= [ (Cir +Car) rire + Ar’ (Ci2+Cop) J (17) 


which, for our values, gives us 


rife A (B,+ Be) (1 — a) 


2 od Dee Tee en 


If we again assume that the free energy is expended by the cell in driv- 
ing the reaction B — £ in region JJ and define efficiency as before, the 
following expression results 


_ rari C (wor +RT InC) —Co (woe + RT In Cr) ] 
7 Artd,(Cot+ A) (uz— we) ; 


(19) 


One interesting fact immediately emerges: 7 does not approach infinity as 
the outer region goes to zero. Instead n approaches a finite value 


_7s7i(C (wo +RT InC) —Co (uot RT In Cr) ] 


rari (C (uo +RT In C) —Co(uo+RT In Co) ] ‘ 
ae Gaia es) O ; (20) 
where 
_ ay [irs (Bit B2) (A +Co) + A (aitas)] +A (01+ a2) (Bi +62) 
eee ETE eN TE PET en oe RAS 


and a; = Arai, a, = Ar'ap. 

We can see that this is so because a2 in (14) approaches ~ as Ar? be- 
comes zero. Hence in the autocatalytic case, as distinct from the case of con- 
stant reaction, the efficiency does not become arbitrarily large as the outer re- 
gion becomes infinitesimal. 

We also have to consider the flux. The flux is here simply proportional 
to Cy + Coo = Bo + Co minus the medium concentration, which is also 
By + Co. Hence we again find, as in the case of constant reaction, that the 
flux across the membrane is always zero regardless of the value of A. This 


518 ARTHUR BIERMAN 


result of course holds only if we include both B and £ and is then a neces- 
sary consequence of the imposed steady-state condition. 

We shall now consider the various possibilities in evaluating (18) and 
its efficiency. 

1) Case A > 0. Then, using (16), Bo must be positive and greater than 
A. Also, depending on whether a; is greater or smaller than 1, C is less or 
more than the medium concentration Cy + Bo. While 7 is stivisacly not 
limited to 1, but depends on the value of the various parameters, there is 
no physically reasonable parameter, which, by assuming large or small 
values, can make 7 arbitrarily large. For example, if wz — us, would be- 
come extremely small, 7 would become infinite, according to (19). On the 
other hand, ug — ws must be sufficiently large to prevent the reaction 
B — E from occurring spontaneously in region 7. Hence a physical limit is 
placed upon this parameter by the physical situation. 

2) Case A < 0. Here By can be zero, but —A must be less than Co. 
If a, is greater than 1, then (18) can be written as follows 


(a, — 1) (— A) (Bi+ Bs) rife (21) 


nes Gerke ata en een) eae 


Here then C is greater than C, and is preserved because of the zero flux. 
It can be seen now from (19) that 7 can be made arbitrarily large as Cy 
becomes very small or as —A approaches Co. Hence in the autocatalytic 
case the efficiency becomes very large if the concentration of region II is made 
very small. 

If B is not a complex, its concentration is ignored in the calculation of 
C and 7. The flux is now proportional to A. If A > 0, the metabolite flows 
out of the cell; if A < 0 it flows into the cell. The average concentration 
inside the cell is now 


rired; (Cot A) 


Cal a Sb Bal Bia 
Fe bE ol 


(22) 


1) For A > 0, the flux is out of the cell and hence C must be smaller 
than Co. From (22) this requires that 


2 
117 20,Co 


0: ge ee eR he ah Pe ee 
7,74 (8, +B, + a@,) a ri? 5, 


(23) 


Hence » cannot become arbitrarily large. 
2) For A < 0, the flux is into the cell and C must be greater than Co. 
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But we can see from (22) that this is not possible. Here apparent active 
transport cannot even arise. 

To summarize: Concentration differences are preserved by this model 
in the case of autocatalytic reactions. In particular, if B is a complex and 
Ciz = Co — A’, A’ > 0, the efficiency can be made very large for small 
(Orn 

A second-order chemical reaction. The last case which we shall treat here 
involves a carrier metabolite model whereby the two substances react by a 
second-order reaction. Specifically, we shall assume a transported sub- 
stance S and a carrier R which form the complex (RS) in region J, the 
interior region, and come apart in region JJ to give us R and S separately. 
We will also assume that S is present in the medium at the constant con- 
centration S» while R is not preserved in the medium. Instead, R is being 
constantly replenished inside region J at a constant rate w. In terms of 
actual substances S might be Na and & a polyphosphate molecule, a co- 
enzyme, ATP, or a nucleic acid molecule which is synthesized by the cell 
in a certain region, e.g., the microsomes. There is some likelihood that 
some of these polyphosphates are of crucial importance in ion exchange in 
the cell; RNA, for example, has been shown to be a significant factor in the 
calcium accumulation of cells (Steinbach, 1952). 

As stated before, we take the chemical reactions to be of second order, 
i.e., the rate of formation of (RS) in I is y:(S)(R) and the rate of dis- 
appearance of (RS) in IJ is y2(R)(S). From (11) we then find the following 
equations 

(RS)? = yi (R) 7 (S) Pir — y2 (RS) “Poe ; 


(RS) ¥ = yi (R)7(S) Pax — 2 (RS) "Par » 
(S)7 = y2 (RS) VPs — 1 (S)7 (R) Past So 5 


(24) 
(S) = ye (RS) UP2s — ¥1 (S) 7 (R) Past So , 
(R)7 = 2 (RS) "Por — "1 (R)7(S)"PirtwPrie » 
(R)¥ = 2 QS) Pen = 91 (R)7 (S) 7Psr+wPsr 
where 
ee 1 = riers 
Bit Bs | Ar8 (ar-az) ’ Oy 
1 es rire 
ihc tey aS” 


and w is a positive constant, indicating that R is being supplied at a fixed 
rate to the cell from an unspecified source. 
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To simplify the problem we impose the condition of zero flux upon the 
system, i.e., we let S// = So, where the Roman numeral indicates the cell 
region. This results in 


yo _ (S)7 (R) Pas 
‘1 (SR) “Pos j Bo 
Substituting (26) in (24) results in 
(RS)\' = (5) (pyr bie PenP ss 
Pos : 
FLAC EL 
Pog c 
(S)7 = yi (S)7 (R)? [Pas —Pis] + So , (27) 
(R)P= ys (R)(S)! gE 4 win, 
(R)# = 9, (S)"(R)! FEE tw Pan 
It follows immediately from (25) that 
P3_ rife Par—PorPas _ ¢ 
P» Ar Pos * on 
Pia PirPase 1 ay Toad 
Pigs et Bl bse oe ee a 
Hence (27) gives us 
x Pg 
a5) 4B, gynag 
I 
(S)t = Sy eee (S) = Sy, 
2 (28) 
(R)i = rife ek a (R)7 (S)7 
hy lactectam moon Birt Bor ’ 
(R) i mene rife 1 


Ar® air+aer 


Solving for RX results in 


1 rir 
Pye ae 1? 2 Bir+B 
Laie Eresomerciers: @2R) laste orl <2?) 


We can now use (29) to solve for (S)!, which we can do from the following 
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quadratic equation 


: Birt Bor w wrir: D 
gyi sy) 172 eee 
( ) ce ) l yt TBs pp are Dis S| 


_ 5, (Biz + Ben) _ 
yal 


considering that 
Bir ie Bor ms Dir 
Bist Bes = Dis’ 


if we let kh; = ~ 

Rather than actually solve (30) we shall only consider the physical con- 
ditions imposed here by the appearance of active transport. This requires 
that the average concentration of (S) over the whole cell be different from 
So, although no flux takes place. Therefore in the calculation of (§) we 
must consider not only the concentration of free S but also that present in 
bound form, i.e., (RS). This means that 


1 
ace asf rari [(S)7+ (RS)2] +Ar* [(S)#+ (RS)#] } 2 So. (31) 

3.4 
We established previously that (RS)! = 0 and (S)”¥ = Sy. It therefore 
only remains to calculate (RS)/. This we can do easily from (28). It gives 


us 
SRI BS Bare dt fs DE 
ca ieee CD, 


This in turn tells us that Sy) must be greater than (S)’. Introducing (32) 
into (31) gives us finally 


(RS)7= Som ee (s2) 


DS 1554 = (sit 0. (33) 
Dey 
This shows that § > So for 
Pin Shs | 
1g ’ 
and § < S> for 
Dir 
1. 
Deo 


Since R is a carrier, we can assume that 


Dir 


a Me 
Dis 


Hence S > So. 
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To summarize: We have developed a simple model based on enzyme 
localization which gives rise to an apparent active transport. Under proper 
conditions the “efficiency” of this model can be made very large, as long 
as we ignore the work done in preserving the localization. We have shown 
that this model works for three simple situations: 1) constant reaction, 
2) autocatalytic reaction, and 3) a second-order reaction. For all three, 
certain limitations were imposed upon the cellular parameters, but none 
of them of such a nature as to rule the model out. 


The author wishes to thank Drs. H. D. Landahl and G. Karreman for 
their many helpful suggestions in the preparation of this paper. 
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A previously derived iteration formula for a random net was applied to some data on the 
spread of information through a population. It was found that if the axon density (the only 
free parameter in the formula) is determined by the first pair of experimental values, the 
predicted spread is much more rapid than the observed one. If the successive values of the 
“apparent axon density” are calculated from the successive experimental values, it is noticed 
that this quantity at first suffers a sharp drop from an initial high value to its lowest value and 
then gradually ‘‘recovers.” 

An attempt is made to account for this behavior of the apparent axon density in terms of 
the ‘‘assumption of transitivity,” based on a certain socio-structural bias, namely, that the 
likely contacts of two individuals who themselves have been in contact are expected to be 
strongly overlapping. The assumption of transitivity leads to a drop in the apparent axon 
density from an arbitrary initial value to the vicinity of unity (if the actual axon density is not 
too small). However, the ‘‘recovery” is not accounted for, and thus the predicted spread turns 
out to be slower than the observed. 


The apparent axon density. The assumption of a completely mixed popu- 
lation implies that contacts between all pairs of individuals are equi- 
probable at all times. It was shown in previous papers (e.g., Rapoport, 
1951) that under this assumption the spread of state in “ordinal time” or 
“by removes” will be given by the following iteration formula 


Ptt+- 1s [1—2{4)] (le FO], Gt) 
which can also be written 
[1—x (¢+1)] ew =1-—-x(0). @) 


Here ¢ takes discrete integral values and represents the ordinal time, 
that is, the number of “removes” or times the information was trans- 
mitted, P(é) is the fraction of mew knowers at the ‘th remove, 


x(t) = >) Pj) 


7=0 


BYR) 
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is the total fraction of knowers at the ‘th remove, and a the “‘actual axon 
density,” is the number of tellings per knower during the entire process, 
assumed independent of ¢. We shall continue to use the neural net term 
“axon density” in this sense, which suggests that a “axons” or contacts 
‘Gssue” from each individual (on the average), and that this number is 
characteristic for a particular process. 

In-an experiment conducted by the Washington Public Opinion Labora- 
tory (1952), 33 simple messages were made to diffuse through a population 
of 184 school children under contest conditions. It was possible to trace 
through how many hands (removes) each message has gone. The experi- 
ment can be thought of as one in which a single message diffuses through a 
population of 6072 individuals, a population large enough for the applica- 
tion of our probabilistic method. Accordingly equation (1) was tested by 


TABLE I 


t 0 1 2 3 4 5 6 7 8 9 10 


PiU)icales 2030" 191 9) S85) 73191) 20025 006 

P(j) obs: .030_<:191) 191-156.) 1076 = 063.) 2035) 2 OSes 00S se 0U seme UOe 
x(t) calc. .030 .221 .806 .997 .999 .999 

“(t)lobs.. ©.030" 221 2412" 2568: 615" 738-173. (Ol S02 CUS nad 
a(t) hga ate Ai OU Gh) SD I8L S22005 9228 2 © 2 30 s206- Tor 2290 Raieay 


the data of that experiment. The value of P(0) = x(0) was fixed by the 
conditions of the experiment: each child was the starter of a message. 
Hence P(0) was fixed at 184/6072 = .03. The value of a was obtained 
from the observed value of P(1). The predicted ordinal time course could 
thus be compared with the experimentally observed one. This comparison 
appears in the first four rows of Table I. 

Obviously equation (1) predicts an ordinal time course of a spread 
which is much more rapid than the observed one. Attempts to account for 
this discrepancy led us to examine the so-called “apparent axon density” 
of the experimental process. The apparent axon density of any real ordinal 
time course of the spread of information is defined as the function a(é), 
which, when substituted for the actual axon density as in equation (1), 
will account for all the points in the ordinal time course. In other words, 
we shall describe any ordinal time course by the equation 


“P+1) = [1—%()] [1— e=O PO], (3) 


where a(#) is a function to be determined. Obviously if all the P(#), and 
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therefore the «(#), of the process are known, a(é) can be determined by 
solving for it in equation (3). We have, in fact, 


1 1—x( 
© 0 = 5G Me eaE i 

The fundamental assumption of the completely random process is 
equivalent to the assumption that the expression on the right side of (4) 
is a constant (i.e., the actual axon density a). Therefore, given any actual 
process, the constancy of the expression on the right of (4) is a test of the 
random net hypothesis. 

The determination of a(/) in the experiment referred to gave a set of 
values which appears in the last row of Table I. We observe that a(t) 
starts at a rather high value (7.2) at ¢ = 0, drops sharply for ¢ = 1 to its 
lowest value, then rises steadily for several successive values. For ¢ > 7, 
a(t) becomes erratic. However, as is seen from Table I, P(¢) becomes ex- 
tremely small for ¢ > 7, and, as can be deduced from (4), a(#) becomes 
excessively sensitive to small fluctuations in P(t) when x(é) is close to 
unity. It follows that the values of a(t) in that range are not reliable. At 
any rate, it appears reasonable to assume that the behavior of a(#) in the 
experiment at hand is characterized by a large initial drop and a subse- 
quent steady “‘recovery” in its significant range. 

One could attempt to explain this behavior of a(t) by psychological con- 
siderations. For example, one could attribute the initial high value of a to 
a “start effect,” that is, the enthusiasm of initiating the process. If, as is 
the case in some message diffusion experiments, only a portion of the pop- 
ulation are the “starters,” they can be expected to be more strongly 
motivated to spread the information than those who get the information 
“second hand.” This can conceivably explain the initial drop in a(¢). In the 
case considered, however, this explanation is not convincing, inasmuch as 
every one of the 184 children started one of the 33 messages. The same 
individuals, therefore, were involved in all the removes. It is hard to be- 
lieve that the motivation of the same individuals fell so sharply immedi- 
ately after the start that the axon density was reduced by a factor of 5. 

In addition, we are still faced with the opposite effect, namely, the un- 
mistakable steady rise of the apparent axon density. One could, of course, 
also explain this effect by psychological considerations The experiment 
was conducted under contest conditions, where prizes were offered both 
for giving and for receiving as much information as possible. It is con- 
ceivable that as the population became more and more saturated with 
knowers, and, therefore, as it became more and more difficult to find non- 
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knowers to give information to, the total number of contacts tended to 
increase with each remove. 

In this paper we will not make any psychological hypotheses. We will 
try instead to explain at least the initial drop in a(t) by some assumptions 
concerning the structure of the population, which places certain constraints 
on the possible contacts, so that they are no longer equiprobable. 

The finite acquaintance circle. The first constraint will be contained in 
the assumption that the contacts of each individual take place only within — 
his acquaintance circle. We will assume first that an individual has on the 
average g acquaintances randomly chosen from the population. We will 
assume, moreover, that g is large compared to a. Following the reasoning 
in the derivation of equation (1) (Rapoport, Joc. cit.) it is easy to see that 
under the modified assumption of the finite acquaintance circle we shall 
have 


PGA Stee [1-G-2)"] (5) 


If q is fairly large, so that the power in the second bracket on the right 
of (5) can be approximated by an exponential, (5) reduces to (1), and no 
essential modification has been introduced. If, however, g, although large 
compared to a, is small compared to NV, we can reason along a different 
line. 

Let us fix our attention on an arbitrary individual A at the tracing of 
the (¢ + 1)th remove. We seek the probability that on that remove A does 
not receive the information from an arbitrarily chosen individual B among 
the q individuals in his acquaintance circle. This can happen in either of 
two mutually exclusive ways: either B did not become a knower on the 
ith remove or he did become a knower on the ¢th remove, but his contacts 
among his own acquaintances did not include A. The probability we seek 
is the sum of the probabilities of these two events, that is, 


1—P(i) +P (4 --)'. (6) 


Now if all the states of the acquaintances of A are independent of each 
other, and if g is small compared to the total population, so that sampling 
with replacement can be assumed for any sample of individuals not 
greater than q, then the probability that A did not receive the information 
on the (¢ + 1)th remove will be given by 


q 


{1 -P( +P(o(1-2)t (7) 
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Then the probability that A did receive the information on the (¢-+ 1)th 


remove will be 


1—J1-P(y +P (1-2) Pa1-f1-Pw mh, (8) 


m=1—-(1-2)", (9) 


Expression (8) corresponds to the expression 1 — e~?? in the completely 
random case. We therefore write for our modified equation representing 
the spread of information by removes 


Fei) = [12x (4) ] (1 — {1 —P) mm}. (10) 
Solving for a(#), as defined by (4), we have 


where 


a (£) = py og lt —P Cs) m). (11) 


Since P(t) < 1, and m < 1, we can expand the right side of (11) to 
obtain the series 


gm+ 3qm*P (t) +iqm'P (t)2+.... (12) 
We note further that if g is large compared to a 
bo oe 
7 (13) 


so that 
(14) 


2 
a(t) = of 1+ 5 + aie ein) 
Since P(t) < 1, and a/g «1, the series converge rapidly. Thus a(f) de- 
pends (but very weakly) upon P(d), rising and falling with it. As qg be- 
comes very large a(#) tends to a, as, of course, should be the case and as is 
clear from (5). 

Although this approach to the finite acquaintance circle case still gives 
no substantially new result, it does lend itself rather readily to the imposi- 
tion of a socio-structural bias, which we will now discuss. 

The dependence of probabilities. So far the assumption underlying the 
whole argument was that the probabilities, 1 — P(t)m, that each of the ¢ 
individuals in A’s acquaintance circle did not inform A on the (¢ + 1)th 
remove were all equal, that is, the associated events were all independent. 
Another way of saying this is that our knowledge that the first acquaint- 
ance did not inform A did not affect our assumption about the state of the 
second acquaintance, etc. If we drop this assumption of independence, the 
compound probability that none of A’s g acquaintances informed A can no 
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longer be represented by (7). We must write instead of the gth power a 
g-fold product 


SF 1 -P.W ml, (15) 
k=0 


where the P,(#) are conditional probabilities to be determined below. 
Using expression (15) instead of (7) in equation (10), and solving for 
a(t), defined by (4), we now have 


q—1 
a(t) = > log [1 —P; (4) m]. (16) 
k=0 


—1 
Pb) 
For m <1 (and a forteriori P,(t)m < 1), the logarithm in (16) can be 
well approximated by — P;(#)m, and we have the simplified form of a(¢), 
namely, 


m <A 
a (t) 10) al (17) 


For the special case of the completely mixed population all the P; are 

equal, and 

a(t) =qmaa, (18) 
as, of course, should be the case. It now remains to determine the P, on 
the basis of certain assumptions which we wish to make about our socio- 
structural bias. 

The conditional probability P,. Let us compute the following conditional 
probability: given that an arbitrary selected acquaintance of A did not in- 
form him on the (¢ + 1)th remove, what is the probability that this ac- 
quaintance was not a new knower on the ‘th remove? 

We apply Bayes’ Rule 


b(i|k) =P ELH? He) 


DS (El Hi) pA, a 


which gives the probability of a hypothesis H;, given the occurrence of the 
event E, in terms of the probabilities of the event E, given each of the 
possible hypotheses H;, and the probabilities of these hypotheses. In our 
case we have two hypotheses, namely, 

H,: he is a new knower on the /th remove, and 

H;: he is not a new knower on the ‘th remove. 
Therefore 


P(Hi) =P(t); bp (Ha) =1—P (i). (20) 
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Our event E stands for “he did not inform A.” Hence 


1\2 
E = _— = — _ = 
p(E|Hy) =(1-2) =t—m; pl) =1. (21) 
The probability we seek is 
| _ 1—P) 
p (H2|E) =5—; 4100 E (22) 


whence the complementary probability, namely, that the acquaintance of 
A who did not inform him is a new knower on the ‘th remove, is 


P(t) (1—m) 


P, (2) =1-— )(H,2|E£) + bP ne 


(23) 

The assumption of transitivity. Let us now think of our population of V 
individuals as composed of NV subsets or neighborhoods, each composed of 
a “representative” member of the population and his g acquaintances. Of 
course, the neighborhoods are overlapping. The over-all density of new 
knowers in the population on the ‘th remove is by definition P(é). How- 
ever, if we select a subset of the neighborhoods according to some criterion, 
we can expect that the density of new knowers in that subset will be dif- 
ferent from P(¢), depending on the criterion we use in the selection of the 
subset. Let the criterion now be the following: we select all those neighbor- 
hoods in which the representative individuals were not informed on the 
(¢ + 1)th remove by an arbitrarily selected individual in the correspond- 
ing neighborhoods. According to our argument above, the fraction of new 
knowers among those non-informing individuals is not P(é) but slightly 
less, namely, P;(#), given by equation (23). 

The “‘assumption of transitivity” which we now make is the following: 
P(t) can be taken as the density of new knowers in the subset of neighbor- 
hoods characterized by the criterion above (i.e., those neighborhoods in 
which none of the arbitrarily selected first acquaintances of the represent- 
ative individuals informed the representative individuals on the ¢th re- 
move). This is tantamount to making the following two assumptions: 

1. The density of new knowers in a subset of neighborhoods can be 
taken as the density of new knowers in a set of individuals, each ar- 
bitrarily selected from each neighborhood. 

2. If instead of choosing a given set of individuals, to define a subset of 
neighborhoods, we substitute for each of the individuals an acquaintance 
of his, we get a subset of neighborhoods practically identical with the 
original one. 

The validity of the first assumption depends on the size of the subset. 
If the number of individuals in it is sufficiently large, then the density of 
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new knowers in a sample (1/¢)th the size of it can be taken to be the den- 
sity of new knowers in the whole subset. The second assumption implies 
that the neighborhoods (acquaintance circles) of two individuals who are 
acquainted are very strongly interlocking. If the population were divided 
into mutually exclusive cliques, within which all the individuals were ac- 
quainted with each other, then the second assumption would be satisfied 
exactly, because in that case any individual in a clique would represent 
the same clique. Our model, however, differs from the exclusive clique 
model, as will appear in the discussion below. 

We see that our assumptions imply that on the ‘th remove the density 
of new knowers in neighborhoods, characterized by the non-telling by the 
first arbitrarily selected acquaintances to the representative individual, is 
P,(). By similar reasoning, the density of new knowers in neighborhoods 
characterized by the non-telling of the first two arbitrarily selected ac- 
quaintances will be 


_ P(t) A—m) 
P2(t) = i=in (24) 
and, in general, 
p, 0) Sf @O tems (25) 


1 —Py—-1\(t) m 


Hence the probability that the representative individual is not informed 
by any of his acquaintances is 


q—1 
I] [1 —P, (¢) m], (26) 


and our quantities P;, are identical with the P,(é) of expression (15). 
- We can now show by induction that 
fF (i). $* 
1—P(t) +P (ft) s*’ 
ats s = 1 — m. For suppose that (27) holds for &. Then, according to 
25), 


P,(t) = (27) 


ns P ()es® P(t) s® =I 
Pete eee eae a. 
*) 1-P() FPO z[1 1—P) +P 1 s) | 
(28) 
P (é) gkt1 


~T—PW +P se 
which establishes the induction. 
To obtain the expression a(#) as given by (17), we seek 


ql 
ay Belay 


k=0 
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Since m <1, and s = 1 — m, s is nearly unity, and its successive powers 
differ slightly from each other. We can therefore approximate our sum by 
the integral 


a Pt) s* <teal 

LPH EPH HEM ys BAP = 991. 9) 

Recalling that s = 1— m ~1 — a/q, we can further approximate the 
right side of (29) by 


1 
Se eel SP hed re) Ni, (30) 
thus obtaining the desired expression, namely, 


a (t) CONE S09 8 Se at ie (a1) 


eee. 
P(t) 
Equation (31) holds for ¢ = 1. For ¢ = 0, we must take a(é) as a in equa- 
tion (1), because the initial knowers, being selected entirely at random 
from the population, determine a set of entirely randomly selected neigh- 
borhoods, so that the very first step of the process must be supposed to be 
governed by randomness. We thus have the apparent axon density given 
by two equations, namely, 


120) 


Ah Kee any ee ee eB 


(32) 


a (t) ==—~ log [1 —P (#) (1-—e)], i> 0. 


P a 
If P(O) is sufficiently small and a not too large, P(1) will also be a small 
fraction. In that case a(1) will be approximately equal to 1 — e~, or, for 
a > 2, quite close to unity, where it will remain during the remainder of 
the process, since the successive P’s must then become even smaller. In 
part III of this series, some indication will be given on the closeness of 
approximation by (31) to (17). 

The large drop from a(0) to a(1) is thus accounted for by the assump- 
tion of transitivity, but not the subsequent rise of a(#). A modification of 
the assumptions toward that end will be undertaken in part II. 

The spread of state by removes. Equation (4) gives the ordinal time course 
of the process if an explicit expression is substituted for a(¢). In our case, 
this is given z (31). Thus 

1—<«x(t) 
a0) log Lhe Bibl seest) Lie SUI 8 [aT 
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Recalling that P(t-+ 1) = x(t + 1) — a(t), and rearranging, we obtain 
P(t+1) =([P()] [1-«()] [1—eJ, #5 0% (33) 


The “physical meaning” of equation (33) is the following. Suppose all the 
new knowers were “bunched” in clusters and normally talked only among 
themselves. If, however, a non-knower happened to get into one of these 
clusters, he could also be contacted by one of its members with the same 
probability as any other individual in it. Let the average size of these 
clusters be r individuals. Then, if once one gets into a cluster, the prob- 
ability of receiving the information will be 


1-(1-2)'at—e>, (3.4) 
i 


that is, this probability will be practically independent of the size of the 
cluster if the cluster is sufficiently large, so that the approximation (34) 
holds. On the other hand, the probability of getting inside the cluster 
equals the probability of meeting a new knower on the /th remove, that is, 
P(é). Finally, 1 — x(é) is the probability of being a non-knower at ¢. This 
is the meaning of the three factors in (33). 

Discussion. Although one of our assumptions above is equivalent to as- 
suming mutually exclusive cliques, our model is not equivalent to the 
closed clique model. In the latter, the probabilities of becoming a knower 
cannot be taken independently for each remove. If, for example, an indi- 
vidual is not a knower by the ‘th remove, the knowledge of this fact should 
be reflected in assumptions about the state of affairs in his clique. For ex- 
ample, the probability that it had no zmitial knowers is increased thereby. 
This dependence of the distribution of knower densities on ¢ does not ap- 
pear in our equations. 

The kind of situation reflected in equation (33) is, where there is a lim- 
ited mixing in the population in the sense that the “‘acquaintance circle” 
does not stay fixed but changes with /, rather as if individuals moved 
through the population, and their “circles” were defined as simply their 
geographical vicinities. The mixing is limited, however, in that there is a 
certain “inertia,” i.e., new knowers tend to talk almost exclusively to new 
knowers. This happens because a group of new knowers, told by individ- 
uals close together are also close together. It is as if the contacts “lin- 
gered” for a while, so that groups of new knowers being in the immediate 
vicinity of their informants, continued in the vicinity of each other and 
passed their information only to those who happened to wander close by. 

The conclusion that under these conditions the apparent axon density 
becomes approximately unity indicates that the assumption of transitivity 
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is too strong. In practically all the information diffusion experiments per- 
formed, a(0) dropped from various values (ranging from 3 to 12) to a 
value of a(1) ranging from 1 to 2, with only very occasional values greater 
than 2 or smaller than 1 recorded. It appears, therefore, that a weaker re- 
striction on contacts than the one implied by our assumptions should be 
chosen. 

On the other hand, a situation in which a(é) drops from an arbitrary ini- 
tial value to unity and remains there is easy to construct. Suppose all in- 
dividuals are arranged linearly on a line N units long and are denoted by 
positive and negative integers. Suppose the individual denoted by 0 is the 
initial knower. Then P(0) = 1/N. Let each new knower now tell  indi- 
viduals on his right and on his left. Then, because of the complete over- 
lap, the individuals contacted by the extremes on either side will be the 
only new knowers. Hence P(é) will be 2”/N for all t > 0, and 


2nt+ 1 


x(t) = 7 (35) 


Although a, the actual axon density, is approximately 2, the apparent 
axon density is given by 
1 ant ED) 
N N 
«Oy 8 Tan 1) +17 
ee 
N 
for very large N and moderate ¢. In fact, under these conditions a(¢) will 
tend to rise slowly with t. 

If we construct a similar model in two dimensions, P(é) will grow linear- 
ly (as the perimeter of the area of knowers) and a(¢) will exhibit both an 
initial drop and a steady subsequent rise. 

It seems, therefore, that a model combining some features of random- 
ness and some of a spreading “wave” of knowers from a “focus of infec- 
tion” should be able to account for the facts. Some such ‘“‘mixed” models 
will be offered in part II. 

This investigation is part of the work done under Contract No. AF 
19(122)-161 between the U.S. Air Force Cambridge Research Laboratories 
and the University of Chicago. 


Soil (36) 
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The assumption of transitivity treated in part I is modified in various ways to describe an 
information-diffusion process, in which a certain amount of randomness of contact does occur. 
In one model a parameter is introduced which is indicative of a tendency to go beyond one’s 
immediate vicinity to spread the information as the vicinity becomes saturated with knowers. 
In another model the randomness appears in the assumption that new knowers are uniformly 
distributed among the knowers. Two of the equations thus derived, each with two free para- 
meters are in good agreement with experimental results. 


In a previous paper (Rapoport, 1953), hereafter called I, we tried to 
account for the drastic drop in the “apparent axon density” associated 
with the spread of information in a group from its initial value on the first 
remove. We found that the assumption of transitivity roughly accounted 
for the drop. In fact, for a moderate actual axon density a, the apparent 
axon density a(#) became approximately unity and remained at this value 
throughout the process. 

Experimental evidence, however (Washington Public Opinion Labora- 
tory, 1952), indicates that immediately after the initial drop from a(0) to 
a(1), the apparent axon density begins to rise and continues to increase 
throughout the process. In the present paper, we will suggest a number of 
models in which this growth can be accounted for by the introduction of a 
single new parameter. This parameter will be a sort of measure of the 
degree of socio-structural bias in the process. When the value of the 
parameter is the largest, the process will be that associated with the as- 
sumption of transitivity described in I; when the value is zero, the process 
will be characterized by complete randomness, as described in earlier pa- 
pers (Rapoport, 1951; Solomonoff and Rapoport, 1951; Landau, 1952; 
Landau and Rapoport, 1953; Rapoport and Rebhun, 1953). In another 
model “partial transitivity” will be introduced in another way without the 
introduction of a new parameter. The parameter will then have a different 
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meaning, namely, it will represent the “degree of participation” in the 
process. 

Under the assumption of complete transitivity, the conditional proba- 
bility that an individual’s kth acquaintance is a new knower on the ith 
remove, given that none of the k acquaintances gave the information on 
the (t + 1)th remove, is taken to be the a priori probability of being a new 
knower for the (k + 1)th acquaintance. This transitivity of probabilities 
is carried through the g acquaintances. 

Let us now assume that this transitivity is carried through only a frac- 
tion 6 of the g acquaintances, while the probabilities of being a new knower 
for the remaining fraction 1 — 6 is equal to the over-all probability P(¢) of 
the population. This is tantamount to assuming that out of all one’s con- 
tacts 6g are chosen within one’s acquaintance circle, where the assumption 
of transitivity holds, and (1 — 6)g are chosen from the population at 
large. 

Following the same reasoning as in I, we now have* 


P S® 
i ae pose 2p 57 for k<6q (1) 
P.=P, for k>6q. 
Then 
cae 69 = Pskdk 
Zhi if aera aniae nes 
ae 9 ad kee 1 
aiegeen wea Ried 8 ee BS PY 
We thus get 
—1 
a t— — ea — 
a (t) Pw 8! BLL 6 SI 0; mg (3) 
~1i—e “84+ (1-4)a, for’? 1% 


We can see that in this case a(f) is still virtually independent of ¢. How- 
ever, we can now fix a(¢) at any value between 1 — e~@ (for 6 = 1) toa 
(for @ = 0). In the experiment referred to in I (Washington Public Opinion 
Laboratory, 1952), the value of a(1) was observed to be 1.47. This value 
is obtained if we set @ = .95, that is, assume about 5% randomness in the 
process. Of course, under this assumption a() will remain at approximate- 
ly the same value for all ¢ > 1. 


If, however, we assume that @ decreases with f, the growth of a(é) can 


* Throughout these papers whenever variables dependent on ¢ appear without their argu- 
ments, the argument is understood to be ¢. Thus Pk = P(t), 7 = a(#)s ete; 


i yt ah ti Nia 


= 
= 
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obviously be accounted for. In fact, a linear decrease of @ with t, say, 
6 = .975 — .025t, gives excellent agreement with the experimentally ob- 
served values of a(t) for the first several removes. Such an assumption of 
dependence of @ upon ¢ would be strictly an ad hoc hypothesis and would 
not add to the understanding of the process. It is possible, however, to 
give a rationale for the dependence of @ upon x, the saturation of our popu- 
lation with knowers. It is reasonable to suppose that as the population 
becomes saturated with knowers, a tendency to seek out non-knowers 
makes itself felt, especially under the contest conditions of the experi- 
ment.* Hence if we suppose that this active search is reflected in going 
outside one’s acquaintance circle more and more frequently as one fails to 
find non-knowers within one’s acquaintance circle, 9 should decrease as « 
increases. As an approximation, let 


0 = eB , (4) 


Thus when knowers are very scarce in the population, @ has the maximum 
value unity, and we are dealing with the completely transitive situation 


TABLE I 
t 0 1 2 3 4 5 6 7 8 9 
P(t) obs. WO30ra 191) (1901) ASG e107. 4.0636 4035554018) «0115004 
P(t) calc. PO3G (191) 1190), 102 5 dO O08 5053.) 014 ae < 006 — 004 
x(t) obs. OSGmiG22LMGZI2)" SSO8MELGIS= SIS8iw iS IA 91 | 780205806 
x(t) calc. 20305 (228) (411) eo iar 080. ova (90 S042 .810. 7 28t4 


Comparison of observed values of P(#) and x(t) with those calculated from equation (8). The 
values for t = 0 are directly given. The values for ¢ = 1 are used to calculate a = 7.2, those 
for t = 2 are used to calculate 8 = 0.22. These values of the parameters determine all the 
succeeding values. Note that the agreement holds not only for the cumulated fraction of 
knowers x(t) but also for each fraction of new knowers P(?). 


treated in I. The parameter 8 measures the tendency to seek random con- 
tacts when one’s vicinity becomes saturated with knowers. 

Table I gives the values of x and P in the experiment referred to. The 
value of x(1) has already been used to determine the actual axon density 
a, which, we recall, was fixed at 7.2 (cf. I). We now use the value of x(2) 
to determine 8. Since to give a(1) the value 1.47 (see Table I in I) we 
must have the initial value of 6 equal to .95, it follows that 


B= 5-557 1086-99) ~0.23. (3) 


* Prizes were offered for both becoming a knower as soon as possible and for giving as much 
information as possible to non-knowers. 
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N.B. In this range 8 is extremely sensitive to 6, so that the value 0.23 
should not be taken too seriously. As it turns out, 8 = 0.22 gives the best 
fit. 
Let us now find the ordinal time process under our present assumptions. 
Substituting (4) into (3) and setting (as has been shown in I) 
1 1— x (2) 
ol) Ray! ae 


we obtain 


log [1—P (1 — e-2¢**)] + (1— e**) a 
1 1—~-x (ft) 


PHS Tee 


oe 
oy (6) 


which leads after rearrangements to 
“(t+1) =1— (1—2) [1—P(1— e7*) | eae P (7) 


Equation (7) can be considerably simplified in the case under considera- 
tion if we recall that e-« * is always exceedingly small compared to 
unity, because of the large value of a, namely, 7.2, and because Bx < .23. 
Hence this term can be neglected, and we write 


o(t-+1) =1— (1—«) (1—P) e206 (8) 


As can be seen from Table I, equation (8) gives an excellent agreement 
between the observed and the calculated values for both x(é) and P(i). 
Equation (8) contains two free parameters, namely, a and 8, which were 
fixed by the second and third points of the observed curve. 

An alternative socio-structural model. Another approach is suggested by 
applying the assumption of transitivity to all knowers instead of the frac- 
tion of new knowers as in the preceding case. This means that the “‘co- 
hesive forces” operate on the knowers cumulatively. We can imagine such 
a situation if the mixing of the population takes place so slowly that the 
knowers tend to be in each other’s vicinity for a long time. 

Again we apply Bayes’ Rule [eq. (19) of I]. But now we have three 
hypotheses instead of two, namely, 

H,: he became a knower at #, 

Hz: he was a knower at ¢ — 1, and 

Hs: he is not a knower at tf. 


If the event E signifies as before “he did not tell at ¢+ 1,” we have 
1 a 
pE\H) =(1-2); pH) =p ELA, =1; 


(9) 
b(Hi) =P(t); (Hx) =x(¢—-1); p (Hs) =1—-« (2). 


eee ee 
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fHs|B) =1—x,@) = st * = 10) 
(1-2 )P+e@—1)+1-2 


Recalling that P(d) = «(¢) — «(¢ — 1), and that (1 — 1/g)" = 1 — m, we 
write 
1— x (2) x—mP 


x,(4) = ayer LAs) 


enh) 
rs) 1 — mP (t)’ mP 


The assumption of transitivity on x leads to the iteration formula 
xx (t) — mP; (t) 


p41 (t) = 1—mP;() (12) 
whence, by induction, 
k—1 
x—-1i+ al (1 — mP)j) 
x, = = (13) 


k—1 
BI (1—mP;) 
7=0 

In particular, 


mist |] (1 — mP;) 
Mts (pe et =z Be (14) 


q = 
[] (1 - Pp.) 
k=0 
Solving for the product term we get 


i (1—mP;) = a (15) 


Taking logarithms of both sides gives us 


<< 1—«x 
Slog (1 — mPx) = log : (16) 
k=0 i We Xq 

Remark. The case of complete randomness, where «, = « and P; = P 
for all & cannot be obtained as a special case of (16), because the iteration 
formula (12) does not apply, except in the trivial case where m = 0 or 
P = 0. For these trivial cases, both sides of (16) vanish identically. 

Equation (16) leads to the eae ja for a(t) 


== 07 (2) 
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and to the ordinal time course equation 
t—«“ti+1) =41—37 =a (18) 


To make use of equation (18), we must make some assumption concern- 
ing x,(é). We recall that x,(¢) is the density of knowers in the neighbor- 
hoods where the representative individual was not told the information by 
any of his g acquaintances. We can take two extreme cases: 

1. This set of neighborhoods is the same for all #, i.e., they are the 
neighborhoods where the process never started; 

2. These neighborhoods are indistinguishable from others at any pre- 
vious ¢, that is, the density of knowers in them up to ¢ — 1 inclusive is the 
over-all density «(¢ — 1). 

In the first case, x,(¢) = K, a constant. In the second case x,(#) = x(¢— 1), 
since if nobody passed the information at ¢ + 1, it follows that in those 
neighborhoods P(é) = 0, if in all other respects the g neighborhoods are 
indistinguishable from the others 

We are thus led to two equations, namely, 


(1 = 2)F 


1—x(i+1) = ixk (19) 
and 
i 2. tie 
1— zx (é+1) 5, RCT EY (20) 
We note that equation (20) may be written as 
1 — sii-- te) pi ie sas 
i—2(@) ~— fowlter oa 
a constant for ¢ > 1. We can furthermore evaluate K, by setting 
_1=#(1) , > P(t) 
On the other hand, the left side of (21) can be written 
_P(t+1) 
beeper net (23) 
Hence equation (20) leads to 
whe eD 2 
P(t+1) “Taos f SY) = KP). (24) 


Both K and P(1) depend on a, whence the dependence of the ordinal 
time process on a, which is not explicitly seen in equation (20). It is easy 
to see that the asymptotic value of « is unity for this case. Thus, although 
the spread is slower than the observed spread, eventually «(#) as calcu- 
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lated from (20) overtakes the observed values toward the end of the 
process 

To compute the ordinal time course for equation (19), we must deter- 
mine K. This can be done by using «(2) and x(1) in equation (18). We have 
(cf. Table I) 


0.41 = 1— (0.78)?(1 — K) 1=1— (0.61) (1—K)-!, K=-—0.03. 


In our case, negative «, is meaningless, but the value obtained is sufficient- 
ly close to zero to be taken as zero. This implies that there were never any 
knowers in the neighborhoods where there were no tellers. As can be seen 
from Table I, this assumption leads to a much more rapid spread of infor- 
mation than is observed. 


TABLE II 


t 0 1 2 3 4 5 6 7 8 


x(t) from eq. (19) .030 .221 .394 .633 865 .982 1.000 
x(t) obs. O30: o221 eet 568 se" 075 . 738 791 .802 
x(t) fromeq.(20) .030 .221 .374 .497 .596 .676 . 740 Bool 


Comparison of observed values of x(#) with those calculated from equation (19), assuming 
Xq(t) = O and those calculated from equation (20), assuming xq(t) = x(t — 1). 


On the other hand, equation (20) above leads to a slower spread than is 
observed, as one would expect if the event “‘no telling at ¢” could occur in 
an arbitrary neighborhood (cf. Table II). Actually, this event is most likely 
to occur in those neighborhoods where there are either too few knowers or 
too many old knowers (that is, where the saturation has already occurred). 
This situation could perhaps be described by a suitable parameter. We 
will, however, make a different assumption, as our final approach to the 
problem. 

The equidistribution of new knowers. We shall assume that in equation 
(12) we can take 

P,(t) = a (25) 
This is equivalent to assuming that in every neighborhood k the ratio of 
new knowers to the total number of knowers is the same. Strictly speak- 
ing, this assumption is justifiable only where the neighborhoods are indis- 
tinguishable, i.e., in the completely random case. But this is precisely why 
we are introducing this hypothesis—to weaken the assumption of transi- 
tivity. Thus the hypothesis is somewhat analogous to the one made above, 
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where the parameters 6 and 6 were introduced to provide partial random- 
ness in the process. Then equation (12) becomes 


eet) 


i (26) 


mnam@) 


which leads by induction [cf. the argument following (26) in I] to the 


formula 
x sk 


“eT aaa Sh 
where now s = 1 — mP/x. Substituting (25) into (27), we have 
D ss 
en ee a 


Following the reasoning above, we again approximate =P; by an integral 
and obtain 


iy (ee —_ dk = Fogg boelt— xb s9)]. (29) 


Approximating log s by —mP/x, a(t) by its expression (17) in I, and 
S$, by e~*’’ , we finally obtain for the apparent axon density 


a(t) = iy log {i—-s Ge rsT ; (30) 


and for the ordinal time course 
1—«2i-1l) = CL 2) [a Cl =e eFi) ] (31) 


which can also be written as 
Pti-+- 1) see (i) bets ke eels) (32) 


Note the relation between equations (32) and (33) of I. In the previous 
case, new knowers arose when a non-knower met a group of new knowers 
(who “traveled in bunches”) and was not missed by all of them in the 
exchange of information among them. Equation (32), on the other hand, 
implies that a new knower arises when a non-knower meets a group of 
knowers (who in this case “travel in bunches”) and is not missed by all the 
tellers (new knowers) in the bunch, as is indicated by the last factor in 
(32). Thus equation (32) is the natural analogue of (33) in I, where now 


iss Li . 
the “cohesive force” operates on all knowers instead of only on new 
knowers. 
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Obviously equation (32) should give a considerably faster spread than 
(33) of I, since it is much easier to meet an arbitrary knower (who drags 
new knowers in his retinue) than a new knower. This is, in fact, the case: 
the spread calculated from equation (32) is faster than the observed 
spread, whereas that calculated from equation (33) of I is slower than the 
observed spread. 

We still have not introduced our second parameter, however. This we 
will now do by assuming that a fixed fraction of the population do not 
participate in the process at all. That is to say, contacts are wasted on 
them. They are not registered as knowers, nor do they pass the informa- 
tion on. Under this assumption, the basic equation of the random net is 
modified to 


fe ete ge). Pee ORK], (33) 
where f is the fraction of the population participating in the process. Equa- 
tion (4) of I should then read 
f—«(t) 


a (ft) =5a PCR (34) 
and, accordingly, equation (26) becomes 
Gh es *) [ho (1 — et) 1, (35) 
or 
PU 1 alsa) et) (36) 


Obviously the parameter f can be used also in correcting the ordinal 
time courses derived from equations (19) and (20). They now become re- 
spectively 


P(it+1) =x(f-—<«) Cave) 
and ae 
ee AG jp 

Pita lene a oan ek (38) 


It is obvious, however, that (38) will give a spread even slower than that 
obtained from (20). Equation (37) fits not badly with the value of f used in 
fitting equation (36). Equation (37) is, of course, attractive, because of its 
extreme simplicity. However, in its generalized form it necessitates three 
free parameters, a, f, and K, and this makes the goodness of fit less sig- 
nificant than in the case of equations (8) and (36) where only two free 
parameters are used. 

To estimate f for equation (36), we first attempted to obtain its value 
by fitting the value of «(2). However, the value thus obtained was greater 
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than unity, which is, of course, absurd. Hence we let f = 1 for the first 
remove, i.e., we supposed that every one participated on the first remove. 
This gave the value .393 for «(2), which is sufficiently near the observed 
value .412. We therefore estimated f from x(3) so as to get a calculated 
P(3) as closely as possible to its observed value of .156. We find that if we 
set f = 0.81, we get P(3) = .157. This value of f may have real signifi- 
cance inasmuch as just about 0.19 of the cards were missing at the end of 
the contest, which may just reflect the fraction of non-participants. The 


rao 1 2 3 4 5 6 7 8 9 


Ficure 1. Comparison of observed values of x(é) (circles) with those calculated from 
equations (37) (uppermost curve), (36) (middle curve), and (38). 


calculated values of x(t) from equations (36), (37), and (38) are compared 
with the observed values in Figure 1. 

Figures 2 and 3 show the comparison of the values of P(é) and x(t) ob- 
tained from equation (8) with those experimentally observed. 

Summary. We have deduced two models based on certain transitivity 
properties of the acquaintance relation, modified by certain parameters, 
which give good agreement with observed results in the ordinal time 
course of the spread of information. 

It must be stressed that the “acquaintance relation” as defined by this 
transitivity property turns out to be a rather fluid one. It should be inter- 


.20 


FIGURE 2. Comparison of observed increments of knowers (circles) with those calculated 
from equation (8). 


t 0 1 2 3 4 5 6 7 8 9 


Ficure 3. Comparison of observed fractions of knowers (circles) with those calculated 
from equation (8). 
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preted rather as defining a “vicinity” of an individual at a particular 
“time” rather than a definite set of individuals. Thus the assumption of 
transitivity is equivalent to the introduction of mixing into the population 
of which the infinitely rapid mixing implied by the completely random sit- 
uation is a special case. The mixing implied by our assumptions of socio- 
structural bias is rather analogous to an imperfect shuffling of cards, where 
after each shuffling the cards that had been together tend to stay together 
for some time longer. 

It would be interesting to see to what extent either of our models is 
applicable in further experiments and whether experiments can be de- 
signed to test either model and possibly to decide which one is the better. 
A discussion of such experiments will be undertaken in part III. 

This investigation is part of the work done under Contract No. AF 
19(122)-161 between the U.S. Air Force Cambridge Research Laboratories 
and the University of Chicago. 
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Previous work (Macey, 1952) in the application of the one-factor theory to the heart is 
extended. The rate of production of the excitatory state is assumed to be linear. Two possible 
mechanisms are indicated whereby such a situation might arise. Assumptions are made re- 
garding the mode of action of the chemical mediators on the heart, and an equation is derived 
relating the heart rate to the frequency of nerve impulses traveling along the cardiac nerves. 
This result compares favorably with the experimental findings of A. Rosenblueth and F. A. 
Simeone (1934). Other experimental results are interpreted in terms of the theory. 


Many of the controlling mechanisms of the circulation of the blood in- 
volve reflexes which tend to maintain a relatively constant internal envi- 
ronment. The formulation of a quantitative theory of these processes re- 
quires a knowledge of the response or output of each element in a given 
reflex as a function of its input (frequency of nerve impulses or other 
stimuli). In this paper we are concerned with a particular effector element, 
the heart, and, in view of the foregoing, we shall seek an expression relat- 
ing the heart rate to the frequency of impulses impinging on the heart via 
its motor nerves. 

Previously (Macey, 1952; hereafter referred to as I), the one-factor 
theory was applied to the heart. It was assumed that the relation between 
the excitatory state ¢ of the heart and the intensity S(¢) of the stimulating 
agent could be written as 


a= mi eS cle be: (1) 


where F is some function of S(t), and # is a constant. Aside from this equa- 
tion and the postulate that excitation occurs whenever e reaches a con- 
stant threshold of magnitude /, the nature of e was left unspecified. It was 
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shown under fairly general conditions that certain experimental phe- 
nomena associated with the heart could be qualitatively accounted for if 
it was assumed that the effect of acetylcholine is to increase k, and the 
effect of adrenalin to decrease &. 

In this paper we shall assume that ¢ represents the concentration of a 
chemical substance ©. The attribution of the autorhythmicity of the heart 
to a chemical substance is an idea that appeared frequently in the earlier 
literature (Howell, 1905). For the results to be obtained here, this assump- 
tion requires neither a specification of the nature of the substance nor its 
site of action. Thus & or its precursors might conceivably be identified 
with any one of a number of organic extracts and metabolic products, or 
inorganic ions which have at one time or another been proposed as the 
source of the heart’s spontaneous rhythmicity (Wiggers, 1934; Howell, 
1936; Paes, 1949). 

Whatever the substance may be, it seems reasonable to assume that it 
arises as a result of a chemical reaction. In the simplest case, it would be 
formed directly from a precursor P which is present in great abundance. 
Such a situation would not suffice because, due to the large value of & (as 
shown in I), e would approach its asymptotic value too rapidly to account 
for the interval between beats. If the concentration of P is not excessive, 
we must take its source into account. With this added complication, the 
simplest case would have P formed at a constant rate. This would occur, 
for example, if P were the product of a saturated enzyme-substrate com- 
plex. 

Even if the above objection is now overcome, we are still faced with the 
problem of reducing ¢ each time the heart beats, in such a way as to make 
€ periodic. One way of satisfying this requirement is by assuming that any 
P that has accumulated during each cycle along with the P that is formed 
during excitation is immediately destroyed at that time. For then, during 
the finite period of excitation, no € would be produced, and the large 
value of k indicates that € would be destroyed very rapidly. Thus at the 
end of each excitation the heart would start its resting metabolism with 
negligible concentrations of P and €. With no other interference, it fol- 
lows that ¢ will be periodic. The change in the metabolism of P during 
each excitation is possibly caused by one or more of the many sudden 
changes which accompany excitation, for example, changes in local ion 
concentrations. 

Since € is periodic, we may define our time scale by setting e(z) = 0, 
i.e., we set ¢ back to zero at the very first moment of each excitation. Let 
the length of time during which there is no production of © be denoted by 


EXCITATION OF THE HEART 549 


6. With these conventions, the system is described by the following equa- 
tions: 


[P] =0 ae 
ol = 9— «IP for #25; | se 
«=0 tori 6 : 

oe = x [P] — he for t= 6, y 


where q is the constant rate of production of [P], « is a rate constant, and 
x|P] is equal to F[S(é)] of equation (1). 
The solution of (2) is 


[P] =2[1— en] (4) 


For simplicity, let « be small enough to approximate (4) by 
[P] = q(i—6) forbs. f= 7, (5) 


where T is the interval between beats. Then the solution of (3) becomes 
e=tq[u—0) fae) ], 5<t<T. (6) 


The value of & has been estimated at the order of 200 sec.— (I). J. R. 
Di Palma and A. V. Mascatello (1951) have analyzed the refractory period 
in the isolated heart muscle of the cat and found a period of the order of 
0.1 sec., during which the interjection of any arbitrarily large stimulus 
fails to elicit a response at any subsequent time. This suggests that during 
this period ¢ cannot be formed, and we will accordingly identify this dura- 
tion with 6. Using these orders of magnitude for k, we may approximate 
equation (6) by 


e=ZFg(t—6), §<t= Tad) 


The results obtained in this paper depend essentially on equation (7) 
and are independent of any particular mechanism involved in obtaining 
this expression. For this reason it seems worthwhile to consider an alterna- 
tive interpretation which leads directly to the formal equivaleut of (7). 

In a series of experiments on isolated heart muscle, E. Bozler (1942) 
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found that each conducted impulse was preceded by a period during which 
a non-propagated electrical potential gradually rose. The magnitude of the 
rise became smaller as the recording electrode was moved away from the 
pacemaker. These experiments led Bolzer to suggest that these electrical 
variations or “pre-potentials” are responsible for the normal rhythmicity 
of the heart. Later A. Rosenblueth and J. Garcia Ramés (1947) found 
these “‘pre-potentials” lacking in some of their spontaneously beating 
preparations, and present in some of their inactive preparations. They 
concluded that these “pre-potentials” were merely “epiphenomena”’ and 
were not responsible for the normal excitation of the heart. More recently, 
M. H. Draper and S. Weidmann (1951) have recorded cardiac resting and 
action potentials by means of an intracellular electrode. They found a 
slow drift in potential during diastole and have suggested that it is con- 
nected with the heart’s autorhythmicity. 

According to the theory, whenever the stimulating agent is an electric 
current, the magnitude of F[S(é)] in equation (1) is proportional to the 
current intensity (cf. I). The electrical drifts recorded by Draper and 
Weidmann can be very closely approximated by a linear function of time. 
Therefore, starting with the empirical fact that these drifts exist, together 
with the assumption that they constitute the stimulus which causes the 
heart to beat, we may substitute a linear function for F[.S(é)] and arrive 
immediately at the formal equivalents of (6) and (7). In this case, 6 might 
be interpreted as the time taken up by the action potential. 

Accepting either interpretation, equation (7) implies that the magni- 
tude of a minimal stimulus will depend on its time of occurrence within the 
cycle. That is, within any cycle, the earlier a stimulus is applied, the 
greater its minimal intensity. This has been actually observed in toads by 
Yamaguchi and Fujioka (Ishikawa, 1924). In contrast, J. C. Eccles and 
H. E. Hoff (1934), working with cats’ hearts, have found that the in- 
tensity of a minimal stimulus does not change during the latter portion of 
the cycle. They concluded that despite the careful precautions taken, the 
crucial region somehow escaped the effects of the stimulus. If the pre- 
potential theory is correct, it seems apparent that the magnitude of a 
minimal stimulus should be smaller the later it is applied within the 
cycle. 

Let us now turn our attention to the influences of the chemical media- 
tors on the foregoing system. According to our model, acetylcholine in- 
creases k. Hence it is involved in the destruction of ©. The fact that acetyl- 
choline is effective in such minute doses (Clark, 1933) indicates that it 
probably takes part in some enzymatic process. Elaborating on the sug- 
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gestion of J. Welsh (1948, 1950) that acetylcholine acts as a co-enzyme, 
we shall assume that it is involved in the following manner: 


ky 

V+X° = X* (8) 
k 

X*+ €43X*-+ products , (9) 


where V represents acetylcholine (the transmitter of the vagus nerve), X° 
an inactive form of an enzyme which destroys &, X* the active form of 
X°, and ki, k-1, and ko are rate constants. 

Adrenalin acts by decreasing k, thereby antagonizing the action of 
acetylcholine. This could be accounted for by assuming that adrenalin, 
which we shall denote by S (sympathetic transmitter) competes with V for 
X°. However, we shall introduce an intermediate stép in this process in 
order to include the possibility of a situation similar to that implied by the 
sympathin theory of Cannon and Rosenblueth (Rosenblueth, 1950). 
Briefly, Cannon and Rosenblueth have been able to explain many of the 
diverse effects of adrenalin by assuming that it combines with one or both 
of two receptor substances. The products consist of a sympathin E£ and/or 
sympathin J, depending upon which receptor substance is involved. The 
excitatory effects of adrenalin are then attributed to sympathin £, the 
inhibitory effects to sympathin J. The possibility that S reacts directly 
with X° and that their product represents sympathin £ cannot be given 
credence because sympathin £ from other sources, e.g., the liver (Rosen- 
blueth, 1950), has the same effect as adrenalin itself. Consequently, we 
shall complete our system by the following reactions: 


ky 

S+Y=E (10) 
ky 

EAsX tes Ze; (11) 


where Y is a receptor substance, E an intermediate product (possibly 
sympathin £), Z the inactive product of Z and X°, and hy, k-2, ks, and 
k_3 are rate constant. 

There may be many more intermediates than are indicated in (8), (9), 
(10), and (11). These equations are intended to represent the over-all 
process and hence their mathematical consequences may be regarded as a 
possible approximation to a more complex system. The inclusion of the 
intermediate E has been suggested by the sympathin theory, but our re- 
sults are not necessarily dependent on the validity of this theory (for criti- 
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cal reviews of the sympathin theory, see Rosenblueth, 1950; Tainter and 
Luduena, 1951). 

One further complication must be introduced due to the fact that when- 
ever acetylcholine or adrenalin is applied locally to a tissue, its concentra- 
tion will diminish, because it tends to diffuse into the surrounding me- 
dium, and also because it is destroyed enzymatically. A. J. Clark and 
J. Rayentés (1938a, 1938b, 1939a, 1939b) studied this problem in the 
frog’s heart and came to the conclusion that the rate of removal of both 
acetylcholine and adrenalin is proportional to their respective concentra- 
tions. 

Let II, and II, denote the rates of production of acetylcholine and 
adrenalin respectively from sources other than those indicated by (8) and 
(10), and let a, and a, be constants of proportionality. Then as a result of 
(8), (10), and (11), we have 


ale »— 4, [V] tka [X*] —ki LV] (X9 (12) 
eee ind [X°) —k1[X*] (13) 
ols! =0,—«, [S] + &-2[E] — ke [Y) [5] (14) 
CIEL = bs {E] — bY] (SI (15) 
Ed aehy [E] [X°] — ks [Z] (16) 


with the additional material balance equations 
é,= [X°] + [X*] + [Z] (17) 
col Y] + [E] + [Z], (18) 


where ¢; and ¢ signify the total concentrations (free and combined) of X° 
and Y respectively. 

We shall consider a particular case of the steady state where c. >> ¢, > 
[X*]. This latter inequality will be satisfied whenever the total V (free and 
combined) is much less than c. In a steady state, expressions (12), (13), 
(14), (15), and (16) are all equal to zero. Eliminating [¥] from (15) and 
(18) yields 
fe (c= 121) [51 

=F 


[E] = ; 
1+ ree [S] 
—2 


(19) 
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Since ¢ >> c = [X°] + [X*] + [Z] > [Z], we have the approximation 


ES = (20) 


The inequality c > [X*] justifies the approximation of (17) by 
fet = 6 — [X*] (21) 


Solving for the steady state [X*] from (13), (16), (20), and (21), and let- 
ting K; = k;/k_;, 7 = 1, 2, 3; we have 


61K, [V] 


C2 K2Kz3 [S] (22) 
i+Kis1 | * 


[X*] = 


and from (22) and (9), the value of the decay constant k becomes 


al — koi Ki [V] 
k = ko exo = eoK2Ko(S} 4 (23) 
1+ Ke [S] 


From (12), (13), (14), and (15), we have 


he eee 19) aa Sy (24) 


a, as 


The sources of V and S are nerve impulses and/or direct applications 
(e.g., by injection into the blood stream). In addition there is evidence 
that vis liberated by the heart in the absence of nerve stimulation (Heub- 
ner and Runcan, 1939; Runcan, 1940; Bulbring and Burn, 1949). We 
shall represent this synthesis in the simplest manner by 


keg 
O=V, (25) 


where ky and &_4 are rate constants, and Q, the precursor of V, is assumed 
present in sufficient quantity to consider Q normally constant. 

Bulbring and Burn (loc. cit.) have investigated the action of acetyl- 
choline in relation to acetylcholine synthesis and, among other things, 
have suggested that: 1) “the normal rhythmic contractions of the heart 
depend on a synthesis of acetylcholine”; 2) “the synthesizing power is 
high in freshly excised auricle which beats well, and, in addition, the action 
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of acetylcholine which depresses this beat also depresses the synthesis in 
the fresh auricle”; 3) “the synthesis is low in the auricle which has ceased 
to contract, and the action of acetylcholine which starts the beat again 
augments the synthesis also”; and 4) “when the contractions have been 
restarted by acetylcholine, the synthesizing power approaches the value 
for fresh auricle.” 

To interpret these results in terms of the present theory, first note in 
equation (23) that k may be made arbitrarily small merely by making 
[V], the concentration of acetylcholine, sufficiently small. Hence the pres- 
ence of an adequate supply of acetylcholine is essential for the main- 
tenance of a normal rhythmic beat—for, if & is sufficiently small, the decay 
of € will be so slight that ¢ will always be greater than h. (Physically this 
might possibly correspond to a constantly depolarized membrane.) If 
acetylcholine is added to the normal system, & will increase [cf. equation 
(23)], the rate will be slowed (cf. I), and the synthesis, as represented by 
(25), will be retarded. 

The results pertaining to auricles which have ceased to beat need not 
concern us here. However, they are not necessarily inconsistent with the 
theory. For example, one plausible way to account for them is to assume 
that acetylcholine activates some enzyme that is essential in the formation 
of Q. If there is a limited amount of this enzyme, it would be completely 
activated by relatively low concentrations of acetylcholine. In hearts with 
low synthesizing power, virtually all of the synthesized acetylcholine 
would be dissipated by the choline esterase. Then upon addition of a gen- 
erous amount of acetylcholine, the enzyme would be more activated re- 
sulting in an enhanced synthesis which in turn would tend to maintain the 
activated enzyme, etc. In this case, & may increase to a level high enough 
to reduce e below /, and the beat will be resumed. These effects do not 
take place in the normal heart, because its synthesis is sufficiently high to 
maintain the enzyme in a completely activated state. 

Let v, and », represent the frequency of nerve impulses traveling down 
the vagus and sympathetic nerves respectively, and let £, and &, represent 
the rate of administration of V and S respectively from other sources. 
Then, using (25), and assuming that each nerve impulse produces a con- 
stant amount of transmitter, so that the rate of production of the trans- 


mitter can be considered approximately proportional to their frequency, 
we may write 


Il, = mMyry + Eyt+ ky [Q] ~~ h~1V] (26) 
Il,= m,v.+ & , (27) 
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where m, and m, are constants. Equations (26) and (27) together with (24) 
yield 


1 
(V] = (maint Eat Be (OI) (28) 
[S] =— (m+ £) (29) 


Thus the steady state concentration of a transmitter is proportional to the 
frequency of impulses responsible for its liberation. This result is in agree- 
ment with the experiments of A, Rosenblueth and R. S. Morison (1934) on 
the cat’s nicitating membrane which led them to conclude that in the . 
steady state, the amounts of adrenalin-like mediator liberated is a linear 
function of the frequency of stimulation. 

At excitation, ¢ = T, e = h, and (7) yields 


peda i: (30) 
Kq 


If R denotes the heart rate (R = 1/T), then by combining (23), (28), (29), 
and (30) we arrive at the expression 


(vets) 
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The variables ¢, and ¢, have been included merely for the sake of com- 
pleteness. In what follows, we shall consider only those cases with ¢» = $s 
= 0; ie., where R is a function of v, and », alone. 


a2 
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There are five constants in equation (31). Two of these, B and 6, repre- 
sent quantities whose ranges are restricted. This leaves three parameters 
whose ranges are restricted only by the fact that they must be positive. 
It should be noted that 6 may be a function of [V] and [.], but in this first 
approximation it has been regarded as constant. 

The discussion will be simplified if we define o by the relation 


Vv 
= ET. F. oe 
i a1V,+ a2 (52) 
and note that o is a monotonic increasing function of v, and consequently 
of [.S]. Then (31) becomes 
ee Boo 


1 + 60 _ a3Vy : 
A. Samaan (1934) reported that small doses of adrenalin sensitize the 
heart to vagal impulses, whereas in larger doses it tends to reduce this 
sensitivity. The expression 
aR a3(B+o) 
av, (1+ 60+ 43»)? 


is a measure of sensitivity to v,. Differentiating (34) with respect to ¢ 
yields 


(33) 


(34) 


“o ek _ 43 (a3v,+ 1 — 26B— 6c) 
do\ wf — (sepa 


Since 6 is of the order of 0.1 sec., 26B is generally less than 1, so that, for 
very small o, 0/dc0(—0@R/d»,) is positive. If the asympotic value of 
is large enough, as o increases, 0/80(—0R/d»,) will pass through zero 
and then become negative. That is, under these conditions, the sensitivity 
to vagal impulses as a function of the concentration of adrenalin possesses 
a maximum which could possibly account for Samaan’s observation. 

If 6 and B are known, the parameters in (31) may be determined in the 
following manner: from (31) with vy, = 0, we have 


(35) 


B 
R(, 0) ‘a ery (36) 
or 
1 weed a3 
Ry, 0) =RtR-: (37) 


Thus the slope of the predicted straight line plot of 1/R(v», 0) against 
v, is equal to a3/B. To determine a; and a, we note that 


R(0,%) == (38) 
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which, when solved for o, yields 


R(0,».) —B 


ie BESY AGT I 7) 


Expression (39) enables us to calculate the values of corresponding to 
each R(0, v,) and hence to each »,. By rearranging (32), we obtain 


Vs 
Sei aa: ao 9 (40) 


and thus by plotting »./o against v., we obtain a, and a from the slope 
and intercept respectively. 
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FREQUENCY OF STIMULATION OF 
THE VAGUS NERVE (IMPULSES/SEC) 


FicureE 1 


This method has been applied to the data of Rosenblueth and Simeone 
(1934), and the ensuing results are illustrated in the figure. The value of 
B, the basal heart rate was obtained directly from the data. The value of 
5 was arbitrarily selected as 0.10 sec. Perhaps if greater care is exercised in 
selecting a value of 6, the experimental points can be more closely ap- 
proximated, but in view of the highly speculative basis of this paper, our 
purpose is merely to indicate a reasonable approximation. 
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From considerations other than those presented here, Rosenblueth 
(1932a, 1932b) arrives at expressions of the following form. 


Vv 
— '0) = ~—— (41) 
B RY; ) Hive fee 
Vv 
Ri, (42) 
EAE: Mast M4 


where ji, iz, #3, and ws are constants. All of Rosenblueth’s data in these 
publications are checked against expressions (41) and (42). An elementary 
computation will show that equation (31) formally includes (41) and (42). 
However, it should be mentioned that when dealing with R as a function 
of v, and y,, Rosenblueth and Simeone (Joc. cit.) have proposed an expres- 
sion which includes (41) and (42), but is not formally equivalent to (31). 


I am indebted to Drs. George Karreman and Clifford Patlak for criti- 
cally reading this paper, and particularly to Professor H. D. Landahl for 
his many helpful suggestions. 
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James J. Nickson, M.D., Eprror. Symposium on Radiobiology. The Basic Aspects of Radiation 

Effects on Living Systems. 1952. 465 pp. New York: John Wiley and Sons, Inc. 

This book is a result of a symposium held at Oberlin College in 1950. The purpose of the 
symposium was to give a broad view of fundamental concepts and to coordinate various ap- 
proaches to radiobiology in order that current problems might be more clearly formulated and 
that the present knowledge might be evaluated for what it is as well as for what it lacks. It 
was planned that the basic aspects of each topic should be presented together with the neces- 
sary background information. To do this without too much repetition or without too many 
omissions is no easy task. That the various authors were as successful as they were is a tribute 
to their cooperative effort as well to the effort of the editor and the members of the Editorial 
Committee, P. Morrison, M. Burton, E. S. Guzman Barron, G. Failla, and H. M. Patt. 

The twenty-three papers may be divided into five sections dealing with the following topics: 
(1) physical interaction of ionizing radiation with matter, (2) chemical changes which occur as 
the energy of radiation is degraded, (3) biochemical effects of ionizing radiation, (4) effects of 
ionizing radiation on cells, and (5) on the mammalian organism. 

The first topic is covered principally in the first four papers. These are entitled ‘‘Radiation 
in Living Matter: The Physical Process,” by P. Morrison, ‘‘Secondary Electrons: Average 
Energy Loss per Ionization,” by U. Fano, ‘‘Beams of High-Energy Particles,” by R. R. Wilson, 
and ‘‘Neutrons and their Special Effects: Recoil Effects,” by A. K. Solomon. These papers 
attempt to summarize some of the more important physics of the interaction of radiation with 
matter which one would do well to keep in mind while working with biological material. 

The second topic is dealt with primarily in the next five papers. In the paper ‘‘General 
Statements about Chemical Reactions Induced by Ionizing Radiation,” by R. Livingston, 
electronic excitation and ionization are discussed as the two paths of activation. 

In the paper ‘‘Chemical Reactions in the Gas Phase Connected with Ionization,” by 
M. Wallenstein, A. L. Wahrhaftig, H. Rosenstock, and H. Eyring, mass spectrometric data are 
presented together with an interpretation based on quantum mechanical and statistical calcula- 
tions. 

In the next paper ‘‘On the Primary Process in Radiation Chemistry and Biology,” by R. L. 
Platzman, emphasis is placed on primary processes other than simple excitation or ionization. 
Although these other processes are much less probable, they might conceivably play a sig- 
nificant role in certain cases. Special attention is paid to direct nuclear collisions, multiple 
ionizations, energy loss, when a positively charged particle captures an electron then loses it in 
a subsequent collision, and to Auger descriptions. In this latter case double ionization results 
when a K electron of C, N, or O is lost, this being followed by the simultaneous loss of two 
L electrons, one replacing the K electron, the other being ejected. 

In the paper by M. Burton, ‘‘Elementary Chemical Processes in Radiobiological Reac- 
tions,” is given, in addition to the elementary reactions of irradiated aqueous solutions, a dis- 
cussion of the effects on the target theory of free-radical diffusion, ionization transfer, chain 
reactions and local pH changes. In the paper ‘‘Influences of Details of Electronic Binding on 
Penetration of Energetic Charged Particles through Liquid Water,” by R. L. Platzman, is 
given a theory of the stopping-power for a monatomic gas as well as a resume of the data for 
polyatomic gases and for water. 
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The section on biochemical effects is introduced by the paper ‘‘Some Aspects of the Bio- 
chemical Effects of Ionizing Radiations,” by W. M. Dale. This paper, together with the paper 
by E. S. Guzman Barron entitled ‘‘The Effect of Ionizing Radiations on Some Systems of 
Biological Importance,” summarizes some of the important concepts in effects of radiation on 
chemical constituents of cells-enzymes, coenzymes, nucleic acid, proteins. The importance ot 
thiol groups in enzymes as well as in protecting substances is emphasized. The effect of radia- 
tion on cell division and respiration is discussed. These topics are continued in the paper by 
G. Hevesy entitled ‘Ionizing Radiation and Cellular Metabolism,” in which chromosome rup- 
ture is discussed as well as the problem of the uptake of radioactive indicators. 

In the paper ‘‘Some Factors Influencing Cell Radiosensitivity by Acting at the Level of the 
Primary Biochemical Action,” by R. Latarjet, theoretical and experimental results are com- 
pared for the case of virus-infected cells where the multiplicity can be controlled. Similarly the 
differences in the effect of radiation on diploid and haploid cells can be accounted for. Also dis- 
cussed are the role of temperature, oxygen, and peroxidases. 

In the paper, ‘‘On the Localization of Radiation Effects in Molecules of Biological Impor- 
tance,” M. D. Kamen shows how by incorporating P® into phage particles, for example, one 
can obtain information concerning the sensitivity of the O-P bonds from the phage survival. 

The fourth section is introduced by the paper ‘‘Recent Evidence on the Mechanism of 
Chromosome Aberration Production by Ionizing Radiations,’”’ by N. H. Giles, Jr., in which 
particular attention is devoted to the interpretation of the influence of oxygen in increasing 
chromosomal interchanges. In the next paper, ‘‘Physical and Chemical Factors Modifying the 
Sensitivity of Cells to High-Energy and Ultraviolet Radiation,” A. Hollaneder stresses the 
marked effect of oxygen on the sensitivity of E. coli grown aerobically and anaerobically. The 
effect of water content and amino acids are also discussed. 

In the paper ‘‘Gene Mutations Caused by Radiation,” H. J. Muller discusses the relation 
between radiation and spontaneous mutations and the effect on a population of the increase in 
rate due to low level exposures. The implication of the results with ultraviolet light is discussed 
and the target hypothesis in its simplified form is criticized. 

The next two chapters, ‘‘Speculations on Cellular Actions of Radiation,” by R. E. Zirkle, 
and ‘‘The Dependence of Some Biological Effects on Radiation on the Rate of Energy Loss,” 
by C. A. Tobias, are of particular importance to those interested in mathematical theories of 
dose-effect relationships. In the former paper an outline of possible processes is given for the 
purpose of focusing attention on details of the mechanisms and suggesting lines of attack on the 
problem. A mathematical formulation is given and it is shown that rather complex dose-effect 
relationships are obtained. It would appear to the writer that the attempt to make the results 
too general has made it difficult to follow the argument. For example, if 2/s in equation (1) rep- 
resents the number of duplicate chromosomes, then one would expect that homologous gene 
loci would require the same number, m/r, of groups of decisive processes to result in inactiva- 
tion, so that all m/r’s would be equal, as assumed. But in this case equation (1) accounts for the 
mortality due to the defect of one gene only, death by other gene defects being ignored. Recent 
work (Zirkle, personal communication*) suggests that the number of loci may be quite large. 
In the second paper the problem is treated in a somewhat similar manner and the results are 
applied to the problem of change in ploidy with excellent results. The problem of the effect of 
rate of energy loss on the biological effect and the influence of ploidy are both discussed in 
terms of a diffusion model. 

The section on the effects of radiation on mammals is introduced by T. C. Evans’ paper 
entitled ‘“The Influence of Quantity and Quality of Radiation on Biologic Effect.” Some of the 
results on survival at various doses and the effect of the temporal pattern of the exposure to 
x-rays are summarized. The effects of x-rays and neutrons are compared. In the paper ‘‘Some 
Physiologic Factors Related to the Effects of Radiation in Mammals,” H. B. Jones discusses 


*To be published in Arch, Biochem. and Biophysics. 
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the effects of radiation on the blood, decreased uptake of chronic phosphate colloid by the liver, 
changes in rate of exchange of iron in blood and bone marrow, the depression of nucleic acid 
formation and the results of experiments on animals joined in parabiosis. In the paper ‘‘Mam- 
malian Radiation Genetics,” W. L. Russel discusses gene mutations and chromosome aberra- 
tions and estimates the genetic hazard of radiation. 

In the last paper ‘‘Analysis of Mammalian Radiation Injury and Lethality,” by A. M. 
Brues and G. A. Sacher, rates of mortality are analyzed for various species and for various 
exposure patterns. The analysis is carried out under the assumption that there is a single-dose 
injury or impulse function (which also defines a recovery function) whose time integral is 
additive with a normal constant aging process. This function is taken to be independent of age 
so that applications to wide ranges of ages should be made continuously. The misprint in sign 
in equation (1) is not likely to be misleading though the subscript ; which is missing from ¢ 
might cause some confusion, especially since in the short exposition one has to read carefully to 
ascertain the origin of the various ?’s. In this connection it may be pointed out that in (2) the 
fs should be the survival times from birth (or some extrapolated age) or else 8 times the age at 
beginning of exposure is neglected. This interesting method of approach is also applied to 
growth and to death rate from cancer. 

Although most of the material has previously been published the volume serves a very 
useful purpose especially since many references are not easily obtained. The discussions fol- 
lowing the papers also bring out many interesting observations and suggestions. A reader is 
likely to wish that more information were given on many specific topics, but on the whole a 
satisfactory balance has been made to fit the requirements of space. What may have been lost 
by the lack of organization inherent in such a joint effort has been more than compensated for 
by the presentation of the various points of view of the many authors. 


H. D. LANDAHL 
The University of Chicago 
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ERRATUM 


To the paper by G. Karreman entitled “Some Contributions to the 
Mathematical Biology of Blood Circulation. Reflections of Pressure 
Waves in the Arterial System” (Vol. 14, pp. 327-50). 

In the above paper equation (67) on p. 341 should read: 
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Equation (69) on the same page should read: 
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The calculations in the Appendix were actually based on these equations, 
as they appear in this erratum. 
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To the paper by D. M. MacKay and W. S. McCulloch entitled “The 
Limiting Information Capacity of a Neuronal Link” (Vol. 14, p. 131). 

In the above paper the statement that all 2 values of T, should ideally 
be used with equal frequency is true only if pulse-position modulation is 
used. With pulse-interval modulation it is advantageous to use smaller 
values of 7, more often. The improvement that can be effected does not, 
however, affect the order of magnitude of the results, and only rein- 


forces our qualitative conclusions. 
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Swelling stresses, 73 


Tissue, metabolism of, 93 


Viscous liquid, propagation of disturbance 
through, 395 
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